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of  tiie  Office  of  Naval  Research  (Code  VMl),  under  Contract 
No.  NOOOi 4-6 9- A-0 235-0001. 

The  principal  investigator  was  Professor  Albert  I.  King 
and  the  project  engineer  was  Dr.  Priyaranjan  Prasad.  The 
contract  was  monitored  by  Lieutenant  Commander  Kenneth  H. 
Dickerson,  MSC,  USN  of  ONR.  The  guidance  and  assistance 
provided  by  Cast.  C.  L.  Ewing,  MC ,  USN  and  Dr.  D.  J.  Thomas 
of  HAMEL ,  Detachment  No.  1,  Michoud  Station,  New  Orleans, 
Louisiana  are  gratefully  acknowledged. 


(iii) 


SUMMARY 


This  report  reviews  the  results  obtained  under  this 
contract  prioi  to  this  reporting  period  to  set  the  back¬ 
ground  Tor  the  qualitative  and  quantitative  documentation 
of  the  role  of  the  articular  facets  during  +  0,  acceleration . 

Experimental  evidence,  based  on  nearly  4  00  cadaver- 
runs  made  on  a  vertical  accelerator-,  is  presented  to  prove 
the  existence  of  dual  load  paths  in  the  human  spine  -  one 
through  the  intervertebral  disc  and  one  through  the  posterior’ 
structures  in  the  articular  facets.  Utilising  the  above 
£act ,  a  method  is  shewn  and  verified  experimental  ly  to 
increase  the  thr :  .  or  spinal  fracture  due  to  +  G 
acceleration  of  the  spine . 

A  78  degree -of- freedom  mathematical  model  and  its 
exper-i r.ohtal  verification  on  three  cadaveric  spines  is 


present  .-d  to  simulate  the  dynamic  response  of  the  human 
spine  during  +  G  acceleration.  The  dc-sigri  and  use  of  an 
intervertebral  load  cell  to  measure  axial  force  and  moment 
devel  oped  within  the  spine  during  +  G,  accel  -jration  is 
also  presented. 
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CHAPTER  I 


INTRODUCTION 


Several  hypotheses  have  been  advanced  to  explain  the 


anterior  wedging  fractures  of  vertebrae  during  ejections  of 


pilots  from  jet  aircraft.  Most  of  the  proposed  explanations 


consider  the  peak  acceleration  values ,  the  rate  of  onset  of 


acceleration,  and  consider  the  vertebral  column  as  a  single 


structure  or  a  series  of  single  structures  having  failure 


stress  levels  which  when  exceeded  result  in  vertebral  body 


fracture.  These  theories  consider  basically  a  symmetrical 


weight  distribution  in  the  mid-sagittal  plane  of  the  human 


body,  and  hence  do  not  take  into  account  the  rotation  of  the 


vertebral  column.  Vulcan  and  King  [20]  showed  that  the 


vertebral  column  was  subjected  to  significantly  high  bending 


strains  due  to  the  eccentricity  of  the  location  of  the  center 


of  gravity  of  head  and  torso,  with  respect  to  that  of  the 


vertebral  column,  resulting  in  a  forward  rotation  of  head 


and  torso  during  +  G  acceleration.  Due  to  this  bending  of 

z 


the  vertebral  column,  the  strains  developed  in  the  vertebral 


column  during  +  G  acceleration  were  the  highest  along  the 

z 


anterior  aspects  of  the  vertebral  bodies.  Orne  and  Liu  [15] 


studied  the  phenomenon  with  a  multi  degree-of-freedom 


mathematical  model  of  the  human  spine. 


Using  the  bending  theory  of  the  vertebral  column  it 


would  be  possible,  by  altering  the  initial  configuration  of 


J 


2 


the  spine  with  respect  to  the  acceleration  vector,  to  decrease 
the  bending  moments  on  the  anterior  aspects  of  the  vertebral 
bodies,  hence  increasing  the  threshold  of  fracture.  Hence, 
depending  on  whether  the  acceleration  vector  is  before  or 
behind  the  center  of  mass,  anterior  or  posterior  wedge  com¬ 
pression  fractures  respectively  would  occur  if  the  failure 
limit  in  anterior  or  posterior  bending  is  exceeded.  After 
studying  eighty  ejection  vertebral  fracture  cases,  Ewing  [  4  ] 
noted  that  almost  all  were  anterior  compression  fractures; 
only  one  was  a  posterior  compression  fracture.  Hence,  he 
advanced  the  following  hypothesis  [  5  ]: 

"Posterior  compression  of  the  vertebral  column  in  the 
thoraeo-lumbar  area  is  limited  by  the  articular  facets 
of  the  vertebrae  while  anterior  compression  is  not 
limited. " 

This  suggested  a  way  of  preventing  anterior  vertebral 
compression  fracture.  By  proper  hyperextension  of  the  spine 
a  configuration  between  adjacent  vertebrae  would  be  achieved 
that  would  allow  anterior  wedging  fracture  to  occur  only  if 
the  spinous  processes  connecting  the  articular  facets  to  the 
vertebral  bodies  are  fractured  or  if  the  ligaments  binding 
the  articular  facets  together  are  torn.  Also,  by  hyper¬ 
extending  the  spine  it  would  be  possible  to  induce  a  pre¬ 
tension  in  the  anterior  aspects  of  the  vertebral  bodies 
hence  increasing  the  bending  moment  and  axial  force  required 
to  cause  failure  of  the  anterior  surfaces.  Moreover,  by 
hyperextension  of  the  spine,  it  would  be  possible  to  reduce 


the  strains  developed  on  the  anterior  aspects  of  the  vertebral 
bodies  by  shifting  the  center  of  gravity  of  the  torso  with 
respect  to  the  acceleration  vector,  hence  allowing  lesser 
relative  angular  rotation  of  adjacent  vertebrae  as  compared 
to  a  normally  erect  spine. 

Experiments  were  carried  cut  to  verify  the  hypothesis 
and  its  corollaries.  The  following  chapters  describe  the 
experiments  carried  out. 


CHAPTER  II 


DETERMINATION'  OF  THE  OPTIMUM  SHAPE  AND 
OPTIMUM  LOCATION  OF  A  HYPEREXTENSION  DEVICE 

1 .  Introduction 

The  first  phase  of  the  study  was  to  obtain  the  best  loca¬ 
tion  and  shape  of  a  hyperextension  device  which  will  reduce  the 
peaK  strains  experienced  by  the  vertebrae  when  the  device  is 
not  used.  The  criteria  used  to  determine  the  optimum  location 
and  geometry  of  the  hyperextension  block  were  as  follows: 

(a)  The  block  and  its  location  should  cause  the  greatest 
overall  reduction  in  compressive  strain  levels  in 
the  vertebrae. 

(b)  Tensile  strains  in  any  of  the  vertebra  were  con¬ 
sidered  unacceptable,  because  tension  in  the 
vertebral  body  implies  severe  compressive  loading 
on  the  posterior  structures  and  may  cause  hyper- 
extension  fractures . 

c. .  Experimental  Procedure:: 

A  total  of  six  cadavers  were  instrumented  and  subjected 
to  eaudocephalad  acceleration.  Table  i  provides  a  list  of 
the  cadavers  used  and  the  number  of  runs  carried  out  on  each 
specimen.  Selection  of  the  cadavers  were  made  on  the  basis 
of  age,  cause  of  death,  and  body  weight.  Those  over  70  years 
of  age  were  rejected  and  if  the  cause  of  death  appeared  to 
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TABLE  1 


List 

of  Cadavers 

an  d 

N umb  e  r 

of  Runs  Carri 

ed  Out 

Run  No. 

Cadaver 

A r,e  at 

Cause  of 

No. 

Death 

Death 

* 

1  49*1 

59 

pulmonary  embolism 

1-28 

1580 

'17 

chronic  alcoholism 

29-148 

x582 

60 

pneumonia 

'49-50 

U59 

6  3 

o  arc inoma 

51-57 

1536 

68 

pulmonary  edema 

58-80 

66 

pneumonia 

*?  unnumbered  runs 
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have  a  considerable  effect  on  the  strength  of  the  vertebral 


column,  the  cadaver  was  not  used.  Roentgenographic  examination 


of  the  column  was  carried  out  on  suitable  specimens.  Lateral 


and  A-P  views  were  examined  for  pre-existing  fractures  or 


other  abnormalities.  Spines  v.ith  ahnormal  curvatures  and 


those  which  showed  a  high  degree  of  calcium  depletion  or 


were  arthritic  were  rejected. 


(a)  Geometry  of  the  hyperextension  devices: 


The  hyperextension  devices  consisted  of  wooden  blocks 


which  were  fastened  to  the  seat  back  with  two  countersunk 


carriage  bolts.  The  height  of  the  blocks  from  the  seat  pan 


was  adjustable  from  5"  -  20"  above  the  seat.  Four  of  the 


blocks  used  were  6"  wide  with  different  height  and  thickness. 


?nc  block  was  semi-circular,  3"  in  diameter  and  6"  wide, 


Table  2  shows  the  dimensions  of  the  blocks  used. 


’b)  Instrumentation  of  vertebral  bodies: 


The  procedure  for  instrumenting  the  cadavers  consisted  of 


vis '••©ration  of  the  abdominal  viscera  and  the  installation  of 


foi:-type  strain  gages  on  the  anterior  and  lateral  aspects 


of  the  vertebral  bodies.  in  every  cadaver,  the  three  vertebrae 


gat-'.-d  were  T12,  L2  and  Ut .  The  method  of  installation  of  the 


strain  gages  is  described  beiow. 


The  vertebral  bodies  to  be  gaged  were  cleaned  by  cutting 


'.rid  scraping  off  the  ligaments  around  it,  taking  care  that  the 


.surface  of  the  bones  was  not  damaged.  A  minimum  area  of  the 


none  was  exposed  to  ensure  some  continuity  of  the  ligaments 


on  the  bodies.  The  exposed  surfaces  of  the  bones  were  then 
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cleaned  and  dried  using  acetone  and  freon.  Foil  type  strain 
gages,  0.125”  gage  length,  were  then  installed  on  the  cleaned 
surfaces  using  Eastman  910  adhesive  and  catalyst  as  bonding 
agents.  The  installed  gages  were  then  tested  for  resistance 
to  ground.  A  minimum  of  500  megohms  of  resistance  to  ground 
ensured  long-lasting  and  noise-free  gages.  The  installed 
gages  were  then  coated  with  Gagekote  #3  to  ensure  insulation 
from  body  fluids.  The  leads  from  the  gages  were  carefully 
tied  down  by  sutures  to  a  disc  to  prevent  the  gages  from  being 
ripped  off  during  the  experiment. 

One  gage  was  mounted  on  the  anterior  surface  of  the 
vertebral  body  as  close  as  possible  to  the  mid-sagittal  plane. 

The  lateral  gages  were  mounted  in  pairs,  one  on  each  side  of 
the  body,  placed  as  closely  as  possible  in  the  same  coronal 
plane.  The  output  of  the  left  and  right  gages  was  summed  to 
eliminate  the  effects  of  lateral  bending.  All  gages  were 
installed  with  their  sensitive  axis  parallel  to  the  axis  of 
the  vertebral  body,  and  all  leads  were  pre-soldered  to  the 
gage  terminals. 

(e)  Description  of  accelerator  and  instrumentation: 

The  tests  were  carried  out  on  a  vertical  accelerator, 
housed  in  an  8-story  elevator  shaft  of  the  School  of  Medicine,  at 
Wayne  State  University.  The  sled  is  accelerated  over  a  stroke 
f'  8  feet  and  then  gradually  brought  to  rest  by  air  brakes  over 
some  30  to  40  feet.  The  acceleration  pulse  is  approximately 
trapezoidal  in  shape,  the  rate  of  onset  and  the  magnitude  of 
the  plateau  being  variable.  Details  of  the  accelerator  have 
been  described  by  Par  rick  [lb]. 
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I 


'.''he  restraint  system  consisted  of  an  automotive  lap  belt 
under  a  regular  U.S.A.F.  lap  belt  and  shoulder  harness  combina¬ 
tion  and  leg  straps.  The  wrists  were  tied  together  and  anchored 


t.  o 


the  seat  base  by  means  of  a  single  rope  going  through  an  eye 


bolt  in  the  seat  base  and  looping  around  the  U.S.A.F.  lap  belt. 
This  was  done  so  that  no  load  was  transmitted  on  the  arm  rest, 
and  to  prevent,  the  lan  belt  from  creeping  up  on  the  cadaver 
during  pretensioning  of  the  shoulder  harness.  The  head  was 
unrest  ra ined . 

The  electronic  instrumentation  consisted  of  12-channels 
of  bridge  balance  and  carrier  amplifier  units  (Heiland)  and 
a  24-channel  light-beam  recorder  (Visicorder) ,  The  sled  was 
equipped  with  a  50-g  strain  gage  accelerometer  (Statham  A6-50) 
and  the  shoulder  harness  load  was  monitored  by  a  1,000-lb. 
strap  load  cell. 

After  the  cadaver  had  been  placed  in  the  chair,  the  lead 
wires  from  the  strain  gages  were  connected  to  form  diagonally 
opposite  arms  of  a  4-arm  Wheatstone  bridge.  The  other  two 
arms  consisted  of  121-ohm  high-stability  precision  resistors. 
This  configuration  results  in  the  summing  of  the  output  of 
the  1  wo  gages.  The  anterior  gage  formed  a  two-arm  bridge 
with  a  121-ohm  resistor.  The  other  2  arms  were  provided  by 
the  bridge  balance  unit. 

(d)  Experimenr al  runs: 

The  lap  belt  was  always  snugly  tightened  and  restrained 
from  being  pulled  upwards  by  the  equivalent  of  an  inverted 
V-strap  described  in  the  previous  section.  The  shoulder 
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harness  was  pre-loaded  to  35  lb.  tension  before  each  run. 
Although  the  head  was  unrestrained,  its  initial  position  was 
kept  approximately  vertical  by  means  of  masking  tape  that 
broke  once  the  head  started  rotating. 

The  level  of  the  appiJod  acceleration  was  10  g  at  an 
cnsot  rate  of  approx  innately  300  g/sec.  The  low  acceleration 
level  was  used  to  avoid  vertebral  fracture  and  to  allow 
completion  of  a  long  series  of  runs  on  each  cadaver.  The 
acceleration  was  reduced  to  8  g  for  Cadaver  1^71  since  it 
was  the  oldest  among  the  ones  rested. 

In  order  to  compare  the  reduction  in  vertebral  strain 
resulting  from  the  use  of  the  hyperextension  device,  at  least 
two  runs  without  the  block  were  carried  out  on  each  cadaver. 
Cer.erally  the  first  and  last  runs  were  without  the  block.  For 
Cadaver1  1  -M 7 1 ,  a  no-block  run  was  made  after  each  block  was 
tested  to  monitor  more  closely  the  change  in  strain  gage  output 
caused  by  a  change  in  spinal  configuration. 

For  the  convenience  of  identification  of  gages,  the 
following  symbols  will  be  used  in  the  subsequent  sections  of 
this  paper : 

Tlv  prefix  A  denotes  a  gage  mounted  on  the  anterior 
surface  of  a  vertebra,  e.g.  ,  AT  L  2  is  the  gage  on  the  anterior 
surface  of  T12.  The  prefix  D  denotes  gages  mounted  on  the 
latera!  surface  of  the  vert ••tra.  Strain  values  of  DL2 
represent  the  average  output  of  the  two  gages  on  the  lateral 
surface  of  L2  on  either  side  of  AL2.  An  attempt  was  made  to 
install  these  gages  near  the  neutral  axis  of  the  vertebra  so 
that  their  output  indicated  mostly  axial  compression. 
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After  each  run,  a  lateral  roentgenogram  of  the  spine  was 
'.ak'-'P..  Experimentation  was  discontinued  if  fracture  occurred. 
In  some  cases,  the  x-rays  were  developed  after  the  completion 
cf  several  runs.  If  fracture  occurred  in  one  of  these  runs, 
the  data  from  all  pcst-frac!  *.••*  runs  were  discarded. 

Nummary  of  E  xr.o  r  j  n  t  a  1  >a  t  a 
A  'v;tal  of  ~)cl  runs  were  made  -  .ri  six  cadavers .  Of  i  hese , 
i x  runs  on  3  cadavers  constitute  the  bulk  of  the  data  reported 
Huns  Ad-5 7  on  Cadavers  1  A50  >  !53C  represent  incomplete 

•.•ries  due  to  fra p.  cm.  f  ver'.-brsj  nodi-::;.  The  raw  data 
f  typical  runs  are  shewn  ir.  Fi  gures  2.1  -  2.3- 

Tii«  vertebral,  strair.s  .'  jr !  ai  red  lr  those  runs  without 
b  locks  show  the  fa  mix '  a*  .-«  •.  ^  :  •  p- r"  •  i  by  Vulcan  and 

Ki.*n*  [;.!.*]•  An  example  of  these  con  trod  runs  is  given  in 
c’lrur--  .1.1.  B’i  .  1  '  s  i  type  cal  run  v/i  •  h  a  2-1/A"  thick 

block  pposite  ..1 .  Note  t h<  disappearance  of  the  second  peak 
an c  con;  .re  wlt.u  Figure  2.1  to  st  the  reduction  in  strains  on 
.  ga  ••••:;.  On  l.a  :  v-  r  . -V  .  the  ALA  rur  •  was  malfunctioning, 
Vf.CH  1 1 .5  record  i  .•>  mean.tr.  gl  ess .  Figure  2.3  is  an  example  of 
tension  developing  in  the  early  part  of  the  acceleration  pulse 
clock  No.  1  d  *  t  '  ■>)  war.  u.o  d  in  12  runs  at  various 

!  e  cat. !•*•', s.  In  Cadaver  IbcO  and  1  [  '62  ,  when  the  centerline  of 
t •  i •  •  I  i  .  be  war  placed  op:  »  the  f,3-I<A  disc  (Cadaver  1580), 

hi,  .2-53  disc  arid  he  ,  th<j  a:  t<Tior  surface  of  T12  (AT12)  was 
in  tension. 

By  increasing  tne  thickness  of  the  block  to  3"  (Block 
No.  2,  3  '  b  y  6),  it  was  observed  that  AT12  went  into  tension 
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when  the  block  was  opposite  the  L1-L2  disc,  L3-L4  disc  and 
1,5  in  Cadaver  1580,  and  opposite  Til  and  T12  in  Cadaver  1471. 
AL2  went  into  tension  when  the  block  was  2"  below  L5-  In  3 
out.  cl’  i  cases,  a  reduction  in  compressive  strain,  without 
any  vertebra  going  into  tons'  cu,  was  observed. 

There  were  12  runs  with  tne  semi-circular  block  (No.  3, 
o  in.  in  diameter  and  thick).  For  Cadaver  1580  ,  AT12  went 
int  .  tension  with  the  block  opposite  T3 2—1, 1  disc,  L1-L2  disc, 
arid  L3-L4  disc.  There  was  an  increase  in  compression  in  AT12 
and  17712  with  the  bj  o<  V  2  In.  below  1.5,  but  there  was  a  marked 
reduction  in  compressive  strain  in  AL4.  AL2  went  into  tension 
when  the  block  was  opposite  h5.  Hence  all  the  5  runs  with 
this  block  displayed  unuosir  .-tie  cnaracteristics .  The  results 
for  Cadaver  15&2  were  slightly  better.  Opposite  L3  and  L2-L3 
disc,  tension  de  lr  c^d  in  Al’l.:,  but  there  was  a  slight  re- 
■  uv.i.  r.  in  compi  ssiv-  strain  for  all  gages  with  the  block 
opposite  L2  and  T12-L1  disc.  In  Cadaver  1471,  there  was  an 
oV’-rai.  reduction  in  strain  when,  the  block  was  opposite  Til 
.  i  Ti<  ,  out  ATl.i  developed  tension  with  the  block  opposite  LI. 

The  results  from  Block  do.  4  (3  *  4  /  6),  were  also  un¬ 
set  i  s  factory  .  For  Cadaver  1580,  AT12  developed  tension  with 
the  slock  opposite  the  T 12-Li  disc,  L1-L2  disc  and  L3-L4  disc, 
and  il/Tx.--  showed  an  increase  in  compressive  strain  when  the 
o.-ox  '..as  2"  bei-w  Li.  Fur  'adiver  1582,  with  the  block 
crposite  Til,  an  increase  in  compressive  strain  for  all  gages 
was  observed.  However,  opposite  LI,  there  was  a  general  re~ 
lucb  l-.n  in  strain  on  all  gages. 


For  Cadaver  1471,  AT12 
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developed  tension  with  the  block  opposite  Til.  There  was  a 
slight  reduction  in  strain  when  it  was  opposite  the  T10-T11 
disc . 

The  results  from  Block  Mo.  5  (2.25  x  4  x  6)  were  en¬ 
couraging.  In  Cadaver  15t2,  ’d.>jro  was  an  overall  reduction 
in  compressive  strain  for  ail  gages  v/hen  the  block  was  placed 
between  T12  and  L2 .  There  was  a  minimal  increase  in  strain 
v/hen  it  was  located  opposite  Til.  Similarly,  for  Cadaver  1471, 
a  good  overall  reduction  in  strain  was  obtained  with  the  block 
opposite  Til,  T12  and  LI. 

4  •  Analysis  of  Data 

The  data  from  Cadavers  1580,  1532  and  1^4 7 1  were  read  off 
the  records,  converted  to  s-.  compared  with  the  control 

runs,  and  the  percentage  change  in  strain  computed  in  terms 
of  a  reduction  in  compression .  The  reduction  was  tabulated 
for  each  block  and  listed  in  Tables  3  through  7*  The  description 
of  the  results  in  the  previous  section  is  given  in  quantitative 
terms  ir.  these  tables.  In  order  to  evaluate  the  relative 
merits  of  these  blocks,  Table  8  was  prepared  using  the  runs 
in  which  no  tension  was  developed  and  no  increase  in  compression 
was  noted.  This  procedure  reduced  the  number  of  runs  from 
50  to  21 .  The  average  percentage  decrease  in  compressive 
strain  was  computed  for  the  anterior  and  lateral  gages 
separately  and  an  overall  average  was  then  obtained.  It 
can  be  seen  from  Table  8  that  the  2-1/4"  blocks  show  a  larger 
reduction  than  the  3"  ones,  and  in  the  former  group  the  best 
location  is  opposite  LI  for  every  cadaver.  An  exception  is 


TABLE 


#indicates  a  malfunctioning  strain  gage 


dicates  a  malf unct ionin 
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TABLE  8 

Analysis  of  Data  for  Runs  in  Which 
Tension  Was  Not  Developed 


Block 

Run 

Location  of 
Block 

Centerline 

A  ./or age 

Strain  Reduction  (%>) 

No. 

No. 

Ant.  Gages 

Lat .  Gages 

Overall 

1 

~> 

2"  below  L5 

51.00 

27 . 46 

39.23 

4 

L5 

69.53 

42.73 

56.13 

6 

L1-L2  disc 

44.26 

41 . 20 

42.73 

7 

T12-L1  disc 

29.80 

29.73 

29.76 

46 

T12-L1  disc 

57.26 

27.20 

42.23 

66 

Til 

36.60 

40. 30 

38.40 

67 

T12 

54.15 

49.60 

51.87 

6S 

LI 

65.25 

65.50 

65.  37 

2 

14 

T12-L1  disc 

4  3-2o 

43.73 

43.50 

3 

42 

T12-L1  disc 

34.80 

18.75 

26.80 

4  3 

L2 

39.9  3 

34.30 

37.12 

**»  t 

1  L 

T12 

38.60 

40.00 

39.30 

72 

Til 

3i.90 

25 . 00 

28.45 

'4 

34 

LI 

22.86 

27.70 

25.28 

75 

T10-TL1  disc 

16.45 

15.53 

15.49 

5 

36 

L2 

30.10 

30.06 

30.08 

37 

LI 

36.26 

27.30 

31.81 

38 

T12 

2  3.26 

11.64 

17.54 

60 

Til 

38.85 

34.96 

36.91 

61 

Top  of  T12 

52.80 

51 . 26 

52.03 

62 

LI 

59.90 

51.26 

55.58 

found  in  Cadaver  1580  using  Block  Wo.  1.  In  this  case,  a 
large  decrease  of  89%  in  AL2  in  Run  No.  9,  with  the  block 
opposite  L5,  indicated  that  the  LS  location  was  better  than 
Li.  However,  the  other  cae.avers  show  definitely  that  LI  is 
1  he  optimum  location.  Pur4  her  consideration  of  the  overall 
reduction  at  LI  for  Blocks  1  and  lj  shows  that  the  decrease 
is  50. 1%  fer  Block  No.  i  (Runs  6,  96,  and  63)  and  93-7?  for 
Block  No.  9  (Runs  37  and  6?) .  The  difference  is  quite  small 
considering  the  large  biological  variations  among  cadavers 
and  the  small  sample  sire.  The  consistent  indications  that 
Li  is  the  best  location  using  either  Block  No.  1  or  5  is  the 
major  conclusion  of  this  phase  of  the  study. 

0.  I: i  s 0 us slon  a n d  fori c  i  u.~  1  c  .s 

The  primary  objective  <  4>  selecting  the  shape  of  the 
hyperextension  device  and  it::  optimum  .1  cation  has  been 
achieved.  The  data  reported  from  the  three  cadavers  used 
..•voj  shat  the  block  should  be  about  2  in.  thick,  9  to  6  in. 

Led  gat  and  that  l»  ::  -out  Jr,o  should  be  opposite  LI. 

Before  the  start  of  this  series  of  experiments,  it  was 
expected  that  the  lateral  gages  would  show  an  increase  in 
.■■train  when  the  anterior  gaga.;  showed  a  decrease.  This  was 
based  on  the  hypothesis  of  a  more  evenly  distributed  loading 
f  the  vertebral  bodies  when  the  hyper^xtenslon  device  was 
used.  However,  in  ail  cases  the  lateral  gages,  too,  showed 
n  marked  decrease  in  strains.  This  brings  up  the  question 
>f  whether  the  vertical  force  is  being  transmitted  through 
the  block  or  through  the  posterior  structures  of  the  vertebrae. 
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The  segments  of  vertebrae  from  the  cadavers  tested  were  ex¬ 
cised  and  examined  radiologically  to  detect  fractures  of  any 
posterior  structures.  However,  none  was  found. 

The  development,  of  tension  on  the  anterior  surface  of 
the  vertebrae  is  another  problem  which  requires  further  in¬ 
vestigation.  An  examination  of  the  primary  and  secondary 
curvatures  of  the  vertebral  column,  shown  in  Figure  2A, 
reveals  that,  if  the  block  is  placed  between  L2  and  L‘5 ,  it 
forces  the  lumbar  vertebrae  forward  of  the  center  of  gravity 
of  the  torso  which  is  located  about  1/k  to  1/2"  in  front  of 
the  anterior  surface  of  T9  (see  Ref.  [19]).  It  can  be  seen 
from  the  data  that  the  tensile  stresses  occur  at  the  beginning 
of  the  acceleration  pulse,  before  the  upper  torso  has  had  a 
chance  to  begin  its  forward  rotation.  Also,  the  strap  load 
dropped  from  its  initial  pretension  (see  Figure  2.3)  during  the 
early  p,  rt  of  the  acceleration  pulse.  This  suggests  a  slight 
backward  rotation  of  the  upper  torso  during  that  period. 

Forward  flexion  of  the  torso  immediately  causes  the  gages  to 
go  into  compression.  By  placing  the  biock  opposite  LI,  the 
forward  displacement  of  t he?  segment  from  T12  to  I/I  is  just 
enough  to  reduce  the  compressive  strain  on  the  anterior 
surface  without  the  development  of  tension.  With  the  block 
opposite  the  lower  thoracic  vertebrae,  the  center  of  gravity 
of  the  upper  torso  is  pushed  farther  forward,  causing  a 
greater  eccentric  loading  of  the  spine,  which  results  in  an 
increase  in  compressive  strain.  The  question  of  whether  the 
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un  gaged  vertebrae  above  T12  would  go  into  tension  when  the 
block  is  at  its  optimum  location  for  the  T12  to  segment 
Is  still  unanswered.  However,  the  primary  thoracic  curve 
would  tend  to  prevent  tension  of  the  anterior  surfaces. 


23 


TABLE  9 

Variation  of  Breaking  Strength  with  Age 

Median  Breaking  Strength 

(kp  */cni2) 

10  7 
98 

76 

77 

43 

*  kp  is  a  unit  of  force, 

1  kP  =  1  kilogram  force 


Age 

(years ) 

20  -  30 
30  -  H 0 
40  -  50 
50  -  60 
over  60 


ATI* 


Cadaver  1^71,  Block 


M£ 


STRAP  LOA0 


» 
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CHAPTER  III 

EFFECTIVENESS  OF  THE  HYPEREXTENSION  DEVICE 


1 .  Introduction 

With  the  determination  of  an  optimum  geometry  and  the 
location  of  the  hyperextension  device,  it  was  decided  to 
document  the  effectiveness  of  the  device  used.  The  hyper- 
extension  block  that  produced  the  greatest  strain  reduction, 
described  in  Chapter  II  (2^-"  x  6"  x  V),  was  used  opposite 
the  best  location,  i.e.,  the  block  centerline  was  opposite  LI. 
A  total  of  12  cadavers  were  subjected  to  75  runs  using  three 
different  restraint  configurations  or  spinal  modes.  In  the 
hyperextended  mode,  in  which  the  block  was  used,  and  in  the 
erect  mode  with  no  block,  the  shoulder  strap  pretension  was 
set  at  30  lbs.  For  runs  In  the  flexed  mode,  the  harness  was 
loosely  attached  to  the  torso  and  the  cadaver  was  restrained 
essentially  by  the  lapbelt  only. 

The  verification  of  the  effectiveness  of  the  block  was 
carried  out  by  the  following  methods  of  data  analysis: 

(a)  An  analysis  of  the  acceleration  levels  at  fracture 
for  each  mode. 

(b)  A  study  of  the  mode  in  which  fracture  occurred  when 
two  or  more  modes  were  run  on  the  same  cadaver. 

(c)  A  comparison  of  strains  in  the  flexed,  erect,  and 
hyperextended  mode. 
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2 .  Experimental  Procedures 

Of  the  12  cadavers  used  in  this  study,  3  were  run  under 
all  three  modes,  5  were  tested  in  the  erect  and  hyperextended 
modes  and  4  were  run  in  the  hyperextended  mode  only.  Table  1 
contains  the  pertinent  data  on  the  cadavers  used.  Strain 
gages  were  installed  on  the  lower  thoracic  and  the  lumbar 
vertebrae  and  the  runs  were  carried  out  at  increasing  g-levels 
until  fracture  occurred.  If  more  than  one  spinal  mode  was  to 
be  run,  the  order  in  which  the  modes  were  tested  at  any  given 
g-level  was  picked  at  random  for  each  cadaver. 

The  procedure  for  the  preparation  of  the  cadavers,  the 
installation  of  the  gages,  and  the  roentgenographic  techniques 
employed  were  as  discussed  in  the  previous  chapter. 

In  Cadavers  1584  and  1615,  an  anterior  gage  and  a  pair 
of  lateral  gages  were  installed  on  T12,  L2 ,  and  L4.  There 
was  also  an  anterior  gage  on  LI  and  an  additional  pair  of 
gages  on  the  posterior  aspects  of  L2  in  Cadaver  16 15 .  For 
the  next  8  cadavers,  according  to  the  order  listed  in  Table  1, 
an  anterior  gage  was  installed  on  each  of  the  vertebra  from 
Til  to  L4.  It  was  felt  that  the  monitoring  of  the  strain 
along  the  anterior  aspects  of  all  6  vertebrae  was  more  im¬ 
portant  than  the  measurement  of  strain  along  the  lateral 
aspects  of  3  of  the  vertebrae.  In  the  last  two  cadavers, 

(Nos.  125  .and  127),  only  the  anterior  aspect  of  L4  was  gaged 
and  the  abdominal  cavity  was  not  eviscerated.  The  principal 
purpose  was  to  maintain  integrity  of  the  anterior  ligament 
along  the  major  portion  of  the  vertebra.’  column.  The 
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retention  of  the  abdominal  viscera  can  only  affect  the  fracture 
level  adversely  since  some  of  this  weight  is  borne  by  the  spine 
To  identify  the  location  of  the  gages,  the  prefix  A  is  used  for 
anterior  gages  and  D  for  lateral  gages.  Posterior  gages  on  the 
posterior  aspects  of  the  body  near  the  neural  arch  are  denoted 
by  the  prefix  DD. 

The  input  acceleration  was  a  ramp-shaped  pulse  with  a 
rate  of  onset  of  approximately  300  to  500  g-sec.  and  a  plateau 
varying  from  4.5  to  24. 5g.  The  duration  of  the  input  pulse 
varied  from  about  150-350  msec,  and  was  dependent  on  the 
acceleration  level,  since  the  total  stroke  length  of  the 
accelerating  piston  was  fixed  at  8  feet. 

3 .  Summary  of  Experimental  Data 

The  data  will  be  presented  under  three  headings: 

(a)  Fracture  g- level  and  spinal  mode  at  fracture 

(b)  Strain  histories  in  the  various  spinal  modes 

(c)  Examination  of  spinal  segments 

(a)  Fracture  g-level  and  spinal  mode  at  fracture: 

The  acceleration  levels  at  fracture  and  the  spinal  mode 
in  which  fracture  occurred  are  listed  in  Table  2.  The  number 
of  modes  tested  and  the  order  in  which  the  testing  was  carried 
out  and  the  fractured  vertebrae  are  also  given.  It  can  be 
seen  that  when  the  cadaver  was  subjected  to  all  three  spinal 
modes,  vertebral  fracture  always  occurred  in  the  flexed  mode. 
Similarly,  when  the  erect  and  hyperextended  modes  were  tested, 
fracture  occurred  in  the  erect  mode.  Cadavers  125  and  127 
were  not  eviscerated  and  strain  gages  were  installed  only  on 
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the  anterior  aspect  of  L4.  The  anterior  ligament  was  left 
intact  along  the  major  portion  of  the  vertebral  column  in  an 
effort  to  prevent  hyperextension  fractures  when  the  block  was 
used.  Such  fractures  occurred  in  Cadavers  062  and  095  which 
were  run  solely  in  the  hyperextended  mode  and  were  instrumented 
with  strain  gages  on  the  anterior  aspect  of  each  of  the 
vertebrae  from  Til  to  L4.  Cadaver  125  sustained  24. 5g  and 
its  spine  showed  no  fracture.  Due  to  an  accelerator  system 
malfunction,  a  higher  g-level  run  was  delayed  for  a  few  days. 
Unfortunately,  the  body  was  claimed  by  relatives  while  repairs 
to  the  accelerator  were  being  made  and  no  further  data  could 
be  obtained  from  this  specimen.  The  same  procedure  was  followed 
in  the  experiments  on  Cadaver  127.  Compression  fracture  of 
T8  occurred  at  12g  and  there  was  no  evidence  of  hyperextension 
fracture . 

(b)  Strain  histories  in  the  various  spinal  modes: 

Some  typical  strain-time  histories  resulting  from  runs 
in  the  3  spinal  modes  are  given  to  demonstrate  the  difference 
in  levels  of  strain  among  these  modes.  Figures  3-1  and  3.2 
are  examples  of  two  runs  at  ilg  on  Cadaver  1584  in  the  hyper- 
extended  and  erect  modes  respectively.  Fracture  of  T12 
occurred  in  the  erect  mode  at  the  instant  the  shoulder  harness 
load  reached  a  maximum  and  all  strain  gages  indicated  peak 
values.  The  sharp  drop-off  in  the  T12  strain  trace,  shown  in 
Figure  3-2,  is  a  reliable  indicator  of  fracture  which  was 
subsequently  confirmed  by  roentgenographic  examination. 
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Or  particular  interest  is  a  record  of  a  hyperextension 
fracture  of  T12  and  L2  in  Cadaver  099  at  20g.  This  is  depicted 
in  Figure  3- 3-  The  instant  of  fracture  was  about  120  msec, 
after  the  end  of  the  acceleration  pulse.  The  cadaver  was 
undergoing  a  6-9  g  decel eration  when  the  fracture  took  place. 
Figure  3-4  is  a  24 . 5g  fun  in  the  hyperextended  mode.  There 
was  no  fracture.  The  2  strain  gages  on  the  anterior  aspect 
of  L4  were  labelled  as  RAL4  and  LAL4  for  the  right  and  left 
gages.  In  this  cadaver  the  anterior  ligament  w as  intact. 

A  strain-time  history  in  the  flexed  mode  is  shown  in 
Figure  3*5*  It  is  believed  to  be  the  first  of  its  kind  to  be 
reported.  There  are  the  customary  first  and  second  peaks  as 
reported  by  Vulcan  and  King  [19]  although  the  first  peak  is 
not  very  well  defined  in  Figure  3-5.  This  is  due  to  the 
relatively  low  rate  of  onset  used  for  this  series  of  experi¬ 
ments  and  to  the  extensive  effect  of  bending  on  the  strain 
g age  output.  The  phenomenon  to  be  noted  in  the  data  is  a 
third  strain  peak  occurring  at  the  end  of  the  acceleration 
pcise.  Its  magnitude  is  generally  greater  or  equal  to  the 
second  peak.  The  strap  load  cell  indicated  that  the  torso 
did  not  rotate  forward  far  enough  to  cause  any  significant 
force  in  the  shoulder  harness.  However,  when  fracture 
occurred,  there  was  a  sharp  rise  in  the  shoulder  strap  load 
as  shown  in  Figure  3-6  for  a  different  cadaver. 

The  development  of  tension  in  the  vertebral  bodies  during 
the  initial  60  msec,  of  the  acceleration  pulse  occurred  in 
several  runs  with  the  spine  in  the  hyperextended  mode. 
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Generally,  the  magnitude  was  less  than  300p.  It  was  impossible 
to  eliminate  tensile  strains  from  the  vertebrae  of  Cadaver  062. 
The  block  location  was  changed  several  times  to  no  avail.  At 
l8g,  with  the  block  opposite  T12,  a  hyperextension  fracture 
of  LI  occurred.  The  analog  data  are  shown  in  Figure  3-7.  In 
the  erect  mode,  the  vertebrae  were  usually  in  compression 
throughout  the  acceleration  pulse.  However,  tension  did 
develop  during  several  runs  on  Cadavers  002  and  06l.  The 
tensile  strains  were  all  below  700y. 

A  summary  of  the  data  for  the  12  cadavers  used  are  given 
in  Tables  3  through  1*1.  For  the  last  two  cadavers,  (Nos.  125 
and  127),  was  the  only  vertebra  gaged  for  reasons  stated 
previous ly . 

(c)  Examination  of  spinal  segments: 

Some  typical  roentgenograms  of  spinal  segments  removed 
from  cadavers  after  fracture  had  occurred  are  shown  in  Figures 
3.8  through  3-12.  Cadavers  930  and  017  were  run  in  all  three 
modes.  Fractures  resulted  in  the  flexed  mode.  Figure  3-8  shows 
an  anterior  wedge  fracture  of  Til  of  Cadaver  930  and  Figure  3.9 
shows  compression  fractures  of  L2  and  rJl.  An  anterior  com¬ 
pression  fracture  of  Til  of  Cadaver  001  is  shown  in  Figure  3-10 
and  wedge  fractures  of  T9  and  T10  of  Cadaver  1615  are  shown 
in  Figure  3-11*  These  fractures  resulted  from  runs  in  the  erect 
mode.  Figure  3-12  is  an  example  of  a  hyperextension  fracture 
resulting  from  an  l8g  run  in  the  hyperextended  mode.  Cadaver 
062  sustained  a  separation  of  the  superior  end  plate  of  LI. 
Similar  fractures  occurred  in  Cadaver  095  which  was  also  run 
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solely  in  the  hyperextended  mode.  There  is  no  radiological 
evidence  of  damage  to  the  posterior  structures  of  the  vertebrae 
that  were  run  in  the  hyperextended  mode. 

4 .  Analysis  of  Data 

(a)  Acceleration  levels  and  spinal  inodes  at  fracture: 

Of  the  four  cadaver:'  run  exclusively  in  the  hyperextended 
mode,  there  were  2  hypero/tens i on  fractures  at  18  and  2Cg,  one 
compression  fracture  at  12g  and  one  which  did  not  fracture  at 
24.0 g.  Since  the  hyperextcusior;  fractures  are  possibly  due 
to  the  disruption  of  the  ..ulterior  ligament  by  strain  gages, 
it  would  be  a  safe  estimate  to  take  the  average  of  these  four 


vice J erat ion  values  as  the  fracture  lev7  !.  in  the  hyperextended 
mode  for  comparison  with  '.hot  in  the  ot.v'r  modes.  This 
average  is  l8.bg.  The  iv<w  .'-Level  .  n  the  erect  mode  is 
•  o' .  -g  for  6  car'TM  -v  .  in.  t  *'•  - '  'll  xed  mod-',  the  average  of 
1  hr-  o  fracture  levels  is  <) .  Og.  Table  1  ‘j  is  a  summary  of  the 
Of  -rat  tons  at  fracture  ir.  the  j  nodes. 

■-  *  udent  *  s  -i  o t  war:  perf  orm-d  for  the  fracture  levels 


‘.•“tween  •  tie  various  spina  i  modes .  The-  result  is  given  in 
Table  16.  The  difference  in  g-level  betv/een  the  hyperextended 


o.jd:  ar.d  the  other  tv/o  modes  we  re  four  J  to  bo  statistically 
significant  (P  <  c,.0%).  It  should  be  noted  that  the  null 


hypothesis  was  rejected  despite  conservative  estimates  made 
for  the  fracture  level  in  the  hyperextended  mode.  The 
average  age  of  the  cadavers  in  the  hyperextended  mode  group 
was  61-5  years,  while  in  the  erect  and  flexed  mode  groups, 


TABLE  lr3 


Fracture  Level  a  in  the  3  Spinal  Modes 


Mode 

Fracture  Lev-' 
(?) 

;  no.  or 

Cadavers 

Average  Age 
(years) 

Hyperextended 

1 8 .  r, 

/< 

61.5 

Erect 

11  .6 

A 

61.0 

Flexed 

9.0 

3 

5^.3 

51 

it  was  ol  and  5^-3  years  respectively.  There  is  thus  no 
significant  difference  in  age  among  the  groups  to  confuse 
the  data.  Furthermore,  when  a  cadaver  was  run  in  all  or  two 
of  the  spinal  modes,  fracture;  always  occurred  in  the  mode  in 
which  the  cause  could  be  attributed  to  forward  bending  of 
the  vertebral  column.  The  order  of  testing  was  randomized 
so  as  not  to  bias  the  data  toward  any  particular  mode.  The 
result  confirms  the  proposition  that  the  cause  of  fracture 
is  a  combination  of  axial  compression  and  bending,  resulting 
in  the  common  occurrence  of  anterior  fractures  in  pilots  who 
ej  ect . 

(b)  Reduction  in  strain: 

The  effectiveness  of  the  hyperextensicn  device  can  be 
further  demonstrated  by  the  calculation  of  the  percentage 
reduction  in  compressive  strain  as  a  result  of  its  use. 

There  were  32  runs  on  6  cadavers  in  which  strain  data  for 
the  erect  and  hyperextended  modes  at  the  same  g-level  were 
available.  Table  I 7  is  a  listing  of  the  percentage  reduction 
in  strain  for  the  various  vertebra  that  were  gaged  anteriorly. 
The  overall  percentage  reduction  was  ,b%  with  a  range  of 
2b. 1%  for  AT 11  to  65.6%  for  AL3- 

In  order  to  assess  the  statistical  significance  of  this 
reduction,  a  t-test  was  carried  out  for  each  vertebra.  The 
difference  in  strain  between  the  erect  and  hyperextended  mode 
was  computed,  from  which  the  value  of  t  and  the  probability, 

P,  that  the  null  hypothesis  holds,  were  obtained  for  each 
g-level.  Table  18  lists  this  information  for  ail  six 
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vertebrae.  The  number  of  observations  at  each  acceleration 
was  limited  due  to  fracture  at  relatively  low  g-level?  and  to 
the  different  gaging  patterns  used  on  two  of  the  cadavers. 

The  t-test  was  carried  out  only  when  there  are  3  or  more 
pairs  of  data.  In  general,  the  differences  were  significant 
for  the  lumbar  vertebra  (P  <  5 %) ■  For  T12,  the  reduction  in 
strain  is  still  acceptable  since  P  £  b% .  However,  the  observed 
differences  for  Til  indicate  that  the  hyperextension  device 
may  be  beneficial  in  only  4  out  of  5  cases. 

There  was  also  a  consistent  reduction  in  strain  in  the 
lateral  gages  averaging  2H.7%  for  iy  sets  of  data  from  6  pairs 
of  runs  in  the  erect  and  hyperextended  mode.  A  reduction  of 
7.9%  occurred  for  a  pair  of  posterior  gages  placed  on  L2  ir. 
Cadaver  1615. 

A  limited  comparison  could  also  be  made  for  runs  in  the 
flexed  and  hyperextended  mode  and  those  in  the  flexed  and 
erect  mode.  Tables  19  and  20  list  the  percentage  reduction 
in  strain  in  the  anterior  gages  for  these  combinations.  The 
overall  reduction  between  the  flexed  and  hyperextended  mode 
was  53. 7 %  while  that  between  the  flexed  and  erect  mode  was 

33- 1%. 

9 .  Discussion  and  Conclusions 

The  principal  purpose  of  the  experimental  study  was  the 
verification  of  the  mechanism  of  vertebral  fracture  due  to 
caudocephalad  acceleration.  The  existence  of  significant 
bending  stresses  was  first  noted  by  King  and  Vulcan  [  8 ]  and 
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was  reported  in  detail  by  Vulcan  [21].  By  altering  the  spinal 
configuration,  with  a  hyperextension  device  it  was  possible 
to  reduce  vertebral  body  strain  and  raise  the  fracture  g-level. 
The  results  reported  were  found  to  be  statistically  significant 
despite  the  large  biological  variation  in  the  strength  of 
jadaver  vertebrae  and  in  the  degree  of  curvature  of  the  column. 

When  the  fracture  g-levels  in  the  various  modes  were 
analyzed,  it  was  not  possible  to  use  paired  sets  of  data  and 
the  appropriate  values  of  t  were  obtained  from  the  equation 
for  unpaired  data  with  unequal  samples.  In  particular,  for 
hyperextended  mode,  the  ^  fracture  levels  were  obtained  under 
slightly  different  conditions.  The  fact  that  in  two  of  the 
cases  the  anterior  1igament  was  left  intact  while  it  was 
disrupted  in  the  other  two  is  objectionable  from  the  statistical 


viewpoint.  However,  physical  arguments  can  be  used  to  overcome 
such  objections.  Hyperextension  fractures  associated  with  a 
disrupted  anterior  ligament  imply  that  the  g-level  for  com¬ 
pressive  failure  can  only  be  higher  than  that  used  in  the 
computations,  especially  when  one  of  these  fractures  took 
place  during  deceleration  of  the  sled.  Similarly,  when  the 
abdominal  cavity  was  not  eviscerated,  the  fracture  level  could 
be  lowered  somewhat.  Therefore,  the  actual  difference  and  the 
probability  that  it  did  not  occur  by  chance  can  only  be  higher 
than  that  given  in  this  report. 

The  significant  reduction  in  strain  along  the  anterior 
aspect  of  the  vertebral  bodies  is  further  evidence  of  the  role 
of  bending  as  a  cause  for  vertebral  fracture  and  of  the 


effectiveness  of  the  hyperextension  device.  It  is  unfortunate 
that  the  t-test  can  only  be  carried  out  to  9g  and  that  the 
number  of  observations  at  each  acceleration  was  small.  How¬ 
ever,  in  spite  of  this  small  sample  size,  it  was  still  possible 
to  reject  the  null  hypothesis  for  five  of  the  six  vertebrae 
under  study. 

This  series  of  experiments  indicates  that  the  reduction 
in  strains  and  hence  the  increase  in  acceleration  level  at 
fracture  as  a  result  of  the  use  of  a  hyperextension  device, 
cannot  be  attributed  to  any  one  factor.  The  initial  pre¬ 
tension  induced  along  the  anterior  aspects  of  the  vertebral 
bodies  due  to  hyperextension  of  the  spine  causes  the  vertebral 
bodies  to  act  as  prestressed  materials,  hence  raising  the 
compression  level  for  fracture.  However,  this  would  not 
decrease  the  anterior  strains  as  found  in  the  experiments. 

If  bending  only  was  the  cause  of  fracture,  by  changing  the 
eccentricity  of  torso,  we  would  expect  a  redistribution  of 
the  compressive  load  on  the  vertebral  bodies,  hence  reducing 
the  anterior  strain  which  is  a  combination  of  pure  compressive 
strain  and  pure  bending  strain.  However,  this  redistribution 
should  not  affect  the  lateral  strains  on  the  vertebral  bodies, 
as  they  would  be  a  function  of  pure  axial  load  only.  The 
fact,  that  a  significant  reduction  of  the  lateral  strains 
was  found  due  to  hyperextension,  suggests  a  decrease  in  the 
net  axial  compressive  load  on  the  vertebral  body.  This 
decrease  in  the  axial  load  has  to  be  transmitted  to  another 
structure.  One  possibility  is  the  transmission  of  the  load 
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to  the  seat  frame  via  the  hyperextension  device  used.  An 
evaluation  of  this  was  made  by  making  successive  runs,  one 
in  the  erect  and  one  in  the  extended  mode,  on  the  same 
cadaver  at  the  same  acceleration  level,  using  the  seat  pan 
load  cell  measurement  in  the  erect  mode  as  a  control. 

Results  of  these  runs  showed  a  peak  load  cell  value  of 
599  lbs.  for  the  erect  mode  and  592  lbs.  for  the  hyper- 
extended  mode.  This  indicates  that  the  hyperextension 
device  does  not  support  appreciable  vertical  loads  and 
that  the  decreased  portion  of  the  vertebral  body  axial 
load  has  to  be  transmitted  through  a  structure  in  the  spine 
itself.  This  structure  is  the  lamina  via  the  articular 
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Fig.  3-3-  Oscillograph  Record  of  Run  144,  Cadaver  095,  Kyperextended  Mode 
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CHAPTER  IV 

THE  ROLE  OF  ARTICULAR  FACETS  DURING  +  ACCELERATION 

1 .  Introduction 

In  a  typical  vertebra,  the  articular  facets  or  processes 
are  located  near  the  junction  of  the  pedicles  with  the  lamina. 

The  pair  of  superior  facets  spring  upward  from  the  pedicles 
and  face  in  the  general  posterior  direction  while  the  two 
inferior  facets  project  downward  from  the  lamina  and  face 
anteriorly.  A  view  of  a  typical  vertebra  is  shown  in  Figure  4.1. 
The  articular  surfaces  are  lined  with  hyaline  cartilage  and 
form  a  plane  synovial  joint. 

It  is  obvious  from  this  anatomical  arrangement  that  the 
overlapping  facets  perform  the  important  function  of  limiting 
rotation  and  of  preventing  one  vertebra  from  sliding  with 
respect  to  its  adjacent  vertebrae.  The  question  of  whether 
they  are  capable  of  transmitting  compressive  loads  in  the 
longitudinal  direction  of  the  vertebra]  column  has  never' 
really  been  answered.  In  most  texts  of  anatomy,  the  vertebral 
body  is  considered  to  be  the  weight-bearing  structure  of  the 
column  [  3,7  ].  These  references  exemplify  opinions  expressed 
in  1948  and  1970  respectively.  The  facets  have  been  said  to 
carry  no  load  at  all  [  2,6  ].  On  the  other  hand,  Strasser  [  1 8 ] 
and  Nachemson  [13]  have  indicated  that  the  facets  can  support 
a  portion  of  the  load  borne  by  the  spine.  However,  in  a  later 
paper,  Nachemson  [  1 4 ]  retracted  his  earlier  statement  and 
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declared  that  the  facets  carry  no  load.  His  studies  were 
based  on  the  measurement  of  intradiscal  pressures  in  isolated 
spinal  segments  which  were  subjected  to  axial  loads  while  the 
disc  was  tilted  up  to  5° • 

Ewing  et  al.  [5]  proposed  a  hypothesis  for  the  mechanism 
of  vertebral  fracture  during  +  G  impact  acceleration.  The 

4j 

facets  act  as  motion  limiters  preventing  the  posterior  struc¬ 
tures  of  the  vertebra  from  being  displaced  as  much  as  the 
vertebral  body  is  in  front.  In  addition,  the  vertebral  column 
is  subjected  to  eccentric  loading  and  sustains  a  significant 
amount  of  bending  [20].  Hence,  the  head  and  torso  tend  to 
flex  forward,  resulting  in  the  anterior  wedge  fractures 
commonly  noted  among  pilots  who  eject  from  disabled  aircraft. 

It  was  shown  in  Chapter  3  that  the  level  of  fracture  can  be 
substantially  increased  from  10  to  l8g  in  embalmed  human  cada¬ 
vers  by  moderately  hyperextending  the  vertebral  column.  A  2  in. 
thick  block,  6  in.  high  was  placed  on  the  seat  back  with  the 
•enberline  of  the  block  opposite  LI  to  effect  the  hyperextension. 

There  was  a  concomitant  decrease  in  strain  along  the 
anterior  aspects  of  the  lower  thoracic  ar d  the  lumbar  vertebrae 
but  the  seat-pan  load  was  not  altered  as  a  result  of  this 
hyperextension.  These  results  point  to  the  load-bearing 
capability  of  the  motion  limiters,  namely,  the  facets.  In 
this  chapter,  experimental  evidence  of  the  existence  of  a 
dual  load  path  along  the  spine  Is  documented.  These  findings 
provide  a  better  understanding  of  the  mechanism  of  injury 
which  is  also  delineated. 


2. 


Experimental  Methods  and  Equipment 
Several  techniques  were  used  to  deduce  the  load-bearing 
capability  of  the  articular  facets,  since  direct  measurement 
of  force  in  a  limited  space  environment  is  still  beyond  the 
state  of  the  art. 

Strain  gages  were  an  effective  means  of  providing  a 
qualitative  indication  of  facet  load.  They  were  mounted  on 
the  pedicles  and  lamina.  Methods  of  installing  them  on 
vertebral  surfaces  were  de /el oped  to  allow  the  measurement  of 
strain  just  about  anywhere  on  the  external  surface  of  a  verte¬ 
bra.  A  detailed  description  of  the  techniques  used  to  install 
strain  gages  in  vertebral  bodies  has  been  given  by  King  and 
Vulcan  [9].  The  same  methods  are  employed  for  installing 
them  on  the  pedicles,  for  which  an  anterior  approach  is  used. 

To  mount  them  on  the  posterior  surface  of  a  lamina,  a  posterior 
approach  is  taken  but  the  basic  techniques  remain  unchanged. 

Quantitative  measures  of  facet  load  were  obtained  by 
moans  of  an  intervertebral,  load  ceil  (JVLC)  designed  to  fit 
under  a  disc  or  a  vertebral  body.  It  can  measure  both  axial 
force  and  the  eccentricity  of  that  force  with  respect  to  its 
geometric  center.  figure  A.?  shows  the  first  model  of  an  IVLC. 
It  was  almost  1  in.  thick  and  is  believed  to  be  the  first 
transducer  capable  of  actually  measuring  the  load  carried  by 
a  vertebral  body  in  an  intact  vertebral  column.  A  thinner 
version  is  shown  in  Figure  A. 3.  This  load  cell  was  designed  to 
fit  above  a  lumbar  disc,  replacing;  the  inferior  segment  of  a 
vertebral  body.  A  double-b 1 aded  rotary  saw  was  constructed 
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to  cat  slots  of  a  precise  width  across  the  spine  to  accommodate 
the  IVLC  without  damaging  the  neural  arch.  The  lip  shown  in 
Figures  A.  2  and  .  3  enables  a  strap  to  hold  the  IVLC  in  place 
during  the  experimental  run.  Figure  ^ . k  shows  the  IVLC  in 
position  in  the  lumbar  spine  of  a  cadaver.  A  properly  in¬ 
stalled  IVLC  will  not  result  in  a  change  in  length  of  the 
vertebral  column  and,  if  the  thinner  model  is  used,  the 
mobility  of  the  column  is  also  not  affected.  When  an  IVLC 
is  used,  strain  gages  are  installed  on  both  the  anterior  and 
posterior  (lamina)  aspects  of  the  vertebrae  adjacent  to  it 
to  obtain  correlation  of  facet  load  with  strain. 

Other  instrumentation  consists  of  load  cells  on  the  seat 

pan  and  back,  and  for  the  shoulder  harness.  Two  sets  of 

mutually  perpendicular  pairs  cf  uniaxial  accelerometers  were 

mounted  on  the  head  during  some  of  the  runs  to  study  the 

dynamics  of  head  motion  during  +  G  acceleration.  The  sled 

z 

acceleration  was  also  monitored. 

The  fully  instrumented  cadaver  was  placed  in  the  seat  of 
a  vertical  accelerator  for  testing.  This  facility  is  housed 
in  an  3-story  elevator  shaft  of  the  School  of  Medicine  at 
Wayne  State  University.  The  sled  is  accelerated  over  a  stroke 
of  8  feet  and  then  gradually  brought  to  rest  by  air  brakes 
over  some  30  to  40  feet.  The  acceleration  pulse  is  approximately 
trapezoidal  In  shape,  the  rate  of  onset  and  the  magnitude  of 
the  plateau  being  variable.  Details  of  the  accelerator  have 
been  described  by  Patrick  [16]. 


The  restraint  system  consisted  of  an  automotive  lap  belt 
under  a  regular  U.S.A.F.  lap  belt  and  shoulder  harness  combina¬ 
tion  and  leg  straps.  The  wrists  were  tied  together  and  anchored 
to  the  seat  pan  by  means  of  a  single  rope  going  through  an  eye 
bolt  on  the  seat  pan  and  looping  around  the  lap  belt.  This 
was  done  to  prevent  flailing  of  the  arms  and  to  keep  the  lap 
belt  in  position  on  the  pelvis  during  pretensioning  of  the 
shoulder  harness.  The  head  was  unrestrained. 

Electronic  instrumentation  consisted  of  12  channels  of 
bridge  balance  and  carrier  amplifier  units  (Heiland),  a  14- 
channel  tape  recorder  (Ampex),  and  a  24-channel  light-beam 
recorder  (Visicorder) . 

This  chapter  covers  the  results  of  approximately  82  runs 
made  on  six  different  cadavers.  Only  strain  data  were  acquired 
from  the  first  cadaver,  while  both  IVLC  and  strain  data  were 
obtained  from  the  other  five.  Table  1  lists  the  pertinent 
information  on  the  cadavers  used,  the  test  conditions  and  the 
location  of  strain  gages  and  the  IVLC.  The  ability  to  vary 
the  seat-back  angle  from  0°  to  20°  rearward  is  a  recent 
modification  to  the  vertical  accelerator.  This  feature  was 
used  for  the  first  two  cadavers  and  provided  additional  evi¬ 
dence  of  the  load-bearing  capability  of  the  facets.  The 
hyperextension  block  was  used  t.o  change  the  spinal  curvature 
and  thus  the  role  of  the  facets,  if  they  were  indeed  able  to 
carry  the  load.  In  subsequent  discussions,  the  hyperextended 
and  erect  spinal  modes  refer  to  runs  made  wi„h  and  without 
the  hyperextension  block  respectively.  The  notation  for 
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strain  gages  used  on  the  oscillograph  records  are  as  follows: 

A  denotes  an  anterior  gage,  P  a  posterior  gage,  and  D  any 
gage  on  the  lateral  aspects  of  the  vertebral  body  or  on  the 
neural  arch. 

3 .  Vertebral  Strain  Data 

An  exploratory  study  of  facet  load  was  made  on  Cadaver 
20c/’.  A  strain  gage  was  placed  on  the  lateral  surface  of  the 
neural  arch  of  LI,  near  the  body.  Its  sensitive  axis  was 
inclined  at  approximately  45°  with  respect  to  the  longitudinal 
axis  of  the  body.  In  this  configuration,  it  measures  shear 
strain  on  the  arch  and  is  an  indicator  of  any  shift  in  load 
from  the  body  to  the  facets  or  vice  versa.  The  cadaver  was 
subjected  to  a  6-g  pulse  at  various  seat-back  angles  with  the 
spine  in  the  hyperextended  as  well  as  in  the  erect  modes. 

Figure  4.5  shows  the  strain  data  for  a  run  (No.  250)  made  with 
the  seat  back  vertical  (0°)  and  in  the  erect  mode.  The  shear 
gage  on  LI  (DL1)  indicated  a  slight  compressive  strain  initially 
and  then  a  large  tensile  strain  during  the  pulse.  The  peak 
tension  coincides  in  time  with  the  peak  compression  along  the? 
anterior  aspect  of  LI  (AL1).  The  results  for  a  similar  run 
made  in  the  hyperextended  mode  are  shown  in  Figure  4.6  (Run  249). 
In  this  case  there  is  a  slight  compression  at  AL1  and  a  sig¬ 
nificant  compressive  strain  at  DL1 ,  suggesting  a  shift  in  load 
from  the  vertebral  body  to  the  lamina  and  the  facets. 

The  data  for  a  hyperextended  and  an  erect  mode  run,  made 
at  a  seat-back  angle  of  20°  rearward,  are  shown  in  Figures  4.7 
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(Run  243)  and  4.8  (Run  294)  respectively.  In  both  cases  DL1 
was  in  compression  throughout  the  entire  acceleration  pulse. 
However,  in  the  erect  mode  the  peak  strain  was  20  Op  while  it 
was  800p  in  the  hyperextended  mode.  Further  evidence  of 
transfer  of  load  to  the  posterior  vertebral  structures  is 
given  by  the  absence  of  a  second  peak  in  the  anterior  strains 
which  were  evident  Iji  runs  made  with  the  seat  back  vertical. 

The  rearward  inclination  of  the  seat  back  and  the  use 
of  the  hyperextended  block  have  the  same  effect  as  far  as  the 
facets  are  concerned.  In  each  case,  there  is  a  decrease  in 
anterior  strain  and  a  change  in  sign  of  the  shear  gage  output 
compared  to  that  for  a  run  made  at  0°  seat-back  angle  and  in 
the  erect  mode. 

To  document  this  phenomenon  more  fully,  the  results  of 
runs  made  at  a  seat-back  angle  of  10°  are  discussed.  As  shown 
in  Figure  4.9  (Run  247),  hyperextension  of  the  spine  causes  DL1 
to  stay  in  compression  throughout  the  run  and  the  anterior 
gages  do  not  show  a  second  peak.  However,  in  the  erect  mode, 
DL1  changes  sign  from  compression  to  tension  during  the  run 
with  a  concomitant  appearance  of  a  second  peak  for  the  anterior 
gages,  as  shown  in  Figure  4.10  (Run  248).  Peak  compression  of 
AL1  coincides  again  with  peak  tension  of  DL1.  At  a  seat-back 
angle  of  5° ,  the  same  phenomena  are  repeated. 

To  investigate  further  the  role  of  the  facets,  they  were 
removed  by  dissection  above  and  below  LI.  The  lamina  and 
spinous  process  of  LI  were  also  removed.  A  pair  of  6-g  runs 
were  made  in  the  erect  and  hyperextended  mode  before  the 
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dissection  was  carried  out  and  another  pair  was  made  with 
the  same  input  acceleration  after  the  posterior  structures 
were  removed.  The  seat  back  was  vertical  for  all  runs. 
Figures  4.11  (Run  253)  and  4.12  (Run  257)  give  data  fur  runs 
made  in  the  erect  mode,  with  and  without  facets  respectively. 
The  anterior  strains  at  T12,  1,1  and  L2  show  a  marked  increase 
after  the  removal  of  the  facets.  In  fact,  for  AL1 ,  there  was 
a  121%  increase  in  compressive  strain.  Moreover,  the  strain 
pattern  for  the  shear  gage  (DL1)  was  changed  drastically. 

Its  output  dropped  to  nearly  zero  from  a  peak  tensile  strain 
of  1240p  when  the  spine  was  intact. 


A  comparison  of  the  effect  of  facet  removal  in  the  hyper- 
extended  mode  can  be  made  from  Figures  4.13  (Run  254)  and  4.14 
(Run  258).  Run  254  was  made  with  the  spine  intact.  The 
increase  in  AL1  strain  was  90%,  as  a  result  of  facet  removal. 
There  are,  however,  two  inconsistencies  in  the  data.  The 
decrease  in  shear  strain  was  only  18%  whereas  it  was  again 
expected  to  drop  to  almost  zero.  AL2  showed  an  11%  decrease 
in  strain  instead  of  an  increase.  A  possible  explanation  of 
these  effects  is  the  impingement  of  the  vertebrae  against  the 
hyperextension  block,  the  centerline  of  which  was  located 
opposite  LI. 

These  results  led  to  a  more  direct  approach  to  the  problem. 
If  the  posterior  structures  constitute  a  load  path,  a  strain 
gage  on  the  lamina  should  provide  considerable  clarification 
of  the  situation.  A  single  gage  on  the  lamina  of  L2  of 
Cadaver  2062  constituted  an  initial  feasibility  study  of 
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acquiring  useful  data  from  this  location.  There  were  again 
anterior  gages  from  T8  through  L3  and  the  experimental  plan 
was  to  make  runs  at  various  seat-back  angles  in  the  erect 
and  hyperextended  mode.  It  was  found  that  the  output  of 
this  gage  was  less  than  600p  for  a  6-g  run.  However, 
interesting  trends  could  be  picked  out.  In  the  hyperextended 
mode,  the  strain  was  generally  compressive  or  remained  com¬ 
pressive  longer  than  in  the  mode.  In  this  mode,  strains 

that  were  initially  compressive  frequently  became  tensile 
during  the  latter  half  of  the  pulse. 

A  more  extensive  stuiy  of  laminar  strain  was  carried  out 
on  this  cadaver  by  installing  additional  gages  on  the  lamina 
of  T10  and  Til.  The  series  of  runs  was  repeated  using  these 
3  posterior  gages  and  anterior  gages  from  T9  through  L3.  The 
input  acceleration  level  was  raised  to  8-g.  Data  for  the 
erect  and  hyperextended  mode  are  shown  in  Figures  4.15  (Run  283) 
and  4.16  (Run  282)  respectively.  The  seat  back  was  vertical. 

The  difference  in  response  of  the  posterior  gages  is  quite 
obvious.  In  the  erect  mode  (Figure  4.15)  PT10  and  PT11  re¬ 
mained  in  tension  throughout  the  pulse  while  PL2  was  initially 
in  compression  and  went  into  tension  at  about  50  msec  after 
the  onset  of  acceleration.  As  a  result  of  hyperextension,  PL2 
remained  in  compression  for  the  entire  duration  of  the  pulse, 
while  PT10  and  PT11  underwent  a  change  in  sign  from  compression 
to  tension,  as  shown  in  Figure  4.16.  A  marked  reduction  in 
anterior  strain  is  also  quite  evident. 
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Rearward  rotation  of  the  seat  back  did  not  alter  the 
posterior  strain  pattern  which  was  predominantly  dependent 
upon  the  curvature  of  the  spine.  The  response  of  the  F-gages 
at  20°  was  the  same  as  that  at  0°  in  terns  of  the  strain¬ 
time  history  and  the  reversal  in  sign.  The  difference  in 
response  between  the  troracic  and  lumbar  gages  is  due  to  the 
location  of  the  vertebrae  on  different  curves  of  the  column. 
However,  the  qualitative  evidence  from  these  data  point  to  a 
transfer  of  load  from  the  facets  to  the  vertebral  body  as  the 
head  and  torso  rotate  forward. 

a .  Intervertebral  Load  Cell  Data 

The  purpose  of  the  IVLC  was  to  measure  the  load  carried 
by  the  vertebral  body  and  to  compare  it  with  that  borne  by 
the  column,  the  total  spine  load.  The  latter  was  taken  to 
be  proportional  to  the  measured  seat  pan  load.  The  ratio 
used  was  the  weight  of  the  torso  above  the  IVLC  to  the  total 
body  weight.  Any  difference  between  the  spine  load  and  the 
IVLC  output  is  the  load  carried  by  the  articular  facets,  the 
facet  load.  Justifications  Tor  the  validity  of  this  method 
of  deducing  facet  load  are  provided  in  the?  discussion  section 
of  the  paper. 

IVLC  data  were  obtained  from  five  different  cadavers. 

The  initial  thicker  version  was  used  in  the  first  two  cadaver 
while  the  second  thinner  model  was  employed  during  runs  made 
on  the  third.  A  typical  oscillograph  record  of  IVLC  output 
is  shown  in  Figure  4.17  (Run  304).  This  was  a  10-g  run  made 
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on  Cadaver  2062  in  the  erect  mode.  The  IVLC  replaced  the 
inferior  portion  of  the  L3-L-4  disc  and  the  superior  segment 
of  L4.  The  neural  arch  of  L4  was  left  intact.  Figure  4.18 
shows  the  computed  spine  load  and  facet  force  which  is  the 
difference  between  the  spine  load  and  the  intervertebral 
body  force  (Run  30*0-  The  facets  were  in  compression  for 
the  first  125  msec  of  the  pulse.  However,  as  the  head  and 
torso  rotated  forward,  the  facets  unloaded  and  went  into 
tension,  resulting  in  an  intervertebral  body  force  larger 
than  the  total  spine  load.  The  facet  load  and  PL3  strain 
also  show  good  correlation,  with  the  zero  crossover  point 
of  both  traces  occurring  almost  simultaneously.  To  minimize 
the  number  of  figures,  description  of  subsequent  IVLC  runs 
will  be  accompanied  by  plots  of  spine,  facet,  and  inter¬ 
vertebral  body  loads  and  oscillograph  records  will  be 
omitted.  Figure  4.19  (Run  303)  shows  force  and  strain  data 
for  a  10-g  run  made  in  the  hyperextended  mode.  The  compres¬ 
sive  facet  load  was  larger  in  magnitude  and  longer  in  dura¬ 
tion  than  that  for  the  erect  mode  (Figure  4.18,  Run  304). 

The  IVLC  (1st  model)  was  used  again  in  Cadaver  2093. 

It  was  located  between  L3  and  L4.  The  inferior  segment  of 
L3  and  the  superior  portion  of  the  L3-L4  disc  was  removed  to 
accommodate  the  IVLC.  A  2.5-g  run  was  made  in  the  erect  mode 
to  study  spinal  response  near  the  1-g  environment.  The  data 
are  shown  in  Figure  4.20  (Run  323)-  The  facets  were  taking 
about  50$  of  the  total  spine,  but  they  did  not  unload  and 
go  into  tension.  The  strain  pattern  of  PL3  again  followed 


that  of  the  facet  load  very  closely.  At  9  g,  the  proportion 
of  spine  load  carried  by  the  facets  was  about  3 5 % >  as  shown 
ir.  Figure  4.21  (Run  324).  The  spine  was  in  the  erect  mode. 

The  facets  and  PL3  were  in  compression  throughout  the  pulse, 
indicating  that  unloading  never  took  place.  The  location 
and  thickness  of  the  IVLC  may  have  decreased  the  mobility  of 
the  vertebrae  and  the  unloading  phenomenon  may  also  be  de¬ 
pendent  on  the  curvature  of  the  spine. 

A  series  of  runs  was  carried  out  on  three  cadavers  using 
the  thinner  model  of  the  IVLC.  Its  overall  thickness  was 
0.4  in.  and  it  was  placed  above  the  L3-L4  disc  by  replacing 
the  inferior  segment  of  L3.  The  disc  was  virtually  intact. 
These  cadaver  runs  were  used  to  validate  mathematical  models 
of  the  spine,  described  in  Chapter  V,  hence  the  results  of 
all  the  runs  made  are  not  shown  in  this  chapter.  Each  cadaver 
was  run  at  6.8  and  10  g's  acceleration  in  the  erect  and  the 
hyperextended  modes. 

Figure  4.22  (Run  378)  shows  an  8-g  run  made  in  the  erect 
mode  on  Cadaver  2231.  The  facets  unloaded  rather  early  in 
the  acceleration  pulse  (at  about  60  msec),  but  there  was 
confirmation  from  PL3  strain.  For  an  identical  run  made  in 
the  hyperextended  mode,  the  facets  remained  in  compression 
for  the  entire  pulse,  as  shown  in  Figure  4.23  (Run  377).  PL3 
and  PL4  strains  were  also  in  compression  throughout  the  pulse. 
The  total  spine  load  was  the  same  in  these  two  runs,  but  the 
intervertebral  body  force  decreased  by  about  400  lb.  as  a 
result  of  hyperextension.  The  results  of  the  runs  made  at 
6  and  10  g’s  on  this  cadaver  were  essentially  the  same. 
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The  results  obtained  from  Cadavers  2209  and  2^13  were 
essentially  the  same  as  those  obtained  from  Cadaver  2231. 

The  results  of  all  the  eighteen  runs  mode  on  the  above  three 
cadavers  are  shown  in  the  chapter  on  the  experimental  verifi¬ 
cation  of  the  mathematical  model. 

5 .  Discussions  and  Conclusions 

Experimental  results  have  been  presented  to  document 
the  role  of  the  articular  facets  during  +  0  impact.  Qualita¬ 
tive  evidence  in  the  form  of  strain  data  from  the  anterior 
and  posterior  aspects  of  both  thoracic  and  lumbar  vertebrae 
indicate  the  existence  of  a  dual  load  path  along  the  vertebral 
column.  This  led  to  the  development  of  an  IVLC  which  was  used 
to  obtain  quantitative  data  supporting  the  claim  that  both 
tensile  and  compressive  loads  can  be  transmitted  via  the 
facets  or  the  posterior  structures  of  the  lumbar  vertebrae. 

The  complexity  of  the  posterior  structures  and  the 
limitation  of  space  in  and  around  the  joints  of  the  facets 
precluded  a  direct  measurement  of  facet  load.  The  load- 
bearing  capability  was  deduced  by  comparing  the  total  spine 
load  with  that  obtained  from  the  IVLC.  The  spine  load  was 
taken  to  be  proportional  to  the  seat  pan  load.  That  is,  the 
dynamic  response  of  the  torso  above  the  IVLC  was  assumed  to 
be  the  same  as  that  below  it.  Since  the  IVLC  was  placed 
between  L3  and  L4,  it  can  be  argued  that  there  should  be  very 
little  dynamic  overshoot  from  the  lower  torso  and  legs  and 
that  the  spine  load  should  be  the  difference  between  the  seat 


pan  load  and  the  product  of  sled  acceleration  and  mass  of  the 
lower  torso  and  legs.  However,  the  dynamic  overshoot  of  the 
seat  pan  load  is  10?  or  less  so  that  the  spine  load  computed 
by  either  method  does  not  differ  significantly.  Furthermore, 
during  a  major  portion  of  the  pulse,  the  seat  pan  load  cell 
response  is  almost  flat.  Consequently,  the  facet  load  has 
been  estimated  fairly  accurately  by  means  of  the  IVLC  and  the 
seat  pan  load  cell. 

The  use  of  the  hyperextension  block  and  the  variation 
of  the  seat-back  angle  provided  additional  data  confirming 
the  role  of  the  facets  and  aided  in  the  understanding  of  the 
mechanism  of  injury.  On  the  basis  of  the  data  obtained,  it 
is  now  possible  to  assemble  the  research  results  on  spinal 
injury  and  propose  an  injury  mechanism  which  fully  explains 
the  commonly  observed  anterior  wedge  fractures  sustained  by 
pilots  during  emergency  egress.  Previous  work  by  Vulcan  and 
King  [20]  established  the  fact  that  during  caudo-cephalad 
(+  G  )  acceleration  the  spine  is  subjected  to  both  axial 
compression  and  bending.  The  bending  effects  are  due  to  the 
eccentricity  of  the  torso  with  respect  to  the  spine  and  are 
enhanced  by  the  forward  rotation  of  the  head  and  torso.  Sub¬ 
sequent  work  described  in  Chapter  III  showed  that  fracture  levels 
could  be  raised  significantly  by  placing  a  2-in.  thick  hyper¬ 
extension  block  opposite  LI  on  the  seat  back.  The  curvature 
of  spine  was  altered  by  the  block  but  there  was  insufficient 
movement  to  decrease  the  eccentricity  of  the  torso  and  hence 
the  bending  moment  on  the  spine.  Nevertheless,  the  average 


level  of  fracture  was  raised  from  10  to  18  g  and  the  decrease 
in  anterior  strains  was  generally  significant  at  the  95$ 
confidence  level.  The  hypothesis  that  the  facets  act  as 
motion  limiters  led  to  this  documentation  of  their  role  during 
+  G  acceleration.  In  the  erect  spinal  mode,  the  facets  tend 
to  unload  and  go  into  tension,  causing  the  vertebral  bodies 
to  sustain  more  compressive  load  than  the  total  spine  load. 

This  occurs  when  the  head  and  torso  undergo  maximum  forward 
flexion.  The  anterior  wedge  fractures  are  therefore  the  result 
of  eccentric  compression  coupled  with  the  unloading  of  the 
facets.  In  the  hyperextended  mode,  the  facets  relieve  the 
vertebral  bodies  of  some  of  the  compressive  load  and  thus  it 
was  possible  to  raise  the  fracture  level  by  such  a  considerable 
margin . 

In  summary,  the  following  conclusions  can  be  made: 

1.  The  articular  facets  are  capable  of  bearing  compres¬ 
sive  and  tensile  loads. 

8 .  Strain  gages  wore  employed  to  provide  qualitative 
evidence  of  facet  load. 

3-  Intervertebral  body  force  was  measured  in  an  intact 
spine  during  Impact  acceleration  by  means  of  a 
specially  designed  intervertebral  load  cell  (IVLC). 

4 .  From  the  IVLC  and  seat  pan  load  cell  output,  a  facet 
load  history  was  computed. 

5.  A  better  understanding  of  the  in.iury  mechanism  of 
the  spine  has  been  achieved. 
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6.  Hyperextension  of  the  spine  transfers  more  load  to 
the  facets. 

7-  The  proportion  of  the  load  carried  by  the  facets 
appear  to  increase  with  the  decreasing  g-levels, 
suggesting  that  they  may  also  carry  a  portion  of 
the  static  body  weight  when  the  body  assumes  a 
normal  erect  posture. 


Fig.  4.1.  Top  and  Side  View  of  a  Typical  Vertebra 
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Fig.  U.ii.  Oscillogi'aph  Record 


Fig.  4.12.  Oscillograph  Record  of  Run  257 »  Cadaver  2067  -  Facets  Rerr.oved 
Erect  Mode 
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Fig.  4.16.  Oscillograph  Record  of  Run  282,  Cadaver  2062,  Seat-back  Angle  0 
Hyperext ended  Mode 
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CHAPTER  V 


A  MATHEMATICAL  MODEL  FOR  THE  SPINE 

1 .  Introduction 

In  an  effort  to  explain  vertebral  body  fractures  due  to 
+  G  acceleration,  several  mathematical  models  have  been 

L» 

formulated.  In  recent  years  two  models  by  Vulcan  [21]  and 
Orne  [15]  have  rendered  a  better  understanding  of  the  mechanisms 
of  failure  of  the  human  spine.  Vulcan  considered  the  effects 
of  forward  flexion  of  the  upper  torso,  the  head  and  neck,  the 
curvature  of  the  spine,  and  the  effects  of  the  restraint  system 
on  the  spine  during  +  G  acceleration.  He  established  the 
presence  of  high  bending  moments  in  the  spine  due  to  the  forward 
eccentricity  of  the  upper  torso  center  of  gravity  with  respect 
to  the  centerline  of  the  spine.  The  bending  effect  is  greatest 
when  the  chin  contacts  the  chest  and  the  upper  torso  reaches 
the  maximum  flexion  allowed  by  the  restraint  systems.  An 
important  fact  brought  out  by  his  experiments  and  model  is 
the  role  played  by  the  whipping  of  the  head  during  the  accelera¬ 
tion  pulse.  In  his  mathematical  model  he  considers  the  head 
and  half  the  neck  as  one  rigid  body  connected  by  springs  and 
dampers  to  another  rigid  body  representing  the  torso  from  T1 
to  L*J.  The  acceleration  pulse  is  applied  to  L5.  This  re¬ 
sulted  in  a  four  degree-of- freedom  model. 

Orne  proposed  a  discrete-parameter  model  consisting  of 
alternately  rigid  and  deformable  bodies  simulating  the  behavior 
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of  the  vertebrae  and  discs,  respectively,  the  bodies  being 
arranged  In  such  a  fashion  as  to  describe  the  natural  curva¬ 
tures  of  the  spine.  The  mass  of  the  system  is  assumed  to  be 
concentrated  in  the  rigid  bodies,  which  are  capable  of  three 
degrees-of-freedom  in  the  sagittal  plane  of  the  body.  This 
resulted  in  a  5^  degree-of- freedom  model.  This  model  renders 
the  study  of  stresses  at  different  levels  of  the  spine  in 
greater  detail  than  the  Vulcan  model.  However,  the  head  and 
the  neck  are  represented  by  only  one  rigid  link  ignoring  the 
whipping  effect  mentioned  above.  Another  serious  shortcoming 
of  this  model  is  the  exclusion  of  any  restraint  forces,  such 
as  the  shoulder  strap  and  seat  back  reaction  forces.  The 
resulting  response  is  totally  different  from  that  observed 
experimentally  on  cadavers  or  during  ejections  from  aircraft. 

In  his  model,  although  there  is  flexion  of  the  vertebral 
column,  it  is  accompanied  by  a  rearward  movement  of  the  upper 
torso,  resulting  in  rearward  bending  moments  at  the  lumbar 
vertebrae  rather  than  forward  bending  moments  observed  experi¬ 
mentally  . 

In  view  of  the  fact  that  there  exist  two  load  paths  in 
the  spine,  one  through  the  vertebral  body  and  the  other  through 
the  lamina  via  the  articular  facets,  it  was  decided  that  the 
above  two  models  were  inadequate  to  simulate  the  response  of 
the  vertebral  column  due  to  +  Gz  acceleration.  Hence,  a  two 
dimensional  discrete-parameter  mathematical  model  of  the  spine 
was  formulated  with  the  following  requirements: 
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(i)  The  model  should  predict  the  amount  of  forward 
flexion  of  the  spine. 

(ii)  It  should  account  for  the  natural  curvatures  of 
the  spine. 

(iii)  It  should  account  for  the  eccentric  inertial 
loading  on  the  spine. 

(iv)  It  should  predict  head  and  neck  motions  and  their 
effects  on  the  forces  and  moments  in  the  spine, 

(v)  The  load  transmission  at  each  vertebral  level  has 
to  take  place  through  the  vertebral  body  and  the 
articular  facets,  thus  incorporating  a  parallel 
load  path  with  the  vertebral  body  -  a  factor  not 
considered  in  any  of  the  existing  mathematical 
models  . 

(vi)  The  effects  of  the  restraint  and  support  systems 
must  be  incorporated  to  properly  simulate  the 
ejection  problem. 

(vii)  The  model  should  be  able  to  simulate  off-axis 

impacts  in  the  mid-sagittal  plane  and  predict  the 
response  due  to  any  input  pulse. 

2 .  Development  of  the  Model 

Considering  the  requirements  of  the  model  stated  in  the 


preceding  section,  the  following  assumptions  were  made  in  the 
mathematical  development: 

(i)  The  2k  vertebral  bodies,  the  head  and  the  pelvis 
are  rigid  bodies  constrained  to  move  in  the  mid- 
sagittal  plane. 


(ii)  Each  rigid  body  has  three  degrees-of-freedom  in 

the  mid-sagittal  plane,  two  translational  and  one 
rotational . 

(iii)  The  intervertebral  discs  are  massless  and  deforma¬ 
tion  of  the  spine  takes  place  at  the  discs. 

(iv)  The  discs  are  replaced  by  a  system  of  springs  and 
dampers  -  one  spring  and  damper  for  axial  forces, 
one  spring  and  damper  for  shear  forces  and  another 
spring  and  damper  arrangement  for  restoring  torques 
due  to  relative  angular  motion  between  adjacent 
vertebral  bodies. 

(v)  The  facets  and  laminae  are  springs  connected  to 
the  vertebral  body  by  a  massless  rigid  rod. 

(vi)  Each  rigid  body  is  assumed  to  carry  a  portion  of 
the  torso  weight  which  is  eccentric  with  respect 
to  the  centerline  of  the  spine. 

(vii)  The  rigid  bodies  are  arranged  to  simulate  the  spinal 
curvatures  as  closely  as  possible. 


3 .  Kinematic  Preliminaries 

Figure  5-1  shows  the  initial  configuration  of  two  success! v< 
links  (the  vertebrae)  and  a  deformable  link  (the  disc).  It 
has  been  assumed  that  the  axis  of  any  disc  is  coincident  with 
the  axis  of  the  vertebra  immediately  below  it  in  the  initial 
configuration  for  the  erect  mode  as  shown  in  Figure  5.1.  Also 
it  is  assumed  that  the  disc  is  of  uniform  thickness,  which 
is  not  necessarily  true.  The  two  vertebrae  are  shown  as 
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trapezoids  to  simulate  the  change  of  curvature  of  the  spine. 

This  is  a  valid  assumption  because  it  has  been  observed  that 
in  the  lumbar  region  the  vertebral  bodies  are  wedged  posteriorly, 
whereas  in  the  thoracic  level  they  are  wedged  anteriorly.  The 
position  in  the  mid-sagittal  plane  of  the  center  of  the  i— 
rigid  link  (vertebral  body)  is  determined  by  the  three  generalize- 
coordinates  ui ,  w^  and  0  as  shown  in  Figure  5  •  1  • 

At  time  t  (t  >  0)  ,  the  configuration  of  two  successive 
rigid  links  is  shown  in  Figure  5.2.  The  links  have  undergone 
translations  and  rotations  causing  axial,  shear  and  rotational 
deformations  of  the  discs.  The  chord  length  AB^  of  the  disc 
is  determined  by  the  generalized  coordinates  u^ ,  w^  and  0^ 
of  the  centers  of  the  rigid  links.  The  calculations  necessary 
for  computing  the  deformation  of  the  discs  follows  in  the 


order  used  in  the  computer  program. 

Deformation  of  the  disc: 

From  Figure  5.2 

*1  =  (uA  +  d^  sin  0i)  -  (u1_1  -  di_1  sin  0j__1)  (1) 

x~  =  (w.  -  d.  cos  0.)  -  (w.  ,  +  d.  ,  cos  0,  ,  )  (?) 

2  l  i  l  l-l  l-l  i-1 

ABi  =  xi  +  x2  ^3) 

=  tan-'*'  x^/x^  ( ) 


and 
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AB^  sin  (ai 


6i-i> 


AB^  cos  (a^ 


0i-l} 


The  lengths  AC  and  BC  will  be  used  as  the  criterion  for 
the  development  of  forces  on  the  vertebral  bodies.  Hence, 
at  time  t  =0,  AC  will  be  denoted  by  ACo  and  BC  by  BCq. 

The  time  rate  of  change  of  AC  and  BC  have  to  be  found  to 
generate  viscous  forces.  For  this,  we  have  to  first  compute 

the  rate  of  change  of  AB  and  the  angle  a.  The  time  derivatives 

•  •  •  • 

will  be  referred  to  by  AB,  a,  AC,  and  BC,  etc. 

Differentiating  Equation  (3)  with  respect  to  time,  we  get 


2  ? 
X1  +  X2 


( 2xx  x1  +  2x2  x2) 


=  x.  cos  a,  +  xn  sin  a. 
1  i  2  i 


where  and  x2,  obtained  by  differentiating  equations  (1) 
and  (2),  are 


xn  =  (u.  +  d.  0.  cos  0.)  -  (u1_1  -  d1_1  0i_]  cos  0,  J 


l  ii 


x2  =  (wi  +  d^_  0i  sin  0i)  -  (w1_1  -  di_1  0i_1  sin  0i_1) 

cu  can  be  obtained  by  differentiating  (^),  but  a  simpler 
equation  is  obtained  if  we  consider  Figure  5.3. 

In  this  figure  if  v  t  denotes  the  tangential  velocity  of  BA, 


vat  =  AB1  “i 
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vat  "  x2  cos  ai  ~  xi  sin  ai 


ai  ~  ( x2  cos  ai  ~  x]  sin  ai)/AB1 


Differentiating  (5)  and  (6)  we  get 


AC^  =  AB1  sin  (ajL  -  0i_1) 

+  ABi  cos  (ou  -  ei_1)(a1  -  0  ) 


BC^  =  AB^  cos  (ct1  -  0  ) 

-  AB.  sin  (a.  -  -  0^) 


Geometry  of  the  facets: 


From  Figure  5.2 


yl  “  (ui  -  hi  303  9i>  -  K_  1  -  h. 


1-1  -  hl-i  003  0i-B 


y2  (w^  -  hA  sin  0i)  -  (w1_1  -  hi_I  sin  9i_1) 


A'B,i  -  y22 


4>-t  =  tan  y  /y 


2  J  1 


AlA2i  =  A  *  B ’ 1  sin  (^.  -  01_1) 


^2A31  =  A'B’.  cos  (({,.  -  01_1) 
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4 .  Forces  Developed  In  the  Disc 

The  initial  configuration  of  the  disc  is  shown  in  Figure  5.', 
where  AB  is  the  centerline  of  the  disc.  After  deformation  the 
neutral  axis  AB  takes  the  shape  shown  in  Figure  5-5.  The  force..- 
and  restoring  torque  developed  on  the  disc  are  assumed  to  be 
functions  of  the  change  in  lengths  AC,  BC  and  angular  deforma- 

-t-  V*\ 

tion.  The  forces  acting  on  the  i —  disc  are  shown  in  Figure  5-r'- 
The  axial  force,  T7Y.  ,  is  given  by 


T7Y.  =  XK.(ACi  -  AC0.)  +  C.  AC 


(  l;i  ; 


The  shear  force,  T7Xi3  is  given  by 


T7Xi  =  XKs.(BCl  -  BCol)  +  CS1  BCj 


(15) 


The  restoring  torque  is  given  by 


>1  ■  XKT1  ((ei  -  6ol>  -  (ei-l  -  W 


+  ctl<8l  -  h-l> 


(  lh  ) 


where 


=  stiffness  of  the  axial  spring 
=  stiffness  of  the  shear  spring 
=  damping  in  axial  loading 
=  damping  in  shear  loading 
=  stiffness  of  restoring  torque  spring 
=  damping  of  restoring  torque  spring 


i- a  i  ■  -i  iiiBii 


We  have  to  note  that  the  disc  has  not  been  modelled  as 
a  beam  although  the  above  method  is  very  close  to  one.  The 
disc  has  been  replaced  by  a  system  of  springs  that  behave  in 
the  above  manner. 

The  reaction  of  the  i—  and  the  i  +  i—  discs  on  the  i— 
vertebral  body  are  shown  in  Figure  5-7  where  T6Y^,  ToX^  and 
Bi+1  are  the  reactions  of  the  i  +  i—  disc  on  the  i—  vertebral 
body  and  T7X^,  T7Y^  and  are  the  reactions  of  the  i—  also 
on  the  i —  vertebral  body.  Mathematically  they  are  given  by 
the  following  equations: 


T6Yi  XKi+l<Aci+i  “  ACoi+l^  +  Ci+]  ACi+] 


T6Xi  XKsi+l(BCi+l  "  BCoi+l)  +  Csi+1  BCi+l 


Bi+i  XKTi+l{(0i+l  °oi+l)  _  (Gi  "  0oi)} 

+  Cti+I^ei+1  -  V 


and  T7Y.  ,  T7X^  and  are  from  Equations  (1^1)  -  (1(3). 


5 •  Forces  Developed  at  the  Facets 

The  articular  facets  have  been  modelled  by  two  springs, 
one  limiting  rotation  and  the  other  limiting  the  sliding  of 
one  vertebra  over  the  adjacent  ones. 

The  force  resisting  relative  rotation  between  the  i  — 
and  the  i-1—  vertebrae  in  the  erect  mode  is  given  by 


t  ==  xKhi  x  h1  x  sin  {(e._1  -  eoi_1)  -  (0.  -  eoi)}  (20) 


n:t,  intr  at  a  distance  h^  from  the  center  of  the  vertebral  body 
In  a  direction  parallel  to  the  longitudinal  axis  of  the 
vertebral  body,  where 


hi  “  the  distance  of  the  articular  facets  from 
the  center  of  the  vertebral  body 
XKh^  =  stiffness  of  the  spring  resisting  the  relative 
rotation 


xn  the  hyperex o ended  mode,  due  to  the  forced  change  in  curva¬ 
ture  of  the  spine,  it  is  assumed  that  the  facets  have  "bottom' 
out  and  hence  the  lamina  acts  as  another  beam  parallel  to  thM 
dioC.  The  change  in  length  A^A^  is  used  as  the  measure  of 
axial  deformation  of  this  beam,  and  the  force  developed  is 
given  by 
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=  XKh , 


( A1A2 l 


AlA20i^ 


A1A2  =  lenSth  at  time  t  >  0 
A1A20  =  len§th  at  time  t  =  0 


(20 


The  force  resisting  sliding  motion  at  the  facets  between 
the  i—  and  the  i-1—  vertebrae  is  given  by 

Fxl  =  “hi  *  <V3i  -  A2A30i > 

Due  to  the  overlapping  nature  of  the  articular  facets  it  is 
difficult  to  define  the  point  of  application  of  the  above 
iorce.  However,  in  the  model,  it  has  been  assumed  that  this 


. .  ■ . . . . .  .  *•- 
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force  acts  perpendicular  to  the  longitudinal  axis  of  the 
vertebral  body  at  a  point  (d^  +  AC^/2)  below  the  center  of 
the  vertebral  body.  The  reaction  of  this  force  on  the 
vertebra  immediately  below  acts  at  a  distance  d.  ,  above 
the  center  of  the  i-1—  vertebra.  This  assumption  is  justi¬ 
fied  because  the  superior  articular  facets  are  shorter  than 
the  lamina  and  the  inferior  articular  facets. 

With  the  above  forces  developed  in  the  spine  we  can  now 

t  h 

draw  a  free  body  diagram  of  the  i—  vertebra,  (Fig.  5.8). 

This  diagram  does  not  include  the  seat  back  reaction  force, 
shoulder  strap  or  lap  belt  forces  and  the  chin-chest  contact 
force  because  these  forces  act  only  on  selected  vertebrae  and 
are  described  in  the  section  on  auxiliary  equations. 

6 .  Auxiliary  Equations 

(i)  Shoulder  strap  force: 

Tne  exact  model  of  the  shoulder  strap  is  not  possible 
since  r<o  data  is  available  regarding  the  distribution  of 
inis  force  via  the  rib  cage  to  the  spine.  The  mode]  has  the 
capability  of  distributing  this  force  to  various  levels  of 
the  spine,  however  in  the  present  study,  the  development  of 
the  shoulder  strap  force  is  based  on  the  movement  of  Ti  and 
the  entire  force  is  transmitted  to  that  vertebra.  The  shoulder 
strap  force  in  the  model  is  also  assumed  to  exert  a  restoring 
moment  on  Tl  to  restrict  its  rotation.  This  is  a  valid 


assumption  due  to  the  friction  generated  between  the  belt 
and  the  shoulders.  The  shoulder  strap  is  considered  to  have 
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a  soft  spring  for  the  first  three  inches  of  deformation, 
after  which  the  stiffness  of  the  spring  increases  six  times. 

This  allows  for  the  initial  compression  of  soft  tissues  and 
also  for  an  initial  "hunching"  described  by  Vulcan  [21]  and 
observed  in  high  speed  movies  of  the  experiments. 

The  forces  developed  in  the  shoulder  strap  are  given  by 
the  following  equations: 

The  initial  length  of  the  shoulder  strap  is  given  by 
^io  =  uoi 

where  uQi  is  the  distance  of  the  center  of  T^  from  the  seat 
back. 

The  length  of  the  strap  at  time  t  >  0  is  given  by 

Z .  =  ./(  w  .  -  w  . ) ^  +  u. 2 

i  v  oi  i  i 

where  vi  =  initial  height  of  T1  from  the  seat  pan 

and  w^  =  height  of  T1  from  the  seat  pan  at  t  >  0 . 

The  force  generated  in  the  strap  is  now  given  by 

T4  =  XKSi(5,i  -  JL  )  (2 

This  force  is  set  to  zero  if  SL .  <  l  .  or  if  u.  <  u  .  signifying 

l  oi  i  oi  &  J  F;> 

that  T1  has  moved  towards  the  seat  back. 
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(ii)  Lap  belt  force: 

The  lap  belt  force  is  generated  when  the  pelvis  ,  the 
first  rigid  link  in  the  model,  moves  away  from  the  seat  back. 
One  end  of  the  lap  belt  is  assumed  to  be  fixed  at  the  inter¬ 
section  of  the  seat  back  and  the  seat  pan  while  the  other  end 
is  attached  to  the  center  of  the  first  link.  Hence,  the  lap 
belt  generates  a  vertical  and  horizontal  force  at  the  pelvis 
given  by 


where 


and 


strapu  =  XKS^  x  (^u^2  +  w^2 

2  2  -1 
Ulo  +  wlo  ^  *  cos(tan  w1/u1) 


'ay} 


strapu 

XKS, 


un 


w, 


lo 


lo 


horizontal  component  of  lap  belt  force 
stiffness  of  lap  belt 

distance  of  the  center  of  pelvis  from  seat  back 

distance  of  the  center  of  pelvis  from  seat  pan 

distance  of  the  center  of  pelvis  from  seat 

back  at  t  =  0 

distance  of  the  center  of  pelvis  from  seat 
pan  at  t  =  0 


2  x  2 

•1  +  W1 


strapw  =  XKS^  x 

I  2  2  -] 

-  >/ui0  +  wlo  )  *  sin(+  an  ’  w]L/u1) 


where  strapw  =  vertical  component  of  lap  belt  force. 
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(ill)  Seat  back  reaction: 

The  seat  back  reaction  force  Is  generated  when  any  jink 
makes  contact  with  the  seat  back.  From  the  x-rays  taken  before 
the  experimental  runs  we  can  determine  the  vertebrae  that  are 
in  contact  with  the  seat  back.  Hence,  if  during  the  run, 
these  vertebrae  tend  to  move  towards  the  seat  back,  a  reaction 
force  is  generated  on  the  vertebrae  involved.  For  the  vertebrae 
not  in  contact  with  the  seat  back  initially,  distances  are 
specified  for  motion  towards  the  seat  back  without  any  seat 
back  reaction.  If  the  movement  of  a  vertebra  is  beyond  the 
distance  specified,  a  seat  back  reaction  force  is  generated  to 
stop  the  vertebra.  Hence, 


XSF  =  XKSB  x  (h .  +  Ffcp  -  u.  )  u,  <  h.  +  F£p 

1  1  L  1 


=  0 


u.  >  h.  +  F£.p 

ii 


(2'l) 


where 

XSF  =  seat  back  .reaction  force  on  the  i—  vertebra 
h^  =  distance  of  articular  facet  from  the  center 
of  the  vertebra 

F£p  =  distance  allowed  to  move  before  contacting 
the  seat  back 

XKSB  =  stiffness  of  the  seat  back 
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(Iv)  Chin-chest  contact  force: 

The  rotation  of  the  head  relative  to  the  torso  is  impeded 

when  the  chin  contacts  the  chest.  A  chin-chest  contact  force 

is  generated  to  stop  further  rotation  of  the  head  relative  to 

the  torso.  The  relative  angle  between  the  head  and  the  torso 

for  the  chin-chest  contact  varies  from  one  spine  to  another 

and  is  determined  from  movies.  In  the  model,  the  criterion 

for  chin-chest  contact  was  chosen  to  be  relative  angle  between 

T1  and  the  head.  Hence,  if  6^  is  the  angle  moved  by  T1  and 

0  is  the  angle  moved  by  the  head,  the  relative  angular  move- 
n 

ment  of  the  head  and  T1  is  given  by 
Rdisp  =  9h  -  0T1 

and  the  chin-chest  contact  force  is  given  by 

HC  =  Ch1  x  (Rdisp  -  0CH1 )  (26 

Ch1  =  spring  constant  Ibs/rad 

OCh^  =  relative  angular  movement  between  head  and 
chest  without  contact 

Ch .  is  considered  to  be  constant  for  .2  radians  and  is 
increased  three  times  for  relative  rotation  greater  than  .2 
radians.  Hence,  essentially  we  have  a  bilinear  spring,  the 
initial  part  being  soft  to  simulate  the  soft  tissue  compres¬ 
sion.  It  should  be  mentioned  here  that  no  experimental 
determination  of  the  chin-chest  contact  force  has  been  made, 
hence  the  constant  Ch^  is  selected  arbitrarily  to  stop  the 


where 

and 
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head  In  the  model  based  on  the  observed  angular  rotations 
by  Vulcan  [21]. 

The  reaction  force,  HC,  is  assumed  to  act  at  a  distance 
1.5"  from  the  center  of  gravity  of  the  head  parallel  to  the 
S-I  axis  of  the  head.  An  equal  and  opposite  reaction  force 
is  exerted  on  Tl. 

(v)  Reaction  from  the  hyperextension  block: 

The  hyperextension  block  is  placed  opposite  LI  and 
covers  T12  and  L2 .  Hence  in  the  hyperextension  mode  a  seat 
back  reaction  is  generated  at  T12  ,  LI  and  L2  similar  to  that 
given  in  section  (iii).  The  stiffness  of  the  block  has  been 
assumed  to  be  the  same  as  that  of  the  seat  back. 

7  .  Equations  of  Motion 

1“  In 

A  free  body  diagram  of  the  i—  vertebrae  is  shown  in 
Figure  5-8-  For  the  sake  of  clarity,  the  auxiliary  forces 
described  in  section  6  have  been  omitted  but  they  have  been 
utilized  in  the  computer  program. 

t  h 

Resolving  the  forces  on  the  i— —  vertebra  parallel  to  ir 
and  axis,  we  get: 

ZFui  =  (T6Xi+pX2-FXl)  cos  0i  -  (T6Y. +T52-T51 )  sin  6. 
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ZFwi  =  ( -6Yi+T52-T51)  cos  6;.  +  (T6Xi  +  FX2-PXl)  sin  6. 

-  T7Y1  cos{e1  -  (eol-0ol_1) > 

-T7Xl  sinie,  -  (eol-eol_l)}  (27) 


Taking  the  sum  of  moments  about  the  center  of  gravity, 


EKGi 


T7Yi  cos{0ol-eol_1)  e,  +  -  TBY.-e, 

-  T6XJ.d1  -  T7Y,  sin(e0l-90l_1)  d, 

-  T7Xj  cos(eol-eol_1)  d. 

-  T7X1  sinO^-e^^)  ei  -  B. 

+  (T51-T52)(ei+hi)  -  FX2-di 


-  FXl(d1+ACi/2) 


(28) 


The  auxiliary  forces  and  moments  are  added  to  the  above 
if  applicable  to  the  vertebra  in  consideration. 

Now,  acceleration  of  the  center  of  gravity  is  given  by 

3-g  =  X  I  +  Y  J  +  I  +  J  +  e^  0^  cos  6^  J 

•  2 

-  0^  cos  0^  I  -  e^  e  sin  6^  I 

'  2 

-  0^  sin  0^  J 

*'  *  p 

=  {X  +  u^  -  e^  0^  cos  0^  -  e^  0^  sin  0^}  I 

+  {Y  +  wi  +  ei  0^  cos  0^  -  e1  0^2  sin  0^  J 

(29) 

where  I  and  J  are  unit  vectors  parallel  to  the  ui  and  w.^  axes, 
and  X  and  Y  are  the  horizontal  and  vertical  acceleration  of  the 


sled. 
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Using  Newton's  laws  of  motion: 


‘  2 


EFui  =  mi^X  +  ui  ~  ei  9i  oos  ej  “  ei  9j_  sin  e1}  (3 


EF  .  =  m  (Y  +  w.  +  e.  0.  cos  0.  -  e.  0. 2  sin  0.}  (3 

W1  1  1  1  1  1X1  X 


mGi  IG  9i 


where  m .  is  t^e  mass  supported  by  the  i—  body  and  I,,  is  th» 

1  U 

polar  moment  of  inertia  about  the  center  of  gravity. 


8 .  Solution  of  the  Equations  of  Motion 

Equations  30-32  are  non-linear,  second  order  differential 
equations.  For  each  body  we  have  three  such  equations,  hence 
there  are  a  total  of  J8  equations  to  be  solved  simultaneously 
for  i  =  1-26.  This  is  achieved  by  using  an  IBM  supplied 
numerical  technique  on  an  IBM  36 0  digital  computer.  The 
method  used  is  known  as  Hamming's  Predictor  Corrector  method. 

To  use  the  above  routine,  we  have  to  reduce  each  second  order 
differential  equation  to  two  first  order  differential  equations. 
Hence,  we  now  have  a  total  of  156  first  order  differential 
equations  to  be  solved  simultaneously  with  given  initial  con¬ 
ditions,  i.e.,  uQi,  wQi,  60l,  uQi,  wQi,  0Ql.  The  accelerations 
x  and  y  are  input  parameters  corresponding  to  the  sled  accelera¬ 
tion  in  the  x  and  y  directions. 

To  integrate  numerically  using  the  above  method,  we  have 
to  have  the  highest  order  derivative  on  the  left  hand  side 
of  the  equation.  Hence,  we  have  from  Equations  (30-32) 


130 


9i  “  ^Gi^Gi 


EF  .  .  2 

u.  =  — —  +  e.  8.  cos  6.  H-  e .  6.  sin  6.  -  X 

l  m.  li  ill  l 


EF  .  .  ? 

w.  =  — —  +  e.  0. 11  sin  0.  -  e.  0.  cos  0.  -  Y 

i  ni .  li  l  ii  l 


9 .  Choice  of  Parameters 

(i)  Geometry  of  the  spine: 

The  model  requires  the  initial  coordinates  of  the  center 
of  each  vertebral  body,  the  length,  the  angle  the  longitudinal 
axis  of  the  vertebra  makes  with  the  vertical,  the  distance  of 
the  articular  facets  from  the  center  of  the  vertebral  body 
and  the  thickness  of  the  disc  coincident  with  each  vertebral 
body.  These  measurements  were  made  with  the  help  of  x-rays 

taken  before  each  run  with  the  cadaver  placed  in  the  accelerate 

It  should  be  noted  that  the  trapezoidal  shape  of  the  vertebral 

bodies  is  an  idealization  in  the  model,  since  in  many  spines 

the  intervertebral  discs  are  wedge  shaped,  especially  at  Lrj. 
However,  the  curvature  of  the  spine  has  been  simulated  by 
making  the  rigid  links  trapezoidal. 

(ii)  Mass  and  moment  of  inertia  of  each  rigid  link: 

The  model  requires  the  distribution  of  the  mass  of  the 
whole  body  to  the  rigid  links.  Liu  et  al.  [10]  have  estimated 
the  mass  distribution  of  segmented  cadaveric  trunks.  The 
mass  distribution  used  in  the  model  is  consistent  with  their 
data.  The  three  different  cadavers  used  were  weighed  before 


the  experiment.  After  the  experiments  the  torso  was  cut  at 
the  level  of  L3  and  the  two  halves  were  weighed  again  to  giv*3 
the  total  mass  above  L3-  Three  further  segments  were  made  of 
the  upper  torso  -  the  first  consisting  of  L3  to  T8,  the  second 
from  T8  to  C7  and  the  third  from  C7  to  the  head.  The  three 
segments  were  weighed  again,  their  centers  of  gravity  were 
determined  using  a  load  platform  which  utilizes  three  load- 
cells  at  fixed  distances  from  one  another.  A  trifilar 
pendulum  was  used  to  determine  the  moment  of  inertia  of 
these  segments.  Based  on  the  above  measured  data  the  mass 
and  moment  of  inertia  at  the  various  vertebral  levels  were 
estimated.  The  mass  of  the  arms  were  assumed  to  be  distribute 
to  the  first  five  thoracic  vertebrae. 

(iii)  Physical  constants: 

Each  rigid  link  in  the  model  is  associated  with  a  disc 
requiring  three  spring  constants  and  facets  requiring  two 
spring  constants.  lienee,  a  total  of  130  spring  constants  are 
required  for  the  model.  Theoretically  all  1 30  constants  can 
be  different  and  in  an  ideal  situation  they  have  to  be  ex¬ 
perimentally  determined  for  each  specimen.  Very  Tittle  data 
is  available  in  the  literature,  and  in  the  case  of  the 
articular  facets  nothing  is  available.  Vulcan  [71]  has 
measured  the  stiffness  in  compression  of  the  lumbar  and  the 
lower  thoracic  vertebrae.  Markolf  [ll]  has  measured  the 
axial  stiffness  of  discs  from  L^t  to  T8  and  has  found  values 
ranging  from  7000  Ibs/in  in  the  lumbar  region  to  19000  Ibs/in 
in  the  thoracic  region.  No  data  is  available  for  the  upper 


thoracic  level. 


Markolf  also  gives  the  rotational  stiffness 
of  discs  from  700  in-lb/rad  to  2400  in-lb/rad.  However,  thes**1 
measurements  were  made  with  no  pre-load  and  at  very  small 
deflections .  Most  biological  materials  possess  a  sirmoidal 
load-do  flection  characteristic  -  i.e.,  the  stiffness  increases 
with  increasing  loads.  Hence,  the  above  values  for  rotational 
st  iff n'.  ..s  of  t  he  discs  appear  to  be  much  too  low  for  our 
application.  Hence ,  in  the  lumbar  region  a  rotational  stiff¬ 
ness  of  6000  in-lb/rad  was  assumed.  In  the  thoracic  region 
the  rib  cage  Imparts  more  rigidity  to  the  spine,  hence  a 
rotational  stiffness  of  12000  in-lb /rad  was  assumed.  In  the 
cervical  region  a  rotational  stiffness  of  2400  in-lb/rad  was 
assumed  based  on  the  values  given  by  Vulcan. 

The  constants  used  in  the  model  are  shown  in  the  Appendix 
or.  input  data  to  the  model. 
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Free  Body  Diagram  of  the  i—  Vertebral  Body 


CHAPTER  VI 


RESULTS  OF  THE  MATHEMATICAL  MODEL 
AND  COMPARISON  WITH  EXPERIMENTAL  RESULTS 

Three  cadavers  were  use a  to  experimentally  verify  the 
mathematics  !.  .  Each  cadaver  was  run  at  6  ,  8  and  10  g's 

in  the  erect  and  hyperoxtcndod  modes.  Figures  6.1  through  6.3 
show  the  results  of  the  experimental  aid  model  runs  on 
Cadaver  No.  2209  in  the  erect  mode.  The  model  predicts  the 
first  ana  second  peaks  of  the  intervertebral  ax  it  1  force  -.1 
[■3  at  6  g's.  At  8  and  10  g's  the  first  peak  in  the  axial 
force  is  lower  than  that  measured  experimental Jy .  In  the 
experiments  we  can  see  tint  considerable  amount  of  tensile 
force  is  generated  at  the  *'acets  at  the  time  of  the  first 
peak.  However,  in  the  model  enough  relative  rotation  between 
L3  and  L9  has  not  taken  place  to  generate  tension  at  the  facets 
at  the  time  of  the  first  peak.  Hence,  the  model  predicts  a 
lower  first  peak  in  the  I n'-'-r vertebral  load  at  L3 .  It  should 
be  noted  here  that  no  attempt  was  made  to  curve-fit  the  ex¬ 
perimental  data  with  the  model  at  the  three  acceleration  levels. 

Figures  6.9  through  6.6  show  a  comparison  of  the  experi¬ 
mental  results  and  the  model  results  at  6,  8  and  10  g's  in  the 
h.yperex tended  mode  for  Cadaver  No.  2209.  Both  the  model  and 
ine  experiments  indicate  a  reduction  of  the  intervertebral  load 
at  L3  due  to  the  facets  carrying  compressive  loads  for  a 
longer  duration  when  compared  with  the  erect  mode.  There 


is  rood  correlation  between  the  experiments  and  the  model 
results.  The  second  peak  in  the  model  appears  earlier  than 
that  in  the  experiment.  This  indicates  that  tne  stiffness 


con. .-t ants  chosen  in  the  mode. I  resulted  in  a  higher  frequency 
content  when  compared  with  the  cadaver  used. 

fi  cures  (.-.7  if  rough  6.9  show  a  comparison  of  the  model 
result.,  v/ith  the  experiment al  results  in  the  erect  mode  of  runs 
waac  on  Cadaver  No.  2231-  Trr  resuits  are  similar  to  those 
obtained  on  Cadaver  Mo.  2209-  The  correlation  between  the 
xporimental  and  the  model  results  is  very  pood  in  this  iever. 
An  :nt  ere  sting  point  tc>  is  the  increase  of  only  29  Jhs. 

in  the  peak  intervertebral  load  measured  when  changing  the 
ae  ••deration  lu’-el  from  P  ’•  to  in  g's  whereas  the  change 
from  6  g's  to  8  g's  is  200  1  bs .  The  model  predicts  a  change 
of  292  lbs.  rron  6  so  8  g's  and  a  change  of  only  75  lbs.  from 
8  to  10  g's.  Hence,  it  can  be  said  that  the  model  has  simu¬ 
lated  the  three  runs  very  accurately. 

The  experimental  and  the  theoretical  results  of  the  runs 
made  in  the  hyperextended  mode  are  shown  in  Figures  6.10 
through  6.12.  There  is  good  correlation  between  tne  model 
and  the  experimental  intervertebral  force,  but  the  facet  load 
predicted  by  the  model  at  the  first  peak  is  higher  than  the 
experimental  facet  force.  This  is  because  the  model  has  a 
Larger  dynamic  overshoot  at  the  first  peak  when  compared 


with  the  cadaver. 

The  experimental  and  model  results  for  Cadaver  No.  2 N 1 3 
are  shown  in  Figures  6 . 1 3  through  6.18.  There  is  a  good 
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corn0  1  at  ion  for  the  6.5  and  8  g  runs  between  the  model  and  the 
experiments.  However,  during  :;he  10  g  experimental  runs  in 
the  erect  and  hyperext ended  modes,  the  head  of  the  cadaver  hit 
Ui'  rails  vn  the  sled ,  '  s  a.  resuit  of  which  the  .  ‘‘~l  leak  in 
the  axial  force  did  not  bu.iii  up  ns  predicted  by  the  model. 

In  r.he  hype rox:  d  m. :d-;  Uv-  model  nas  a  higher  frequency 

content  t  han  the  ob.-. -rved  results . 

Co:. cl  us  Lons_: 

Over  the  range  of  measurements  made,  there  is  a  good 
>.  -rrelut  ion  between  the  experimental  and  the  theoretical 

:lts .  Due  to  th»j  '  i  of  .i  nst  rumer.t  at  ion  for  the  experi¬ 
ments,  attempts  were  not  mad*  to  verify  the  model  at  di  f ferent 
■.  :  »  -bra l  D5’ ho .•! ; >/  he  in-vivo  measurement,  of  axial 
force  has  been  made  f '  r  *  : :  -  *  first  time.  With  more  refinements 
in  ’he  experimental  •  •••  ?.h:>  imms  ,  measurements  at  various 
vertebral  levels  can  be  made  to  verify  the  mathematical  model 
further.  Due  to  a  lack  of  data  on  the  various  stiffness 
v a  'ues  needed  by  the  m  i-  ?  ,  a  perfect  match  between  the  ob¬ 
served  and  the  predicted  time  history  of  the  axial  force  has 
not  been  achieved.  However,  the  model  is  general  enough  to 
taa:  any  constants  if  the  data  is  available.  Some  optimiza¬ 
tion  techniques  are  available  for  curve  fitting  and  "hunting" 
for  constants.  These  techniques  if  used  in  conjunction  with 
tr  >:  model  may  lead  to  better  understanding  of  the  physical 
constants  involved. 

The  model  does  bring  out  the  importance  of  the  initial 
curvatures  of  the  spine  as  can  be  seen  by  the  difference  in 


I  Hi 

the  response  of  the  spine  between  the  erect  and  the  hyper- 
'.•x tended  modes. 

The  effect  of  the  seat  back,  restraint  system  and  the 
fti'i  -  '.■!)  °  s  t  contact  fore'-  -  n  be  studied  with  t  V  ]  r.  of  the 
model.  The  results  of  ..he  xu'ia .  :  _n  of  the  above  parameters 
have  not  been  Inc .  .  !  Ir;  this  study  because  the-  above 

pare  me'-  -rs  were  n  "l  ■viii^d  Jn  the  experiments,  but  in  si.mu- 
1 ai  ion  runs  made  with  the  model  these  parameters  do  alter 
the  response  of  the  spine.  lienee  in  future  experimental 
work,  the  above  parameters  should  also  be  studied. 


Cadaver  2209  in  the  Erect  Mode 
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Fig.  6.14.  Comparison  of  the  Model  and  the  Experimental  Results  of  8g  Run  on 
Cadaver  2413  in  the  Erect  Mode 
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.15.  Comparison  of  the  Model  and  the  Experimental  Results  of  lOg  Run  on 
Cadaver  2-413  in  the  Erect  Mode 


CHAPTER  VII 


CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

On  the  basis  of  experimental  and  theoretical  data  it  is 

now  possible  to  explain  the  mechanism  of  anterior  wedging 

fracture  of  the  vertebral  bodies  during  +  G  acceleration. 

z 

A  combination  of  bending  and  transference  of  load  from  the 
articular  facets  to  the  vertebral  body  results  in  high  com¬ 
pressive  strains  at  the  anterior  surface  of  the  vertebral 
body  resulting  in  fracture  when  the  structural  limits  are 
exceeded. 

Qualitative  and  quantitative  evidence  has  been  presented 
to  document  the  load  bearing  capabilities  of  the  articular 
facets.  This  is  the  first  time  that  such  a  study  has  been 
done.  The  discovery  of  a  load  path  through  the  articular 
facets  has  led  to  the  fabrication  of  a  device  to  increase  the 
threshold  of  fracture  of  the  spine  during  +  G^  acceleration. 

This  is  of  practical  interest  because  a  recent  report  on  in¬ 
juries  sustained  by  pilots  during  the  Viet  Nam  war  states  that 
almost  one-third  of  the  pow's  sustained  injuries  involving 
fractures  of  the  spine  during  emergency  ejections  from  aircraft. 

A  mathematical  model  has  been  formulated  and  experimentally 
verified  to  predict  the  response  of  the  vertebral  column 
during  +  G  acceleration.  This  is  the  most  comprehensive, 
experimentally  verified  2-dirnensional  mathematical  model  of 
the  spine  at  the  present  time. 


Future  work  in  this  area  should  involve  the  determination 
of  the  material  properties  of  the  intervertebral  discs  and 
vertebral  bodies  and  the  failure  limits  under  a  state  of 
combined  axial,  shear  and  bending  loads.  The  properties  of 
the  lamina  and  articular  facets  should  be  determined.  There 
is  no  work  available  regarding  the  relative  movement  of  the 
laminae,  i.e.  ,  when  and  under  what  conditions  do  they  "bottom 

The  transference  of  the  shoulder  strap  load  to  the  spine 
via  the  rib  cage  should  be  investigated.  In  the  present 
mathematical  model,  all  the  strap  load  is  transmitted  to  Tl, 
but  this  may  not  be  the  case  because  the  rib  cage  would 
transmit  part  of  the  load  to  various  vertebral  levels. 

In  the  present  study  the  effect  of  the  generation  of 
intra-abdominal  pressure  has  not  been  studied.  It  has  been 
hypothesized  by  Bartelink  [1]  and  Morris  et  al .  [12]  that 
this  pressure  can  reduce  the  reaction  upon  the  vertebral  end 
olate  at  the  SI  and  L5  junction.  However,  in  the  present 
study,  the  abdomen  was  eviscerated  in  all  the  cadavers  used, 
hence  there  was  no  question  of  any  intra-abdominal  pressure 
build  up.  But  in  the  case  of  living  human  subjects  in  an 
impact  situation,  there  might  be  a  build  up  of  pressure, 
hence  opening  up  a  new  load  path.  The  question  that  has  to 
be  answered  is  whether  in  an  impact  lasting  200  msecs,  any 
pressure  will  be  generated  and  if  so,  how  much?  The  force 
generated  will  tend  to  counteract  the  forward  flexion  of  the 
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TP  I'l 
1 
2 

1 
4 
6 
6 
7 
R 
7 
in 
1  1 
12 
1  1 
1  n 
i  s 
16 
17 
1  H 

1  q 

17.25 

27 

21 

2? 

2  1 
24 
26 

27 

2R 

2  R 
10 
3 1 
12 

3  3 
14 

35 

36 

37 
1R 
37 
40 

4  1 
47 
4  I 
44 
4  5 
4  f' 

4/ 

4  R 

4  7 
50 

5  1 
52 
5  1 

54 

55 

56 
67 


C  CALCULATE  INITIAL  LENGTHS 

IMPLICIT  R  E  AL*R  (A  -  H  ,  D-  Z ) 

DIMENSION  Cl  0(2  6)  ,  WO  (26)  ,00  (26)  ,0(26)  ,  H  (2  6)  ,  E  (2  6)  ,  XX  H  ( 25  >  ,  XKT  (  26) 
1CT(26)  ,C  (26  )  ,XK  (2  6)  ,  ABO  (2  6)  , ALP  HAD  (26)  ,A  BDO  (2  6)  ,  Pi  10  (26) 
mi  ENSION  0  (26)  ,H  (26)  ,0  (2  6)  ,XM  (26)  ,  X  I  (2  6 )  ,  XK  S  ( 25 )  ,  XD  1  (26)  ,X  02  (2' 
DIMENSION  X  (156)  ,  DX  ( 1  56  )  ,  PR  MT  (5  )  ,  AD  X  (  1 6,  1  56)  ,  7  DT  ( 2  6)  ,  HD  T(  23  )  , 
ilODT  (26)  ,UDDT(26)  ,WDDT(26)  ,ODDT(26)  ,ALDT  (26)  ,ACO  (26)  ,  BCO(  2  5| 
DIMENSION  DSOO  (26)  ,  DC  CO  (2  6)  ,  H  SOO(26)  ,RC33  (26) 

DIMENSION  A1  A23  (26)  ,A  2A  30  (26) 

D  IMENSION  A  A  (  2)  ,SF(26) 

EQUIVALENCE  (X  (  1)  ,'J  { 1 )  )  ,  (X  (27)  ,  TJDT  (1  )  )  ,  (X{5  3)  ,  W  ( 1 )  )  ,  (X  (  77)  ,  W  D  T 
3(1)),{X(105),O(1)),(X(111),ODT(1)) 

CALL  PLn  MX  (H6. 0) 

N  AM  ELIS  T/ONE  /U 

N  AM  ELIST/TWO/W 

NAMELIS^/THRF  E/0 

NA  MELI ST/FO  MR/D,H  ,  E,XM,XI 

NA MFLIST/F I VE/X KR  ,XKT,CT,  C,  XK  ,  XKS 

N  AM  EL  IST/SIX/PA  CO  ,SCH  1 ,  SCH  2  ,S  CH  1 ,  A  A  ,  TT1  ,  TT  2 ,  IC  1  ,M0  DE 
2  M0DE=1  SIGNIFIES  ERECT  MODE 

NA  MELIST/SE  VFN/0CH1  ,0  CH2  ,  C H  1  ,CH  2,  XKSB  , S F 
EXTERNA  L  DE  R IV  ,  OUTP 

COM  MON/CD  IS  T/D,B,E,  ABC,  ABDO,  ALPHAO,  PHI3  ,03  ,73  ,73  ,  X  0  1,  X  02,  kZ  0,  BCi 

COMMON  /CD  NS  T/  XK  H ,  XKT  ,CT  ,C  ,X  KS  ,X  K,  X  M,  XI 

COMMON/FACE  T/Al  A20,  A2  A30 

COMMON  KOU  NT , LL 

COM  M3N/P  ELI/A BS  BO 

COM  MON /PUL  S  E/PA  CC  ,SC  Ml  ,SC  132  ,S  CH  1,  A  A  ,  TT1  ,  TT2  , 1  Cl  ,M  3  DE 
2  OM MON/HE  AD  /OCT  1 , 0CH2  ,CH1  ,CH2  ,X  KSB,SF 
P  EAD(5,0NE) 

READ  (5,  TWO) 

PEA D( 5,  THREE) 

READ  (5,  FOUR) 

R  EAD  (5 ,  El  VE  ) 

RFhD( 5 ,S  T  X) 

PEA  D(  5,  SEVEN) 

LL- 0 

?OUNr  =  0 

PI=  3.  VI  157  265 

DO  117  1=  1  ,  Ifl 

:  (I)=.  8*DSQRT(XK  (I)  *  XM  (I)  ) 

CT  (  I)  =  . 08  *0  SORT  (XKT  (I  )  +XI  (I) ) 

117  CON  TIN  ME 

DO  2  1=1, 26 
n  (I)=  (PI/  18  0. )  *0(1) 

30(1)  =0  (I) 

M  0(1)  =0(1) 

2  U0(I)=W  (r  ) 

DO  111  1= 1 , 76 
MDT(I)  =3. 

WDT(I)  =0. 

od r  (i)  =o. 

1  1  1  CONTINUE 

nATA  PRMT/3.30,0.  20,0.0001,  .001/ 

DO  112  1=1 ,  156 
1V>  DX  (I)  =1  ./1S6. 

DO  110  1=1,  26 
SOD  =  DS  IN  (00  (I)  ) 


SB  COO=DCOS(OD  (I)  > 

sb  c  soo  (i )  =d  (i )  *soo 

SB  DCOO  (I )  =1  (T  )*COO 

6  1  H  SOO(I|  =H  ( I)  *5  )0 

6  2  HCOO(I)  =3  (I )  *  COO 

SB  110  CONTINUE 

64  AB3  (1)  =  W)  ( 1 ) -  D  {  1) 

55  ALPHAO(  1)  =P  1/2. 

56  A  BO  0  ( 1)  =0  . 

57  PHT0{1)=9. 

6  B  X01(1)=0. 

SO  X  02  ( 1 1  =  ABO  (  1) 

70  ACO  (1)  =ABO(  1) 

71  RCO  ( 1)  =0. 

72  AlA2O(1)  =  0. 

71  A2A3O{1)=0. 

7N  A  B5B0=DSQRT  ((10(1)  *  UO  (1)  +  W0  (1)*W0(1)  ) 

7S  00  1  1=2,26 

7 S  X01  (I)=-(m  (I-  1 )  -  DS  00  (I  -1  ))•*•(  IJO  (I )  ♦  PSOO  (I)) 

77  X02  (E)  =  {HO(I)  -0  (I)  *  DC  OS  (00  (T)  )  )-  (WO  (1-1  )  4-  DCOO  ( I  -  1 1  ) 

7  B  ARO  (I)=DSQPT  (XO  1  (I)  *X01  (I)  +X02(I)  *  X9  2  (T  )  ) 

7  0  AL?  HAO  (I)  =DATA  N2  (X0  2  (T)  ,X01  (I)  ) 

BO  13  Y  01  =-(  U0(  I)  -HCOO  (I)  )  +  ('1 0  ( I  -  1)  -HCOO  (I  -1 )  ) 

8  2  Y  02=  (MO  (I )  -HSOO  (I) )  -  (  MO  (I-  1)-  HS  00  (I-  1 )  ) 

BB  .2  5  AE3D0  (I)=DSQRT  (Y01  *Y01  f  Y02  *Y02  ) 

8 B .  5  P  HIO(  I)  =DATAN  2  (Y02,-Y  01 ) 

8  4  A  1A20  (I)  =  ABDD  (T)  *  DS  IN  (PHIO  (I)  -OC  (T  -  1)  ) 

PS  A  2A  30  (I )  =  AR  DO  ( I )  *DC05  (  PHI  0  (I)  -00  (I-  1)  ) 

86  14  ACO  (I)  =ABO  (I)  *DS  IN  (A  LPHAO(T)  -00(1-1)  ) 

8  7  BCO(I)  =  ABO  (  I)  *DCOS  (ALnHAO  (T)-OO(I-I)  ) 

88  1  CONTINUE 

80  CALL  DHPCG  (PRUT  ,X,DX, 156,  THLF,DERTV,OUTP  ,AIJX) 

00  PRINT  123  , 1  HLF 

01  120  FORMATC  END  OF  SIMULATION.  I  HLF=  *,13) 

02  CALL  PLOT  (0,  0,  099) 

OB  STOP 

04  END 

85  SUBROUTINE  PER  I  V  (T  ,  X  ,  DX  ) 

06  IMPLICIT  REAL*8 (A-H,0-Z) 

07  DIMENSION  X  (156),  DX  (156)  ,  XI  (26)  ,Y 2(26)  ,A  I.DT  (26),  AC  (26),3C(26)  , 

OR  D  ACDT(  2  6)  ,  BC  DT  (26) 

00  DIMENSION  DELU(26) 

00  C  CALCULATE  R.  H.  S.  (S  UB  .  DERI  V) 

0  1  DIM  EN  SION  (IDT  (2.6)  ,V  DT  (2  6)  ,  CDT  (26)  ,  UDDT  (  26)  ,W  DDT  (  26)  ,  OD)  T  (26  )  , 

0  2  «ABDT(26),  AB(26)  ,ALP  HA  (26)  ,  PHT  (26) ,  ABD  (26)  ,  PH 5  (  26  ,2) 

OB  a),  U  (26)  ,  W  (26)  ,0  (26),  RX  (2  6)  ,RY(  26)  ,X10T  (2  6)  ,X2DT  (26),  SICE  U{  26) 

04  Dll  ENSIOY  SO  (26  ), CO  (26)  ,  D  SO  (2  6)  ,  DO  (26)  ,  HSO  (26  )  ,  i3C  0  (  26)  f  3  C  (  2  6) 

05  DIMENSION  A  1  A2  (  26  )  ,  A  ?A  B  (2  6) 

06  DIMENSION  B ETA (  26)  ,  SHTRO  (26  )  ,  SMTP  ( 26)  ,S H  RAP  (2  6) 

07  COM  MON/CD  IS  T/D  (26) ,  It  (26)  ,  F(26  )  ,  ABO  (25)  ,  A  BDO  ( 26)  ,  A  L  PH  AO  ( 2  6)  , 

08  i)P  HIO(2  6)  ,00  (26)  ,U0  (  26  )  ,  MO  (26)  ,  X 01  (2  6)  ,  X 02  (26)  ,  ACO  (  26  )  ,BCO  (2  6) 

00  C  OM  MON/CO  NS  T/XK  H  (2  6)  ,XKT  (26)  ,CT  (26)  ,C  (26)  ,XK5(26)  ,  XK(25)  , 

10  a)XM  (26)  ,XT  (2  6) 

1  1  COM  MON  /DIST/AB  ,ABD,ABDT,  ALPHA  ,  PHI ,  X 1  ,  X2  ,  X  IDT  ,  X2DT  ,  A  LDT  ,  AC ,  BC , 

12  i>  ACDT,  BCDT,S  1GFU 

1  i  COMMON/FACF.  T/A1  A20{26)  ,A2AB0  (26) 

14  COM  MON/ FACET? /A  1  A2,  A2  AB 

15  COMMON/PEL 1/ABS BO 

IS  C OM MON /PEL  VIS/S TRA” 

17  C  0MM0N/S3  0L/SHR  AP 
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l(n2 


1003 


176.26 
1  76  .6 


COM  HON/HE  AD/OCH  1 ,  OCH  2,  CHI  ,  CH2  ,  XKSR  ,  S  F  {  26 ) 

DO  101  1=1,26 
7  (I)=X  (I) 

ODT  (I)  =  X(I+26) 

D  X  { I )  =  0  D1  (I) 

W  (T  )  =X(  1+52) 

HOT  (I)=X  (I*  78) 

DX  (I+52)  =  WDT(I) 

0  {I )  =  X  (  13  4) 

ODT  (I)  =  X(I  +  130) 

DX  (1+104)  =ODT ( I ) 

TF(T.  EQ.  0  . )  GO  TO  1Q01 
GO  TO  107  3 
I  DO  1002  1=1  ,156 
DX  (I)  =0. 

CONTTN0  E 
RETURN 

PI=3. 14159265 
DO  110  1=1,26 
SO  ( I)  =  DSI  N  (0  (I )  ) 

CO  (I)  =  DCO  S  (0  (I)  ) 

DSO(I)  =D  (I)  *SD  (I) 

DC 0  (I)  =D  (I)  *CO  (I) 

HSO  (I)  =  H  (I)  *6  0  (I) 

HCO  (I)  =  H(T)  *CO  (I) 

Z ONTINUE 
A  1A  2  ( 1)  =0  . 

A2A3  (1  )  =0  . 

AB(1)  =  W  (1  )  -  D (  1) 

ALPHA  (1)=PI/2. 

ABD  (1)  =0. 

P  HI  (1 )  =P  1/2. 

X 1 DT( 1 )  =0. 

X  2DT  (1)  =♦  W  D  T  ( 1 ) 

A  BDT  (1 ) =  +WDT  (1) 

AC(  1)  =AB  (1) 

BC{ 1 )  =0. 

A  LD  T  ( 1 }  =0  . 

A CDT  ( 1 )  =  W DT  (1) 

DCDT  ( 1)  =0  . 

XI  (1)  =0. 

X  2  (  1)  =AB  ( 1 ) 

ABSB=D3QRT  (U  (1)  *0(1)  +  W(  1)  *W  (1  )  ) 

S  AN  G=  DATA  N2  (W  (1)  ,U  (1)  ) 

A  BS  T=  ( W  (1  )  *  WDT  (1)  +U  (1)*  'JDT  (1)  )  /  ABS  B 
STR AP=0. 

I  F  (  ABSB.  LE.  ABSRO)STRA  F=0. 

STR  APU  =  -STR  AP*  DCO  S  (SA  NG) 

STPAPW  =  -STRAP*  OSIN  (SANG) 
no  2  1=  2,  26 

XI  (  T)  =-  (0(1-1)  -DS3  (T-  1)  )♦  (0  (I)  ♦DSO(I)  ) 

X2(I)  =  (W  (T)  -DCO  (I))  -  (  W  f  £-  1)  +  DCO  (I-  1 )  ) 

A  B  ( I)  =DSQRr  (XI  (I)  *X1  (I)  +  X2  ( I)  *X2(T)  ) 

BETA(I)  =DATA  N  2  (  (WO  (I)  -W  (I  )),U  (T  ) ) 

SHTRO  (I)  =  DSQR  T (  UO  (I)  *  UO  (I)  ) 

S  HTR  (I )  =  0SQRr  (  (WO  (I )  -  W(  T)  )  *  (  WO  ( I)  -  W  (I)  )  +11  (I)  *0  (I )  ) 
SHDT=  (0  (I)  *UDT  (T)  -W  (I)  *WDT(I)  )  /SRTR  (I) 

SHR  AP  (I)  =XKS(T)  *  (SHTP  (T)  -  S  HTR  C  (I)  ) 

IF  (SHTR  (I  )  -  SB  TRO  ( I)  .GF.2  .)  SHR  An  (I)=XKS  (I)  *?..  +  6.  *X  K  S  (T ) 
ft*  (SHTR  (I)  -S  HTR  0  (I)  -  2.  ) 
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mm 
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1  76 

177 

178 

179 

180 
1  82 

183 

18  3 

184 
1  88 
1  86 
187 
1  88 
1  89 
190 

19  1 
192 

19  3 
1  94 
1  98 

196 

197 

1  98 
199 

20  0 
201 
202 
203 

2  04 
208 
206 

207 

208 
2  09 
210 
211 
2  12 
2  1  3 
214 
218 
216 
216 

217 

218 
2  18. 

21  9 
2  20 
221 
2  22 
222 

22  3 
2  24 
228 
226 
2  27 
228 
22  0 
2  30 

2  JO 

2  3  3 


r  P  (SHRAP(I)  ..  LE.  0.)  S  HR  A  P  v  <  )  =  0. 

IP  (0  (I).  LE.30  (I)  )  SHPAP  (T)  =0. 

T  F(  SHT3  (I).  LF.  S  HT  RO  (I )  )  SH  RAP  (I)  =0. 

ALPHA(I)  =DATAN2  (X2(I)  ,X1  (I)  ) 

22  Y  1=  (0  (1-1  )  -HCO(I-I)  )  -  (0  (I)  -  BCO  (I)  ) 

Y2=  (W  (I)-HSO(I)  )-(W(I-1)  -  HSO  (  I-  1)  ) 

21  PHI  (I)  =DATA  N2  (Y  2,-Y  1) 

.25  ABO  (I)  =  DSQRT(Y1*Y1+Y2*Y2) 

C  HATE  OF  CHANGE  OF  AB 

24  Y  IDT  (I)  =-  (CJ D T  (I  -1)  -D  (1-1)  *0  DT  ( I-  1)  *CD  (I-  1 )  )  ♦  (HOT  (I )  +0  (t) 

3*0DT(  I|  *CO  (  I)  ) 

X2DT(I)  =MWDT(I)  ♦D(T)  *n!)T  (I)*SO  (T))-  (»nr  (1-1)  -0(1-1)  *Dnr(I-  1) 
1*SO  (T-  1)) 

A  RDT  (I )  =+  DCOS  (A  LPHA  (I))  *Y1DT(I)  +  DSIN  (ALPHA  (I)  )*X2DT  (I) 

AC  (I)  =  A8(I)  *  DS  IN  (ALPHA  (I)  -0  (1-1  )) 

BC  (I)  =  AB  (I)  ♦DCOS  (  ALP  HA(  I)  -0  (I-  1  )  ) 

A  LOT  (I)  =  (DO  OS  (ALPH  A  (T  ) )  ♦.<  2DT(  I)  -DSIN  (A  L  ?  H  A  (I )  )  *  X  IDT  ( I)  )  /\ B  ( l ) 
ACDT(I)  =ABDT(I)  *DSIN  (ALPHA  (I)  -0  (T-1) )  +  I3C(I)  *  (  AL  DT  ( I) -ODP  (I-  1)  ) 

3CDT  (I)  =ABDT(I)  *DCOS(  ALPHA  (T)  -0  (I- 1 )  ) -AC  (I )  *  (  AI.DT  (I )  -3  DT  ( I-  1)  ) 

A  1 A  2  (I)  =A  BD  (I )  ♦DSIN  (PHI  (I )  -0  (1-1)  ) 

A 2  A 3  (  I|  =A  BD  ( I)  ♦  DC  OS  (PHT  (I  )  -0  (I  -1)  ) 

2  CONTINUE 

DO  3000  r  =  1,26 
3000  H  C  ( I)  =3  . 

TDE  3P  =  DA3S(0  (13  )-0D  (1  8)  ) 

RD IS P=D  AB S ( 0 (26 ) )  - TDI Sn 

IF(RD1'SP.  GR.OCH  1)  HC{26)  =C  HI  ♦  (RDISP-OCH  1 ) 

I  F(  RDISP.  GE.0CH2)  HC  (26)  CH  1*  (0CH2-DCH1 )  +  CH2*  (RDISP 
3)-OCH2) 

HC  (18)=H:  (26) 

H  C  ( 17 )  =  .  3  ♦  HC  (26 ) 

H  C( 1 6)  =0.  *HC ( 26  ) 

DO  30  01  1  =1  1,15 
300  1  HC (I)  =  0 ♦HC  (26) 

C  C  AL  COL  ATIOS  DP  BENDING  'lOflF  NT 
DO  100  1=1, 2b 
F  LP=0 . 

T  F  ( I.  ED.  5)PLP=3 .9081  49 
I  F  ( I.  EQ.  6  )  FLT=  3.  7  64956 
r  F  (  T  .  EQ.7)  PLP  =  3. 70520  2 
IF  (I.  B  0  - 1  5)  FLP  =  4 .39696/ 

.25  IF  ( I.  EQ.  16)PLP=  4. 504419 

I  F  ( I.  EQ.  1 )  FLP=  3. 

I  F  ( I  .  EQ  .1 )  GO  TD  1 

25  ORE  L1=  (0 (T  -  1)  -00(1-1)  )-  (0  (I)-OO(T)) 

T  1=  XKT  (I)  *ORE  LI  +  CT(  T)  *  (  OD  T  (I- 1 )  -D  DT  (  T)  ) 

GO  TO  3 

1  T1 -XKT(I)  ♦  (  -0  ( I )  +00  (I))  +CT  (T)  ♦  (-OPT  (I)  ) 

3  IF  (I. EQ. 26)  GO  TO  4 

.25  ORE  L2=  (0  (I)  -00  (I)  )  -  (O  (T  +  1  )  -00  (1+1)  ) 

T2=XK  T(  1+  1)  *DR  EL2+CT  (1  +  1  )  *  (ODT  (I) 

ODT  (1+  1)  ) 

GO  TO  5 
'4  T2=  0 . 

5  I  P(U(I)  .GE.  (H(I) +FLP)  )  XSF  =  0. 

TF  (0  (I)  .LE.H  (I)  +  FLP)  XSF=  XKSB*  (H  (T)  +FLP-0  (T  )  ) 

I  F(  I.  GE.  1  8)  X  SF  =  0. 

T3=  0. 

X  K3  T  =  X KS ( I) 

IF  (U  (I).  LE.OO(I)  )  XKS(I)  =0  . 
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2  34 
?  38 
?  38 
2  3  7 
2  38 
2  34 
2  4  4 

241 

242 
24  3 
244 

248  C 

2  46 
247 
14  8 
244 
284 
281 
2  82 
2  83 
284 

2  88  C 

286 

287 

288 
284 
260 
761 
262 

26  3 
264 
268 
266 
267 
7  68 
264 
770 

27  1 

27?  C 

27  3 
2/4 
2  78 

276 

277 

27  8 
274 
284 

28  1 
2  8? 

28  3 
284 
288 
2  86 
287 

287 .28 

2  8  8  c 

2  84 
244 
24  1 
24  •> 


T14  =  SHR  AP  (I)  *DCOS  (BETA  (I)  ) 

T141  =  SH  RAP  (I)  *DSIN  (BPTA  (T)  ) 

T  4=  T  1  4 
T41 =0. 

IP(I.GT.IR)  T4  =  0. 

TP  (I. ST. 18)  T4  1  =  0 . 

IP  (I.  LE.  8)  T4  =  0  . 

I F  ( I .  L  E.  8)  T4  1=  0. 

T42=SHR  AP  (I)  *2.  5 
IF  (I.3T.1  8)  T4 2  =  0 . 

I F (  I.  LE .  3  )  T  42=0  . 

ABOVE  IS  STRAP  FORCE 
X  KHI=XK8  (I) 

T  81  =  X  KH  ff  )*  (A  1A  20  (I)  -  A1  A2  (I)  ) 

T  F(MDDE.  EO.  1)  T  8  1  =-  X  KH  (T  )  *H  (I)  *D  SIN  (OR  ELI) 

I  F  (I .  EQ  .2  6)  GO  TO  6 

T  52  =  X  KH  (T  +  1  )*  (A  1A20(T  +  1)  -A1A2  (1+1)) 

I  F(MO  DE.  EO.  1)  T  52=~  X  KH  (1+1  )  *H  (If  1)  *D5IN  (OR  FL2) 

SO  TO  7 

6  T  82  =  0. 

7  T5=  T51-T  82 

THE  ABOVE  IS  THE  REACTION  AT  FACETS 
I  P(I  .EO  .2  6)  GO  TO  8 

T  6Y  =XK  (If  1)  *(  AC  (If  1 )  -  ACO  ( If  1)  )  fC  (Ifl )  *ACDT  (If  1) 

T6X  =  3  F  (If  1)  *  (XK  (I-M  )  *  (RC(If  1)  -BCO  (T+1)  )  fC  (If  1)*BCDT  (If  1)  ) 
GO  TO  9 

8  P  6Y  =0. 

T6X=0. 

4  '”7Y Y=  XK  (I)  +  (AC  (I)-ACO(I)  )  fC(I)  *ACDT(I) 

T  F  ( I  .  EQ  .  1 )  GO  TO  10 

T  7X  X  =SP  (T  )  *  (  XK  ( I)  *(BC  (I)-RCO  (I)  )  f  C  (I)*BCDT  (I)  ) 

T7Y=r7YY*D:  OS  (00(1)  -00(1-1)  )  -T7XX+DS  IN  (00  (I)  -  00  (I  -  1  )  ) 

I  7X  =T7XX*DC0S  (00  (I)  -00  (I  -  1)  )f  T7YV  *BSIN(00  (I)  -0  0  (I  -  1)  ) 

GO  TO  11 

19  T7Y=r  7Y Y*DC  OS  (0  (I)  ) 

IF  {  AC  (1  ).  GE.  A  CO  ( 1)  )  T7Y  =  0. 

T7X  =T7Y  Y*  DSTN  (0  (T)  ) 

IF  (AG  (1  )  .GE  .  ACO  (1 )  )  T7  X=0. 

THE  ABOVE  APE  FORCES  ON  VERTERRAF  DUE  TO  F>FF03  MA  V  T  ON 
11  SIGNG1  =T1  -T2-T3-T5*  (H  (I )  f  F  (I)  )  f  T4  2 
IF(I.E0.26)  GO  TO  2  3 
S  IG MG2=-  T  6Y  *E  (I)  -T6  X*  D(  I) 

GO  TO  2  5 

23  GTG  MS2  =  HC (I)  *  1. 5 
25  SIGH33  =  T7Y*E(I)  -T7X*n(I) 

TF  ( I.  EQ.  1 )  S  TS  f13  .3=0. 

PX1  =XKH  (I)  *  (A2A  3  (I)  -  A  2A  30  (I) ) 

T  F  ( I .  E  Q  .2  b)  GO  TO  7  0 

p  X2=  X  KH  (Tf1)*(A2A3(If1)  -A2  A30  (Ifl)) 

GO  TO  7  1 

70  FX2  =0 . 

71  STS  HG4=PX  1*  (0(1)  f  AC(I)/2.  )  fFX2*D  (I) 

G  IGMG=S  IGHG  If  SIGHG  2f  SIGHG3-SIGMG4 

0  DDT (  I)  =S  IGHG/XT  (I) 

I  F  ( I  .  EQ  .1  )  ODDT(I)  =0. 

SnMHAriON  OF  FORCES  IN  X-DIRECTION 
F  X 1  U=  F  X  1  *  DC  OS  (0(1)  ) 

I  F(I.  E2.26)  GO  TO  31 
FX2  U=F  X2  *  DCOS  (0(1)) 

FX30=fT6X  *CO(I)  fHC(T)  *S0  (18) 
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GO  TO  32 
n  FX2U=0. 

FX3U=0. 

2  P  X4  U=  -  T  7X  *C  0  (I ) 

S  IGFX=  FX2  U+FX30  +  FX4n-FXin 
IF(I.RQ.I) GO  TO  5 1 
S  TG  MS=-T4+X  SF 
GO  TO  52 

51  SIGMS=STR  APU+XS  F 
2  STGPU(I)  =  SI  GFX+ S IGMS 
SUMMATION  OF  FORCES  IN  Y-DTRFCTION 
FY1=T5 

rP(I.  EQ.  2  6)  GO  TO  41 
17  Y2  =T6Y 
GO  TO  4  2 

41  F Y2=  HC  (T) 

42  FY3=T7  Y+HC{  I)  *DCOS  (0  (26)  )  *DC0S(0(  I)  )  +  HC(  I)  *  DS  IN  {0  (I )  )  *0S  I  N  (0  (2  6) 
SIG  Y -FY1 ♦ FY2-FY  3 

IF  (I.EQ  .1 )  T41  =  STR  AP  W 
62  S  IGFY=S  IGY 
EOT)  AT  IONS  OF  MOTION 

UDPT(E)  =  ( 1  .  /XM  ( T)  )  *  (S  TG  FI!  (I)  -STGFY  *SO(I )  )  ♦ 

3F  (I)  *0D  DT  (I  )*S0  (I)  +  F(  I)  *0  DT  (I)  *ODT  (I)  * 

(dCO  (  I)  -  XDDT  (T) 

IF(I.EQ.I)  IF  DPT  ( I)  =  0  . 

0  MOOT  (I)  =  {1  ./XM  (I)  )  *  (STGFY*CO  (I)  +T41  +  (T  6X-T  7K+  FX2- FX1  )  *S0  (I)  )  + 
t)F(I)  *ODT  (I)  *ODT  (I  )  *S0  (I)-E(T)  *  (000^(1)  )  *C0  (I) 

S-YT)DT(T) 

XKS  (I  )  =XKST 
XKH  (I )  =  XKHI 
00  CONTINUE 
00  200  1=  1,26 
nx  (i^  130)  =onoT(i) 

OX ( 1  +  78) =  WODT  (I) 

0  DX(I  +  26)=IIDDT  (I) 

RETURN 
E  NO 

FUN  CT  ION  YD  DT  (T) 

IMPLICIT  R  EA  L*H  (A  -  H  ,  0-Z  ) 

DTM  FN  SI  ON  SC(2) 

G  OM  MON/PUL  S  E/P  A  CC  ,  S  C  HI  ,  SC  H2 , 5CH  3  ,A  A  (2)  ,T  1  ,T2,IC1  ,  MODE 
SC(  1»  =  386 .4  *A  A  (  1) 

SC  (2)  =3  85  ,4*AA(  2) 

IF ( IC  1  .  E 0  •  1  )  G 0  TO  3 
I F  (T .  GE  .T2)  GO  TO  1 
I  F(  T.  GE.  T  1)  GO  TO  2 
YDDT=T*SC  (1)  /T 1 
RETURN 

2  VDriT=3C  (1 )  -  (SC  (1)  -  SC  (2)  )  *  (T-T1)  /(T2-T1) 

RETURN 

1  Y  DDT  =  SC ( 2 ) 

RETURN 

3  YD  DT  =  0 . 

RETURN 

END 

FUNCTION  XDDT(T) 

IMPLICIT  RE  A  L  *  R  ( A  -H  ,  O  -Z ) 

COM  MO N /PULS  E/PA  CC  ,SC  H  1,  SCH  2  ,S  CH  3,  A  i  (2)  ,  T1  ,  T2  ,  IC1  ,MODE 
TP  (I  Cl  .  EQ  .0)  GO  TO  3 

IF  (T.  LE.5CH  1)  XD DT  =-  f  P  ACC/ SC  HI  )  *  T*  33  6.  4 
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154 

155 
I1)*) 

157 
154 
154 
160 
161 
16  2 
36  1 

164 

165 

166 
1  67 
164 
150 
17n 

171 

172 
171 
174 
17  5 
176 

i  n 

377.25 
174 
1  /') 

140 

141 

142 
34  3 
3  04 
3  45 
3  46 
14  / 

3  44 
340 
300 
30  1 
102 
101 

104 

105 
1  06 

3  07 
100 
100 
400 
40  1 

4  0  2 
4  01 
4  04 
40  5 
4^6 
4  07 
4  04 
4  00 

410 

411 
4  12 


I  F{  T.  GT.  5  CHI)  XDDT  =-  PACC  *38  6.  4 

TF  (T.GT.SCH2)  XDDT=-  (PACC-  (  (PACC/  (SCH  3-0CH2)  )  * 

9  (T-  SCH  2) )  )*3B6.4 
RETURN 
1  XDDT=0. 

RETURN 

FND 

SUBROUTINE  OUTP  (T  fX,DX,THLF,  NDIM,  PRMT) 

IMPLICIT  REAL*B  (A-H,0-7.) 

DIMENSION  X  (1  56)  ,  DX  (156)  ,  FY  (26)  rS  MG  (  26)  ,  X  1  (26)  ,  X2  (2  6) 
DIMENSION  FXT  (26)  f  FXB  (2  6)  ,  FYT  (2  6) 

DIMENSION  0  DT  ( 26)  ,  WDT  (2  6)  ,  ODT  (26)  ,  UDDT(  26)  ,WDDT(26)  rODDT  (26  )  , 
SABOT  (26),  AB  (2  6)  ,  ALPHA  (2  6)  ,P  HI  (26)  ,A BD  { 26)  ,  RHS  (2  6, 2  ) 

#,U(  26)  ,»  (26),  0(26),  RX  (2  6)  , RY(26)  ,  XI  DT(25  )  ,X20T  (26)  rSI3FT  (26) 
«,  PRMT{5)  ,  AL  DT  (2  6)  ,AC  (26)  ,  BC  (2  6)  ,  ACDT(2  6)  ,  BCDT  (25  ) 

COM  MON/CD  I  ST/D(  26  )  ,H  (26)  ,  ?  ( 26  )  ,  AB  0  ( 26)  ,  A  PDO  (2  6)  ,  ALPHAO(25>  , 
3PHI0  (26)  ,  00  (2  6)  ,UO(26)  ,  HO  (2  6)  f  XO  1  (26)  ,X0  2  (26)  ,  A  CO  (2  6)  ,  BCO  (16) 
COMMON  /CO  NST/XKH(26)  ,  X  KT  (26),  CT  (26),C(25)  ,XKS(26)  ,XK(2S), 
SXM(26)  ,  XT  (2  6) 

COMMON/DC  ST/AB,  ARD,  ABUT,  ALP  HA  ,P  HI  ,  XI  ,X2  ,X  IDT  ,  X  2D  Tr  A  LD  T,  AC  ,  BC  , 
a  A  CD  T,  BCDT  ,S  IGPU 
COMMON/P  FL  VTS/PLO  AD 
C  OM  MON/FACE  TV /A1  A2  (26)  ,A2A3  (26) 

COM  MO  N/ HEAD /OCH 1 ,  0CH2,  CHI  ,  CH2  ,XKSR,SF(?6) 

COMM0N/FACET/A1A2O  (26),  A2  A10(26) 

COMMON/SHOL/SHRAP  (26) 

REA  L* 4  XP  (2  8)  ,  YP(28) 

DATA  XP  (27)  ,XP(28)  ,  YP(2  7)  ,YP  (24)/3.  ,1.  ,0.  ,1./ 

DO  101  1=1,26 
11(1)  =X(Il 
ODT  (I)  =X  (If  26) 

U  DDT (I) =D  X  (1  +  26 ) 

W(I)  =  X(T*  52) 

WDT(I)  =X(  T*  78 ) 

WOOT  (I)  =DX(If  78) 

0  (I)  =  X  (If  104) 

ODT  (I)  =  X(If  130) 

X  P(  I)  =U ( I )  /3. 

Y  P  (I)  =N  (T  )  /3.0 
101  ODDT  (I)=DX  (If  110) 

P  1=  3.  14  150  265 
COMMON  KO  UN  T,  LL 
KOU  NT  =  KOUNTf  1 

C  DETERMINE  CORRECT  PRINT  INCREMENT 
I  F(T.  LE.  LL*. 00 5) RETURN 
CALL  PLOT  (2  .0,0.0  1) 

CALL  LI  NE  (X  P,  YP,26, 1  ,1  ,0| 

C  CALCULATE  SEAT  PAN  LOAD 

SLC=(  ( ABO  (1)  -  AB  (1)  )  *XK  (1)-C  (1  )*ABDT(1)  )  *DSIN(  AL?  HA  (1)  ) 

IF  (ABO  (1)  .LE.AB(I))  SLC=0. 

C  CALCULATE  SHOULDER  STRAP  LOAD 
STR  AP  =  S  HR  A P  (18) 

DO  100  1=1,26 

T  7Y  =X  K  ( I )  *  (  AC  ( I )  -ACO  (I)  )  fC(I)*ACDT  (I) 

T7  X=SP  (I)  *  (X  K  (I )  *  (BC  (T)  -BCO  (I) )  f  C ( I)  *BCDT  (I)  ) 

I  F(  I.  EQ.2  6)  GO  TO  6 

T  6Y=XK  (If  1)  *(  AC(Tfl)  -  ACO  (1*1)  )  *C  (If  1 )  *A CDT  (If  1 ) 

T6X  =  SF  (If  1)  *(XK  (If  1)  *  (  BC(  I  f  1 )  -BCO  (If  1) )  fC  (If  1 )  *BCDT  (I  fl )  | 

GO  TO  7 
6  T6r=0. 
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n  1 1  ?6  x  =o . 

4  14  7  FY  (I)  =T7Y 

4  IS  FXT  (I)  =T6 X 

4  1f>  FXB(I)=T7X 

4  17  FYT  (T)  =T6  Y 

4  IS  IP(I.EQ.1)GO  TO  3000 

416.76  ORE  LI  =0  (I  -1 )  -00  (1-1)  -  (O  (I)-00  (I)) 

4  VI  T1=XKT(T)  *0REL1+CT(I)  *  (OPT  (1-1) -ODT  (I)) 

4  70  0,0  TO  300  1 

471  70^0  T  1  =  XK  T  ( I)  *  ( -0  ( I)  *00  ( T  )  )  KTT  (I)  * ( -ODT  (I ) ) 

4  7:1  300  1  SMC  (I)  =T1 

4  3'.: .  1  DIMENSION  T  5 1  (2  6) 

4  2?.  2  5  T51  (I)  =  XKH(I)*(A1A2  (I)-A1A20(I)  ) 

4  22.5  I  P  {  MODE  .  E  0.  1)  T51  (I)  =  XX  H  { I )  *  H  ( I)  *DS  I  N  (OPE  Li) 

471  100  CONTINUE 

434  PRINT  6  ,T,SLC, STRAP, PLOAD, 0(26)  , FOUNT 

426  5  FOPMATpX,1  T= '  ,  F5 .  3, IX,  'SEAT  L  CAD  =  •  ,  F6 . 1  ,  1  X,  •  STR  AP  =  •  ,  F5.  1  , 

4  26.25  2 IX,  ,LAP  =  ‘  ,P6.1,  IX,'  HEAD  A  N  GLE=  '  ,F  6.  3  ,  IX  ,  »KOUNT='  ,13) 

4  26.6  PRINT  2  10 

4  26.  6  210  FORMAT  (1  HO,  »L  INK  •  ,  3  X  ,  •  A  X I A  L*  ,  6X  ,  '  SHE  A  R«  ,  5  X,  •  MOM  ENT  1  ,6  X ,  '  FAE  E  r  ' ) 

4  26.7  PRINT  215 

426  .  H  216  FORMAT  (IX,*  NO.',  IX,' FORCE  (L  P)  *,  2X  FORCE  (LR)  ',1X,  •  (I-J-LB)  *  , 

4  1 6  .  H  1  j)4X,  'FORCE  (LB)') 

476  203  PRINT  20  1  ,  ( I,FY  (I)  ,FX  B  (I)  ,  SMG  (I )  ,T5  1  (I)  ,  1=  1 , 26) 

47S  201  FORM  AT(1  X,I3  ,  1X,F10.  4, 1X,p  10.  4,  1X,ri  0.4,  1  X,  FI  0.4) 

4  26.  1  PRT  NT  211 

4  2  4.25  2  1  1  *OR  MAT  (1H0,  'LINK',6X,  '  O',  0X,»  LPT'  ,0  X  ,’  IJDDP  '  ,  OX  ,  '  W  ' ,  9X , 

426.6  2‘  WOT'  ,7  X,  ' W DDT  1 ) 

424.6  PRINT  212 

424.7  2  1  2  FORMAT  (IX,*  NO.  •  ,  4  X  ,  '  ( I  NS)  •  ,  4  X  ,  •  ( I  N/SEC )  »  ,  3X  ,  '  (I  N/S  EC2)  '  , 

424.4  36 X,  •  (INS)  •  ,4X,'  (IN/SEC)  '  ,2X,'  (IN/SEC2)  ') 

470  PRINT  202,  (1,0(1)  ,(JDT  (I),UDDT(I),B(T)  ,WDT(I)  ,W0DT  (I)  , 

4  11  21=1,26) 

4  12  2  n  2  FOR  MAT  (IX, I  3,  5X,F10.4,1X,F10.4,1X,F10.4,1X,F10.4,1X,F10.4 

4  12.  1  2,  IX,  PI  0.4) 

4  12.  76  PRINT  213 

4  12.6  213  FORMAT  (1H  0,  'LINK' ,4  X,  *0  »,  11X,  'ODT' ,BX,'ODDT' ) 

432.6  PRINT  214 

4  3  7.7  2  1  4  FORMAT(1X,'  NO.  •  ,  2X  ,  '  (RA  I) )  '  ,  6  X  ,  '  ( R  A  0/S  EC)  '  ,  2X  , '( R  A  0/S  EE  2)  ' ) 

4  13  PRINT  2  04  ,  (1,0(1)  ,  ODT  (I)  ,  OP  DT  (I)  ,1=1,26) 

4  16  26  4  FOR  MAT  ( IX  ,T  3,  FI  0.  4,  3X  ,  FI  0 .4  ,2  X,F10. 4  ) 

4  17  LL=LL  +1 

4  VI  KOUNT=0 

4  3  0  0  FT  U  P  N 

440  END 

END  OF  F  r  L  E 
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APPENDIX  II 

DESCRIPTION  OP  THE  INPUT  DATA  TO  THE  COMPUTER  PROGRAM 

The  Input  Data  to  the  computer  program  Is  in  the  Namelist 

format.  The  symbols  used  are  as  follows: 

u  =  An  array  of  the  x  coordinates  of  the  centers  of  the 
26  links  at  time  t  =  0. 

Units  =  inches 

w  =  An  array  of  the  y  coordinates  of  the  centers  of  the 
26  links  at  time  t  =  0. 

Units  =  inches 

0  =  An  array  of  the  angles  made  with  the  horizontal  by 

the  26  links  at  time  t  =  0. 

Units  =  degrees 

D  =  An  array  of  the  half-height  of  each  link. 

Units  =  inches 

E  =  An  array  of  the  eccentricity  of  the  center  of  mass  of 
each  link. 

Units  =  inches 

H  =  An  array  of  the  distance  of  the  facets  from  the  center 
of  each  link. 


Units  =  inches 


XM  -  An  array  of  the  mass  of  the  26  links. 

Units  =  lb-sec^/in 

XI  =  An  array  of  the  Moment  of  Inertia  about  the  center 

of  gravity  of  each  link. 

? 

Units  =  lb-in-sec' 

XK  =  An  array  of  axial  stiffness  of  the  26  discs. 

Units  =  lb/in 

XKH  =■  An  array  of  the  stiffness  of  the  facets. 

Units  =  lb/in 

C  =  Axial  damping  array. 

Units  =  lb/in/sec 

XKT  =  Rotational  stiffness  array  of  the  discs. 

Units  =  in-lb/radian 

XKS  =  Strap  stiffness  array. 

Units  =  lb/in 

9CH1  =  Relative  angle  between  head  and  Ti  at  time  of  chin- 
chest  contact. 

Units  =  radians 

0CH2  =  Relative  angle  between  head  and  Tl  (see  eqn.  ). 
Units  =  radians 

CHI  =  Chin-chest  contact  resistance  stiffness. 


Units  =  lb /rad 
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CH2  =  Chin-chest  contact  resistance  stiffness. 

Units  =  lb/rad 

PACC  =  Horizontal  Acceleration  Input  (see  Fig.  A. 2.1). 

SCH1,  SCH2 ,  SCH3  =  Horizontal  Acceleration  Input  (see  Fig.  A. 2.1). 
AA ,  TT1,  TT2  =  Vertical  Acceleration  Input  (see  Fig.  A. 2. 2). 

Id  =  Parameter  controlling  mode  of  Input  Acceleration . 


INPUT  DATA  FOR  CADAVER  2209  ERECT  MODE 


r;  ;>  ;>  o  i 

1  A  ONE  U  =2.  18,  2.  69  0427, 2. 7  92  220  ,  2.  7  10  2  38,2.  07787,  2.360789,  2.  148  SO  <  , 

2  1  .  96  08  76  ,  1.8  0298  1,  1.  PS  12  3  0,1  .91  740  5 , 2.  0  30  6  07 , 2.  1  4  72  88,  2. .  321  4d0, 

1  2  .  SOU  92 ,2.9  37  0  1  3,3.  31  69  6  7,  3 . 72  83  1  9,  4  .  1 2  7  20  5  ,  4.  09  9  11  3 , 4.  8847  79, 

4  0  .  2  4  02  04 ,5  .  3  9  2  8  1  8  ,  0. 8  3  8  1  68  ,  6  .  0  9  3  3  1  8,  6.71  3  071  , 

5  Ap  NO 

*  AT  VO  W=  2.  19  ,  3.  800049,  0 . 391  53  6,  6  .78  908  3  ,9.  18  20  24, 9.  56066  8,  1  0.  34  7  3  5  9  , 

7  12.0  8312  9,  1  3  .1  73  178  ,  14  .  17  07  27,  10.  1  72  79  0,  1  6  .1  6630  1,  17.  19  3  9  7,17.  969307, 

8  18. 8  44269,  1  9.623  80  6,  20  .3  278  30  ,21.0  1380  9,21 .64  8911,  22 .23029  1,72.  7  6  068  6, 

'1  2  1.2499)  3,23.  734238,  24.202048,24.816  177,26.9  26860, 

in  A  SN  0 

11  ATHPRF  0=-1  2  .,  -7  .,  3.  ,4.  ,8.  ,  10 .  ,  n.  ,  8.  ,  1  .0,  -2  .,  -6  .0  ,  -6 . 5  ,-8 .  1  6.  2  ,  -  2 1  .  2 

1  2  -  2  7.  7,  -  30.  2,  -32.  2,  -  32.  ,-  30.  ,-36.  37  .  ,-  34.  ,-  2  1.  ,  -1  8.,  -1  2  ., 

11  A "N  n 

14  A  7  0  0  P 

10  D  =  .8,6*.s,  .  40,  4*. 4 ,3*.  30  ,  3*.  3,  2*  .  2  0,  0*  .2 , 1  .«, 

17  E=  0.  ,  2*1 .,  1 .3,  1  .4,  1 .7,  1 .8  ,  1.9  ,2*2.  ,.  9,4*  1.,  4*1. ,7*0., 

18  8=2*1.  ,3*1.  1,  3*  1.2,3*1.1,1.00,  .  84,  .94, .9, .97, .9, .3, .70, .72,.  5, 

19  3* .  5, .  4,  .  3  5, 

2n  XM  =  .  17  68  6,  2  *.01  2  0,  .0  03  70  3  7  ,.00  3008  8  , .  0  0  3  4 1  70 ,  .  0029  3022,  .  0066)  79, 

71  .9006399  ,.  00  58787,  3*  .0007377,  .  014847  ,.010  10,.  0  16  309,.  01  72  66,  .01  7  87  , 

21.20  3* .00241  82  ,4*.  00  10  1  1  38  , .  0 2 832 5 , 

22  XI  =.  404,  2*.  07,  4*  .0  10  6  3  ,.  0  6  25,.  0  0  34,  .  C04  6,  .  04  91  ,  .  04  76  , .  0  39  2  )  3  43  , 

7  3  .  0  27  8,  .0  205,  .  01  84,  .  006  6  ,  .004  , .  00  3  , .  0  0  3 ,.  0  020 ,  .  001  6  ,  .  001  3 , 

24  .9013, .2816, 

70  A  R  ND 

76  ap TVS 

77  X«>  10  9  0  9  0.  ,  5*1  0  0  0  0.,  4  *12  0  0  0.  ,8*  140  0  0.  ,  1*18  94  .,4  *0  4  2  4  .,  2  2  1  9  0.  , 

2  8.1  X  K  H=  9 .  ,  17*69  00.  ,  7*  1  2  00.  ,  1  7  00  ., 

?9  C=  740,7*  1  0.,  3*10  .,  2  *20  .,  0*20.  ,3*7.  6  ,4*  1.  3,6.  34, 

30  XK  T=  19  00  0.  ,7*6  00  0.  ,  10*  12000.  ,  8*2400., 

31  C  =20.,  7*10., 1  0*20.  ,7*  1.  ,  2., 

12  XKS=  10  00.,  7*0.,  6  *0  .,  7*0.  ,0.  ,20.  ,8*0.  , 

13  A  EN  13 

34  A. SIX 

90  r\CC=9.,  SCH1  =.  01  ,SCH2=  .9  6, SC  H  3  =  .2,  A  A  =6.,  0.,  TT1  =  .03  ,TT2  =  .04  ,10  1  =  3  , 

10.70  MO  DE  =  1 

9  6  A  PN  0 

37  A  SEVEN 

18  Of  HI  =1  .,  OCH2  =1  .2  ,  ,C8  1=  53  0.  ,C  H2=  1  000.  ,X  KSB  =00^  .  , 

3H  .  20  SP  =  7*.  0,  19*1 

3  J  A  S  N  0 

m  or  piLE 


INPUT  DATA  FOR  CADAVER  2209  hYPEREXTENDED  MODE 


I  76 


'f  °  T*7'Oll  M 

1  60  NF  U=  2.  1  ,2.  7  31448,  3.  0^  2689,  3  .  1  09  282,  3 . 0'S  OR  4  3  , 2 . 9  57  5  38  ,  2.  8520  97, 

2  2.72  2  3  1  3,2.6  2  6  4  4  1  ,2. 5  4  9  /  27,2  .  4  9  2  1  6  3,  2.4  3  9  82  6,2  .  3  9  0  1  06  ,  2  .  36  7  4  33  , 

3  2.  40  75  3  3, 2. 507  56  2,  2.656758  ,2.833  289,3.014729,  3.  1  966  86,  3  . 3  9005  i>  , 

4  3.  57  3  7  75,3.  7  64  HO,  3.  91  9375,  4.044083,4.  332906  , 

5  8  END 

h  F.^WO  w=i  .7,  3.942  85  0,5.  36  6  1  7  7 ,  6.  7  5929  6,  8 .  1  577  68,  9 . 6  54  64  3  ,  1  0 . 84  9  9  36  , 

7  12.093080,  13.188893,  14.1  8586  3 ,15.  184  196,  16.  182  81  6,1  7.1  315  00,  13  .)  30  63  5  , 

0  18.9  29260,  19.77259  8,  20.558441 ,21  . 3  3857  7 , 22.  0  6  62  84 , 2  2 . 74  1  98  9,  2  1  . 3  6  2  51  8  , 

9  2  3  .9  3  203  7 ,24.50  2  487,  25.0  81  84  8,  25  .6  6854  9,27.3  49  9  4  5, 

1n  6° NO 

11  f,v  HP  EE  0  =  -1  7.,  -1  2  .  ,-  7.  ,  1  .5  ,4  .  ,  3.  5 ,6.  5,  5.  ,  5.  ,  3. 5,  3  .,  3  .,  3  .,  -1.  ,-5.  ,-  10.  , 

12  -  1  2.  ,-  14.  ,-  1  4.  17.,  -1  8.,  -19  ., -1  6.  ,-13.  10.  ,- 7.  , 

n  5  PND 

14  6°  OUR 

15  0=1.  34,  6*.  5,  .45,  4*.4  ,3  *.  35 , 7*.  3 , 2*  .  25 , 5* .  2 , 1  .  H, 

16  F  =  0.  ,2*1 .  ,  1.  3,  1.  4,  1.  7,  1  .  8,1  .  9,2*2  . ,  .9,4*  1.  ,4*  1.  ,7*0.  , 

17  9=2*1. ,3*1. 1,3*1  .2,3*1.1,1.05,.  34, .94,  .9,  .  97,  .9,  .8, .75, .72  ,.  5, 

18  3*. 5, .4,  .3  5, 

19  X!  =  .  1763  6, 2*.0  12  5, .  0  03  753  7,  .0035  08  R,  .00  34  17  5,.  0  02  93  022,  .  00664  79, 

20  .9  05  6099  ,.  0053  78  7, 3*.  00  57  32  7,  .  01  4  84  7,  .015  15  ,.  3  16  359  ,.  01726  6,  .017  37, 

2  1  3*.  0  02  41  82, 4*.  00  15  1  1  38  ,.928  325, 

2  2  XT  =  .  4  54,  2*.  0  7,  4*  .015  63,  .0625  ,.  05  34  ,.  0546 , .  049  1 ,  .  04  76 ,  .  03  92 ,  .  )  7 4  3  , 

2  3  .  3  27  8  ,  .3  255,.  0  18  4,.  0  06  6,  .  0  C4  ,  .  003  ,  .00 1 ,.  0  0  25  ,.  00  1  6  ,  .  001  3 , 

24  .  001  3  , .2  816  , 

25  F.END 

-,h  r.wjVR 

>9  XK  =  133010.  ,8*8  00  0.  ,4  *1  00  00.,  5*14  00  0.  ,3*18  94.  ,  4*  542  4.  ,  22 1. 90  ., 

28  XKH  =  0. ,5*b033.  , 1  2* 29  00. , 7*  1200 . ,  1200., 

29  24  0,  7*10.,  3*  15.,  2*  20.  ,5*25.  ,1*  2.  6, 4*1.  3,  6.  3  4, 

1  9  XK  T=  13  03  0.  ,7*6  00  0.,  1  0*12  000.  ,8*2  40  0.  , 

31  CT=2  9.  ,7*10.  ,10*20.  ,  7*  1.  ,  2., 

32  X«S=1O00  ,,  7*0.  ,6*0.  ,2*0.  ,0.  ,  20.  ,8*0.  , 

3  3  F,END 

3  4  R  9  T  X  , 

3  5  PACC  =9  .  ,  SCH1=.01,SCH2=.06,SDH3  =  .2,AA=8.,7.5,TT1=.34R,rr  2- .  058,131  =0, 

3  5.  2  5  MODE  =0 

3  h  F,  8  N3 

3  7  F,  8  EV  EN 

3  8  0CH1  =  1. ,  3CH2=  1.  2  ,,CH  1  =500.  ,3132  =1  500.  ,XKS  P=  159.  , 

39  5  F  =7  *  .  5,  19*1  ., 

4  9  F,  E  ND 

..)  0?  ptl  F 
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INPUT  DATA  FOR  CADAVER  #2231  ERECT  MODE 


L  22  n 

1  89  NF  1=4.2,  5.361  970,  5.43  4295,5.36454  8,5.  17320  3,  4.  951617,  4. 311  328, 

2  4.65  9.310  ,4. 54  40  5 ,4.  4  69  26  5,  4.  42172  8,4  .4  34  802,4. 52  0178 ,4.  69  6826,  4.  97572  3, 

1  5.  2Q  77  39, 5. 600764, 5. 856  19  4,6.  1 0  2  1  0  3,  6.  4  03  2  75,  6  .82  1  03  7  ,7  .  2  23  9  3  3  , 

4  7.  640467,8.075334,  8.472024  ,9.5  16  200, 

5  fiVNO 

8  5  TWO  W=3  .6  ,6.905665,8.  451958,  10. 049521 ,11  .6366,13.  16  9  3 1  3,  1  4.  6  62  6  84  , 

7  1  6. 154  93  P,  1  7.524  88  7, 18  .7  47  58  9,  14.  846283,  2  0.895416,21.39  1312, 

3  2?.  87319  9,  2  3  .8  33  16, 24.  81  1508  ,25.7  17  148, 26.  48053,  27.21  5469,  27. 335555, 

(>  2  3.4  27  17  ,  23.  99  59  87,2  9.561  5  84,3  0.  1  104  53  ,30.  68  2  907,  3  2.7  32  2  03, 

10  8EN0 

11  Fi’HRFT  O  =  -20.,-5.  ,1.  ,5.  ,  10.  ,  5.  ,  6.  ,  5.5,  2*3  .5,1  .,-3. 5, -7.  ,-15.  ,-18.  ,4*- 18. 

12  -  3  7.  ,-  36.,- 3  5.  ,  -38.,  -3  9.  ,-27.  ,-  27.  , 

1  3  6 54n 

14  9  POOR 

15  n  =  2.  36,  .  55,1  *.6  ,  3  *.5  5,.  47  5,2*.  45,4*  .4,  .4  5,  .  3  2  5  ,  .  3  , .  27  5  ,  6*  .  25 , 1 . 8  , 

16  ?=0.  ,2*1  .,  1  .  3,  1 . 4,  1 .7,1 .  3  ,1.  9  ,2*2.  9 , 4  *  1.  ,  4  *  1. ,  7*0. , 

17  8  =  2*  1.  ,3*1.1  ,3*1. 2, 3*1.1,  1.05,.  84,  .94  ,  .9,  .97  ,.9, .3  ,.  75,.  72, .  5, 

18  3*. 5,  .  4,  .3  5, 

14  H  =  .  19  72  2,  .012  5,  .  0 1  2 5, 6* . 004 40  86  ,4  *.  0  10  3 1 ,  5*.  01  6562  5,  3*  .0025  37  , 

14.25  4*  .00  159  59  ,.  02991, 

20  XI  =.4  54  ,  2*.  0  7  ,  4*.  0  1563  , .  0625,  .  0534  ,  .  0546  , .04  9  1, .04  76, .0392, .0340, 

21  .  0  27  8,  .025  5,  .  0  1 3  4  ,  .  0  06  6  , .  0  0  4  ,.  00  3,  .  003,  .  0025,  .0016  ,.001  3  , 

??  .  001  3,  .2  816, 

2  3  8  END 

24  8  17 1 V  8 

25  x;<  =  1  00000.,  2  *8  000.  ,5  *3000.  ,4*  1  031  0. ,  3*16562  3  *2537.  ,4*1595.  ,199  9  1., 

26  X KH=  0.  ,  5*6000.,  1  2  *20  00.  ,7*1200.  ,  120  0., 

.7  7  C=  24  0,7*  10.  ,  3*  15.  ,  2*20. ,  5*25.,  3*2. 6,  4*1  .3,6. 34  , 

2  8  XFT=0., 7*6000.  ,1  0*  12000.  ,8*2 40 0.  , 

20  7T=20.  ,  7*10.  ,1  0*20.,  7*  1  .  ,2.  , 

in  XK  S=  10  00  .  ,7*0.  ,6*0. , 2*0. ,  0., 20. ,8*0.  , 

31  f-FND 

2  4  OS  IX 

35  P\c:  =  0.  ,  s:3  1=.  01  ,  SCH2=.0  6,SCH3  =  .  2,  A  A  =6.  4,6.  ,TT1=.  03,  TT2  =.04  ,  IC  1=0  , 

35  .2  5  MODE=1 

3  6  6  E  NS 

17  83EVEN 

33  07  F31  =  1. 2,  O: 3  2=  1. 4,  G1  1  =30  0.  ,C 82=900.  ,XKSB=  150.  , 

3't  S  T?=7  *.5,  19*  1.  , 

4  0  SEND 

TO  0"  FILE 


I  80 


n  ^.M1!l 

1 


10 
1  1 
32 
{  3 
3  4 
IS 

IS  .  25 


INPUT  DATA  FOR  CADAVER  #2231  HYPEREXTENDED  MODE 


I  A 
3  / 

38 

19 

»  OF  FTC 


‘‘.ONE  tj=4 .2  ,5.  3  75  862  ,5.50  2109,  5.  519  561, 5.4  777  07,8  .3  88755  ,5.  3268  1b  , 

5.  2100  63  ,  5.  0  26  07  8,4  . 786  098  ,4 . 5 67  54  3 , 4.  49  2  61  8,  4 .4  5772  6,  4 .4  35055 ,4.  49780  S, 
'4  .64  3  1  46  ,4.  8  3854  0,4.99903  0,  5.1  53  54  0,5.  37  699  R  ,  5.  766  094 , 6 . 1  8 1  767, 

6. 53  76  9^  ,5.3  53  54  4,7.  110828,7.821566, 

67  ND 

r.r  wo  *  =  3 .6,  6  .890363,  8.433189 ,10. 03  30  36, 1  1.  631  574,  1  3.1  78802,  U  .67  7  463  , 
16.17205  8,  17. 5  34  48  5,  IB.  733  978,  1  9.8  08456,20.8  54  82  6,21.853  973,  2  2.852600  , 

2  3  .85  0  02  1, 24 .367798 ,25.  803604,  2  6  . 5  92  43  8 , 27  .  35  18  68 , 28.  0  3  1  60 1 , 2  8  .  fa  1  1  588  , 
29.1  73  294, 2  9.7756  35, 30 . 40  1  337,  31 .050186,  3  3.23761  0, 

87  NO 

FT  H°  Ep  9=-  2  2.  ,  -7  .  ,  -1 .  ,  0. ,  4  .,2  . ,  i  .,7  .  ,9  .  ,  15.  3,6.  ,  1 .  ,  3 . 5,  -2  .,  -5  .,  -  1  1  . 5  , 

-1  1  .5, -1  1.5  1  1.  5,-31.  ,-  39. ,-32.,  -28.,  -24.  ,-  18.  ,-  18.  , 

6  7  NO 
f.^on  r 

[)=  2.  3b  56,  1*.  6,  3*. 5  5,  .4  75, 2*. 45, 4*. 4, .45,.  325,. 3, .275, 6*. 25, 1.3, 

E-0.  ,2*1.  ,1.3,  1.  4,  1.  7,  1.  8,  1.  9,  2*2.,  .9,4*1  .,4*1  ,  ,7*0.  , 

8=  2*  1.  ,  3*1.1,  3  *1  .2, 3*1  . 1,1. 05, .84, .94, .9, .  97,. 9,  .8,  .75, .72, .5, 

1* .5 ,. 4,  .  15, 

X1=.  19722  ,.0  125,.  0  125,  6*.  00440 86, 4 *.01  03  1,5*.  0  165626,3*.  002537, 

4*. 001  5959,  .  0299  1  , 

Xr=.  454,  2*.  07,  4*  .01563,  .06  25  ,.  05  34  ,.  0  546,.  04  91,  .04  76,  .0192,-0  340  , 

.0  27  3  ,.  )  255,.  0  1  84,  .  0  06  6,  .004,  .  00 3  ,  .003  ,.  00  25  , .  00  1 6, .  001  1, 

.001  3, .28  16  , 

57  NO 
IV  F 

XX  =1  00  00  0.  ,  2*8  00  0.  ,  6*8  000.,  4*!  0  3  10  .,  5*  16  56  2.  ,  3*  2  53  7.,  4  *1  5  95  .,2  999  1  .  , 

XK  8  =  0. ,5  *6  000.  ,  12*  20  00.  ,  7*12  00  .,  12  00., 

"=  24  0,  7*  10.,  3*  15.,  2*  20.,  5*  25.,  3*  2.  6,4*  1.  3,  6.  3  4, 

XK  T=  0  .  ,7  *60  0  0.  ,  1  0*12  00  0.,  8*2  aOO.  , 

CT1  =2  0.,  7  *  10.  ,  10*  20.  ,7*  1.  ,2  ., 

XKS=  1  0  00  .,  7*0.,  6  *0.,  2*0.  ,0.  ,20.  ,8*  0., 

FEND 
5  SIX 

P  A  CC  =  9  .  ,  SC  HI  =.  01  ,  SC  H  2=  .  0  6 ,  SCH  3= .  2  ,  AA  =1  0.  ,  1  0 .  ,  TT1=.06  5,rr2=.  065 , 1  C  1  =0 
nAC7  =  9. ,  SCH1  =.01  ,  SCH2=.7  6  ,SCH3=.  2,AA=6.4,6.,  TT1=.03  ,TT2=  .  0  4 , 13  1=  0, 

“  0  D  F  =  0 
5  1EV7N 

0^81  =  1 .2  ,0  CH2=  1.  4,CH  1=  300.  ,CH  2  =900  .  ,  XK  SB  =150  .  , 

S'7  -7 *  .  5,  10*1  ., 

F7NP 
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INPUT  DAI  A  FOR  CADAVER  #24  I  3  ERECT  MODE' 

;>  m  * 

1  AO  NR  !l  -  4  .  2,  4  .6  57  29  3,  4.  71  91  91,4  .644216,4.  4  398  10,4.  1  17877,  3.871  401  , 

’  1. 72  8191  ,3.  7 16  86  6,  3.  82 6682, 3.9954  87,  4. 32 C004 ,4.750  177, 5.  21 18 3D , 5.  657  0  9 

3  0.  13  04  37,6.5954  3  5,6. 99  49  85,7.  27  5577,  7.  4  793  3  8,  7  .64  7  78  3  ,7  .  783  4  52, 

4  7.  87 6  9  7b,  7. 937099,  7 . 960649  ,7.960649, 

8  5  RN  n 

0  5  I'm  W  =  3. 6, 7.  1  79  60  8,  8. 87  69  50,  1  0.  574.315,1  2  .25  04  4  3,  1  3.  08  2  444,1  5.4  0  02  60, 

I  1  6.95021  5,  1  8 .1  95732,  19.  4519  20, 2  0.55574  0,  21  .62  5565,22.655  15  1 ,  2  3.  58  5  38  8 , 

0  24.  4  23  61  5,  25.2  24  09  1  ,  26. 029480,26.78  4  180, 27.  47871  4,  2  8.1  4  814  8, 23. 3  27  332  , 

"  29.514404,30.2  07  489,  30.90  4663,  31  .6  04  C34, 3  3  . 90  43  22, 

T  AR  49 

11  or  HIRE  0=-8.  ,  -4.,  1  .  ,  5.  ,1  0.  ,1  1.  5  ,9.  ,  3.  5 ,  -  5.  ,  -  5.  ,  -1  4  . ,  -2  1  .  , -2  5  .  ,  -  26  .  , 

1  2  -11.  39  .  ,  -  30.  ,  -  24. ,  -1  8.,  -15  .,  -12.  ,-9.  ,-6.  3.  ,0.  ,  0., 

11  A  r*J  D 

14  A* OH  3 

16  19=  L  36,4*.  65,.  6,  .55 , 2  *  .5 , 4  *.  46  ,4*.  37  5,.  3, 7*.  2  5,  1  .8, 

16  p-0.  ,2*9.  ,6*0.  ,.  9,.  7,.  65,  .5,  .2  ,4  *0  .,7*.  5  ,  1.  , 

17  FI  =2*1. ,3*1.1  ,3*1.2  ,3*1.  1,1.05,.  P4,.c4,  .9,  .97,  .9,. 8  , .75  , .72  ,.  7, 

10  3*.  7,.  7,  .  7, 

10  XI  =  .  2147  ,.  0  125,.  0125,  5*.  0  0493  3  ,5*  .  007764 ,5*.  0  17  36,  2*.  00  2  1 , 2  *  .  0  03  5  , 

00  .  0016, 2*.  001  3, .0  26  1  1  , 

21  XT  =  .  4  54,  2*.  0  7,  4*  .  0  1  56  3  , .  06  25  0  6  34  ,  .  0  546,  .  04  9 1  ,  .  04  76  ,  .  0  3  9  2  ,.  )  1 43  , 

23  .  0  2  7  0  ,  .3  2  5  5,.  3  1  84,  .  0  0  6  6,  .  0  04  ,  .  0  0  3 ,  .  0  0  3  0  0  2  5  , .  00  1  6  ,  .  00  1  3 , 

21  .001  3, .  28  16  , 

3 4  OR  NO 

3  6  A  IV  R 

2  6  YK=1O9O0O.  ,2*8000.  , 5*80 00. ,5*7764., 5*17360. ,2*2100., 2*2500. ,1  600., 

20.26  2*  1309  . ,  261  1  0.  , 

27  y<H  =  0.  ,5*4  000.,  12*20  00.  ,7*  200C.  ,20  0  0., 

3  P  C=  249  ,7*  10.  ,  3*16  2*20  .,  5*25  . ,  3*  2. 6  ,4  *  1 .  3,6.  3  4, 

•’ 0  XX  I^O.  ,7*400  0.  ,10*  10  000.  ,8*2400.  , 

c.  C*  =2  0., 7*1  0., 1  0*  20. ,7*1.  ,2.  , 

M  XKS=  19  00.,  7*0.,  6  *0.,  2*0.  ,0.  ,23.  ,8*  0.  , 

3  3  ARNO 

11  r  ^  t  if 

14  n.\C9 -9.,  SCFF 1  =.  01  ,SCH2=  .06 , 5:H3=.  2,  AA=1  C.,  1  0.,  TT1  =  .  05  ,TT2=.  0  5,i:  1  =  0, 

14.2  6  .10  DR  =  1 

36  ARNO 

16  A  A  RV  RM 

17  OR  HI  =1. 2,  OCH  2=1. 4, CM  1=  300.  ,R H 2=900. ,XKSR=? 60., 

1 1  S"  =  7*.  6,  19*1 

3  0  ARNO 

o5'  ?  r  l  r 


INPUT  DATA  FOR  CADAVER  #24  13  HYPEREXTENDED  MODE 


ci 


7  7  4  1  1H 
1 
X 
\ 

4 

A 

A 

7 

ii 

■I 

1 “ 

1  1 
IX 
1  1 

1  V  X  7 
1  V  A 
1  <  .  6 
1  4 
1  7 

I  f. 

17 
1  -1 
1  • 

■>7 

I I 
•>  y 
?  1 
1  4 
?7 

7f, 

>1 
>  » 

7  0 

yi 

■n 

r> 

i  < 

14 
1  s 

17.25 

v. 

17 
!M 
Vi 
u  - 

O9  FILE 


FINE 

8  =4 .2 , 4  .  7 1  19  41  , 4.  908710,4. 970212,  A.  OOR  14  9  ,4.  964956  ,  4.  93  7202,  4.  R1  A  A  HI  , 
4.776061,4.7  42100,4.782  19  1,4.  9  37372,  '5.207074,7  .7  09657  ,5.9  19  2  57, 

6.  160757,6.4  91811,  6.77  0.170  ,6.938817,7.  0  19  748,  7.  1  1  6378,  7 .1730  3/  , 

7.  217  710  ,7.  241  061,7.24  1061,7.241  06  1  , 

I  END 

A”w  n 

W=  1.  6,7.  1614  68,  8.875491 ,10.554  218,  12.  252675,  1  1.  902  007,  1  5  .4  50»13  ,  16.3  9  79 
13. 2  4S  24  6,  19 .506104,  20.624  32  9,  21 .7  31750,22.8  17  215,  21.  806976,2  4  .7  04524  , 
27.564  18  3,26.4  17  39  1,  27.  2  45  28  5,  27.  9  7715  9,2  8  .3  3  7  7  4)  ,  29.  36  341  0,3  0.061  081  , 

3  9.7  59  62  8,  1  1  .4  59  00  0, 3  2.  1  58  997,  34 .453  984, 

7  END 
AT  HP  FF 

0  =-l  0.  ,-5  .  ,  -  1.  5,-3.  ,1.  5,  1 . 5,  1. 5,  3. 5,  1 .  ,0.  , -5.  ,-12.  ,-  17.  ,  -  17.  ,  -22  . , 

-2  1. ,-21  .,  -15.J-9.  ,-7.  ,-5.  ,-4.  ,-  3.  ,1*0., 

AEND 
A  P‘0  HP 

9-2. 16, 4  *.3  5  ,.  6  ,.  55, 2*.  5,  4*.  46,  4  *.375,  .1,7*.  25, 1.  8  , 
p  =  9.,2*0.,6*O.  ,.9, .7  ,.65,.  5,. 2, 4*0.,  7*. 5,  1., 

8=2*  1.  ,3*1.1,  1*1.2,  1*1.1  ,1.05,  .  84  ,  .94, .9  , .9  7  , .9, .  3  , .7  7,.  72, .7, 

1*  .7,  .7,  .7  , 

V 1  =  . 21  47,  .01  25  ,.  0125  ,7*.  00  49  11,5*.  007764, 5*. 01  736, 2*. 002  1  ,2*.  0  0  25, 

.  101  6,  2*. 001  3,  .026  1  1  , 

n  =  .  45  4,  2*.  0  7,  4*.  01  563,  .  06  25,  .05  34  ,.  C54  6  ,.04  9  1,.  0  476,.  0  3  92,.  03  4  0, 

.  0278,  .0255,  .  0 1  84  ,  . 0  06  6  ,  . 0  0  4  ,.  00  3,  .  C03  ,  .  0025,  .0016  ,.00  13  , 

.00  1  J, .  2  816, 

A  END 
r-  9 1  V  E 

XX  =1  0000  0.,  2  *3  00  0.  ,5  *8  00  0.  ,  5*  7  764  .  ,  5  *1  73 6  0 . ,  2  *2 1  00  .  ,2*2 5 00 .  ,  16 0  0  .  , 

2*  110  0.  ,  261  1  0., 

XKH=0.  ,5*40  0  0.  ,1  2*2000.,  7  *20  00.  ,2000.  , 

C  =  24  0,  7*  10.,  3*  15  .  ,2*20  .  ,5*25.  ,  1*2.  6,  4*1  .1,6.14  , 

X  XT-"  0.  ,7  *4000.  ,  1  0*  10000.  ,8*2400.  , 

CT  =  2  0.  ,7  *10.  ,  10*20.  ,  7*1  .  ,2  ., 

XK  5=1000  . ,7*0.  ,6*0.  ,  2*0.  ,  0.,  2  0.,  8*0., 

60  NO 
67  IX 

[>Ac:c  =  9.  ,  SDH  1=.  0  I  ,5CH  2=.  06,  SCH1  =  .2,  A  A  =5 .5  ,5.  ,TT  1=  .  0  1  5,  TT2  =  .  04 ,  I  Cl  =0  , 

'in  d  i-: = 01 

AE  ND 
AS  EVEN 

DC  HI  =  1.  2,03112=  1.  4,CH  1  =3  0().,CH2  =9  00.  ,  XKSR=  59.  , 

7P=7*.7,  19*1.  , 

6E  NO 


APPENDIX  III 


OUTPUT  DATA  FOR  CADAVER  #24  I x 
FOR  INPUT  DATA  SHOWN  ON  P.  |82 


HYPEPEXTENDED 


MODE 


183 


PLOT  DESCRIPTION  GENERATION  BEGINS 

V  =  O.OPn  SEAT  L  0  A  D=  0.0  STRAP=  0.0  LAP=  0.0  HEAP  A  NO  LF=  0.0  KOUNT=  2 


LINK 

A  XT  AL 

SH  EAR 

HOME  NT 

FACET 

NO. 

FORCE  ( L3  ) 

FORCE  (LB) 

(IN- LB) 

FOFCE  (LB) 

1 

-0.0  ) H 0 

0.0 

-0.0 

0.0 

■> 

-0.  OOP'"' 

0.0000 

0. oooo 

0 .0000 

3 

0. 0001 

-0.0000 

-  0.  oooo 

-0.0000 

u 

-  0.  00  0  0 

-0.0000 

0.0000 

-0.0000 

8 

-  0. oooo 

-0.0000 

-0. oooo 

-r  .oooo 

h 

-0.  ■)  J00 

-0.0000 

0. oooo 

f  .0 

7 

-0. oooo 

-0.0000 

-  0.  oooo 

o  .0 

0 

-  0. 00 no 

-0.0000 

0.0000 

-0.0000 

0 

-  J  .  0  0  0  n 

-0.0000 

-0.0000 

P.0 

1  0 

0. 0000 

0.0000 

-  o.  on o (i 

0 . 0 

1 1 

0. oooo 

0.0000 

-0.  oooo 

0 .0 

1  2 

-  ).  oooo 

-0.0000 

-0.  oooo 

0 . 0 

1  1 

- n. OOOO 

-0.0000 

-0.0000 

0.0 

14 

-0.  oooo 

-0.0000 

-0.  oooo 

0.0 

1  s 

-0.  0 0 0n 

-0.0000 

0.  oooo 

0 . 0 

10 

-0.  0  00  ) 

-o.ooco 

o.  oooo 

0.0 

1  7 

-0. oooo 

-0.0000 

0.  oooo 

0.0 

1  H 

-o. oooo 

-0.0000 

0.000  0 

0.0 

1  1 

0.  ) 

0.0 

0.  0 

o  .  0 

■)(' 

0.  0 

0.0 

0.0 

0.0 

?  1 

0.0 

0.0 

0.0 

0.0 

7  2 

0.  0 

0.0 

0.0 

0 . 0 

2  1 

0.  0 

0.0 

0.  0 

O.o 

24 

0.  ) 

0.0 

0.0 

0.0 

25 

0.  0 

0.0 

0.  0 
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-3.1173 

143 . 3067 

2  0 

-0.2698 

-4.  47.28 

1 37.6878 

21 

-0.2  272 

-5.495  1 

118.2677 

22 

-0 . 2199 

-6.1770 

8  3.0117 

2  3 

-0. 2009 

-6.  4  7  72 

26.8)74 

24 

-0 . 1 454 

-6.  4 1 44 

-47  .861  1 

2  S 

-0.1 408 

-6 .0.24  4 

-138.3677 

26 

-0.  1  )28 

-5.6622 

-211.9927 

r-n. 

09  5  SEAT  LOAD 

=  706.4  STR  AP  =  14.8  LAP  =  0.0  HEAD  ANSLE=-0 

.164  KOUNT' 

LINK 

A  XT  AI. 

SHEAR 

HOME  NT 

FACET 

NO. 

P  0 P 1 : t  (  L  2) 

^ORCE  (LB) 

(TN-  LI) ) 

FORCE  ( Lli ) 

k 


1 

-70b.  47H  7 

0.0 

-0.0 

0.0 

; 

-  1  82.  784  1 

9.1868 

-17.2721 

-1  1  2.992  0 

3 

-  16b.  If'  14 

1 .4693 

-Id.  96  7  8 

-  1  06.2  363 

a 

-  ISO.  2  10  2 

-7  .4  95  0 

-  7.  68  36 

-90. 0690 

S 

-  16  1.  92  4  2 

-9.3265 

3. 7058 

-78.8  17  7 

6 

-171.  4  0  7  S 

-  22.990  2 

2  3.  79  04 

-59.3833 

7 

-  17  2. 27  0  } 

-28.8  34  2 

3  1.  5220 

-30.1 887 

■1 

-  19  9.  4n  16 

-48.9545 

62. 226 1 

-  1  3.  7533 

9 

-Hi.  12  9.7 

-45.8  38  0 

122.  456  3 

-22 .5050 

r 

-  19  1.  OH  29 

-39 .0047 

153.7979 

-  14.  7642 

1 1 

-  17  6.  0626 

-  37.  64  70 

170.9114 

-  10.7  112 

i  / 

-167.4661 

-  25.  3210 

175. 9962 

-  10 .0872 

1  5 

-IS J. 6479 

-5.7231 

16  5.  12  6  2 

-1U.HH3 

14 

-  16  1.  0  144 

13.5493 

12  1.  00  09 

3.  4  87  4 

1  5 

-  1  1  4.  5  2  29 

14.  13  5  9 

9  8 . 3  6  8  3 

C  .898  7 

1  6 

-  1  OH.  6004 

2d  .9048 

69.  59  12 

0.57  3  C 

1  / 

-n  3.  29  4  6 

22  .6056 

45.9 24b 

-1.7257 

19 

-6  8.  79  91 

20.0682 

24. 64  14 

-2. 9674 

1  9 

-  3  2.  34  26 

10.0437 

46.6062 

-1.8552 

.7  0 

-26.6919 

9.090  4 

3  4.4763 

-5.4  06  3 

21 

-2  2.  81  36 

9.0362 

21  .  6  3  3  6 

-6 .392  1 

2  7 

-  1  7.  6106 

8.7255 

10.6129 

-8.  4  56  6 

>.  3 

-10.  36  16 

6.6239 

2  .  36  92 

-  12.  366  9 

>n 

-6. 6900 

6.060  4 

-6.  3897 

-14.2880 

is 

-4 .  4  8 H  o 

6.12  34 

-12.  5656 

-  15.  214  1 

26 

-  6.  60  6  2 

39.6656 

-21.021b 

-  12.9  26  7 

LT  W 

0 

unr 

11  DO  T 

W 

WOT 

W  DOT 

NO. 

(I  N  8) 

(  IN  /SEC) 

(IN/5  EC  2) 

(T  NS) 

( I M  /S  EC) 

(  IN/SEC  2) 

1 

4.  2  0  00 

0.0 

0.0 

3.  592  8 

0.1081 

-5.  3000 

2 

4.  7  77  i 

0.4280 

-6.5479 

7.1310 

0.  7  0  35 

-52.1463 

1 

4.89  19 

2.  0  38  1 

-1 6.4648 

8.8  04  2 

1.  3107 

-97.611  6 

4 

4.9  3  46 

4  .  3  59  4 

-29. 2726 

1  C  .4  8  3  2 

2.0258 

-  149.  6363 

5 

4.  9612 

7.093 1 

-26.  12  1  3 

12.  1  61  2 

2.8169 

-295.  16  52 

6 

4 .8403 

3.  9  48  3 

2. 8846 

1  '  .788  1 

3.  80  52 

-265.1634 

7 

4.92  16 

12.  6  69  7 

6  9.  83  2  3 

15.3118 

4.  9  b40 

-327. 6741 

9 

4 .7416 

14.9817 

17  1.  9  2  39 

16.7339 

6.20  31 

-  336.  7  685 

9 

4.  70  94 

16.6311 

293. 5880 

18.0  58  4 

7.  3695 

-44  0.3  9  18 

1  0 

4  .  7  36  1 

17.  6  50  3 

4  06. 0842 

1  9  .29  3  6 

8.  470  5 

-  493.  651  5 

1  1 

4.8107 

17.9  968 

4  9  0.619'. 

20. 384  4 

9.  5427 

-  545.6042 

12 

6.  06  60 

17.8  119 

548.  108  3 

2  1  .45  5  8 

10.6290 

-  599.  6073 

1  1 

6.  4  1  64 

17.3  032 

5/2.  6  134 

2  2.  4  9  30 

1  1 . 7770 

-650. 4825 

14 

6.  9  16  3 

16.6  367 

56  8.61  51 

2  3 . 4  37  1 

12. 4948 

-669.3779 

1  S 

6 . 2  2  3  1 

16.  1  42  5 

5  3  3. 01  5  3 

24.2856 

1  3.  0859 

-669. 8b80 

16 

6> .  6  6  74 

15.  3070 

4  6  7.  62  97 

25.0952 

1  3. 50  1  1 

-649.  4292 

1  / 

/.  1  0  87 

15.9  85  3 

37  3.  18  42 

25.  908  2 

13.6758 

-610.6764 

1  f; 

7.49  92 

If..  3  92  5 

26  1.4288 

26.668  8 

1  3.  6  355 

-565.1453 

1  9 

7.  76  67 

17.  4  20  3 

1 09. 3485 

27.3  52  1 

1  3.  2  50  7 

-526. 7672 

7  0 

7.9660 

19. 406  1 

-  a  4.  114  3 

28.0070 

12.6942 

-  50  5.  4  52  2 

2  1 

H.  1  6  93 

22.2692 

-241.5851 

28. 6 698 

1 1  .  9  HO  8 

-499. 0470 

7  2 

1.  3  193 

2b.  8  64  9 

-406.0461 

29  .  3  38  9 

1  1.  1  671 

-514.1 789 

2  1 

8.6  103 

30.  116  9 

-547.  1592 

30.010 8 

10.  3477 

-564. 5419 

7  4 

3.6610 

34.85  19 

-6  38. 0  306 

3  0.6  88  4 

9.666  3 

-648.  7058 

26 

8.  7  3  72 

39. 844  1 

-658.425b 

3  1.  3  72  9 

9.2053 

-7  58.  .26  57 

6  6 

9.  16  8  2 

55.  2  309 

-21 8.7602 

3  3.6  40  9 

6.  6728 

-947.1 032 

I.  T  N  K 

0 

0  OT 

ODDT 

NO. 

(PAD) 

(  RAD/C  EC) 

(RAD/5  EC  2) 

1 

-3. 1  74b 

0.0 

0.0 

') 

-0.1 30  1 

-0.7644 

2  .9  865 

212 


) 

-'1.0  17  2 

-  1.  266b 

7.1  106 

4 

J  4  14 

-1  .8417 

-0 . 3400 

S 

0.1  1ft 2 

-1 .6042 

-13.8272 

6 

0.0  10  J 

-  1. 6580 

-37.521 1 

7 

0 .0573 

-1.6267 

-  60. 3690 

8 

0.0  416 

-1.2460 

-86.5697 

1 

-0  .  o  1 4  2 

-0.8621 

-90.0408 

1  0 

-0.0  471 

-0.3678 

' -82 .6209 

1 1 

-0.1  jl  7 

0. 1380 

-  64.9792 

1  2 

-0. 201ft 

0.8509 

-41.5355 

n 

-0.  10  b4 

0.7947 

-  1  3 .5765 

14 

-0.4  07b 

0.7315 

20.6107 

ib 

-0.  ft  04  b 

0.499  1 

57.4680 

ift 

-0. 4040 

0.1433 

91.5702 

17 

-0.4 00  ft 

-0.2836 

1  25 .964  0 

18 

-('.  306  3 

-0.7299 

150.940? 

18 

-0.  1  in  j 

-2.2587 

1  92  .5  370 

20 

-0.2  402 

-  1. 5900 

205. 4280 

21 

-0.2620 

-4.702  1 

190.0000 

22 

-0.240  b 

-5.5901 

147. 3660 

>  3 

-0.  2  12  7 

-6.221  1 

76.0572 

>  '! 

-0 . 1 780 

-ft.  5893 

-18.4360 

2b 

-0 . 172  7 

-6.7122 

-  1  32.0  51  4 

2ft 

-0 . 1630 

-6.7787 

-2  30.0  39  2 

T-8. 

100  SPAT  LOAD 

-  6  6  3.6  STR  AP  =  16.4  L 

AP  -  0.0  HE  A  0  A  NS  LE  =-0.2  01  KOUNT 

LINK 

AX7  AL 

SHEAR 

(10  ME  NT 

FACET 

NO. 

FORCE  (1.3) 

FORCE  (LB) 

(TN-  LB) 

force  (I.B) 

1 

-  6ft  3.  62  04 

0.0 

-0.0 

0.0 

2 

-  16  4.  700  4 

5.  892  9 

-2. 81 46 

-87 .7  25  1 

J 

-14b.  722b 

0.0281 

-7.0877 

-83.4616 

4 

-11b.  162  1 

-7 .9956 

-  1.  4695 

-70.7246 

b 

-  1  3ft.  68  20 

-9.4346 

7. 6454 

-61.0146 

ft 

-  14  3.  0688 

-  20.6066 

2  4.  34  54 

-44 . 28  n 

7 

- IbO. 068 7 

-  24.  4  2  7  6 

31.  8581 

-21 .2829 

2 

-  1 6  > .  6o  34 

-  30.47  18 

5  5.  92  06 

-8.2109 

9 

-  14  7.  Ob 04 

-36.887  1 

10  4.  29  89 

-16.8310 

10 

-  14b.  1112 

-  32.  152  8 

12  9.2898 

-  10.  18  26 

1  1 

-14  0.  10  1b 

-32.6  6  18 

144.  5455 

-l.  .76 2  3 

12 

-  114.  4  IbO 

-2 1 .881  0 

15  2.  b  2ft  ft 

-  8.  7  08  C 

1  1 

-  124.  1  1J  ft 

-  8.  987  5 

14  0 .801 11 

-ft.  2  17  3 

14 

-131.  701) 

ft  .  19  7  3 

119.  C4  4  ft 

(.  .6  80  1 

1  b 

-  1 1  0.  SO  72 

/ .  19  0  8 

10  4.  2  ft  ft  5 

4  .  8322 

Ift 

-^0,  66  64 

1 3.8972 

8  1.  8  ft  40 

b.  1  86  9 

1  / 

- ft  0  .  7  '4  1 

1  7.962  3 

ft? .5616 

2.9  4  40 

1  8 

-4  0.  100  8 

18.  44  3  1 

43 . 5867 

1.  3334 

10 

-3  2.  »S  1  1 

10  .  561  1 

63.36 36 

8.6116 

H) 

-26.  7812 

8  .480  1 

49.  8772 

4.0245 

2  1 

-  2  2.  6  17  3 

9.8398 

3  5.  30  2b 

2.  0  95  6 

22 

-1  7.011ft 

0.8864 

22. 64  86 

-0.0436 

2  1 

-0.802 7 

7.0224 

12.  38  80 

-5.5392 

2  4 

- b . 70 ft  8 

7.238  4 

2.  30  0  1 

-8 . 0442 

2b 

-  3.  10  87 

7.1 705 

-6.  94  38 

-9. ft  04 7 

2ft 

-  1.  7  4  2ft 

38.  42  ^2 

-15.9872 

-in. 2638 

LINK 

U 

UDT 

UDO  T 

W  WOT  WDD  !’ 

NO. 

(I  NS) 

(IN/SEC) 

(IN/SEC2) 

(T  NS)  (TM/5EC)  (TN/SEC2) 

1 

4. 2 000 

0.0 

0.  0 

3.  50  3  3  0.  0764  -7.2  848 

2 

4  .  7  7  0  1 

0 .  3  416 

-25.2602 

7.1  3  19  0.  4276  -58.  4702 

1 

4.4014 

1 .7400 

-88.  7871 

8.8095  0.77  17  -  1  17.  4167 

213 


4 

4.  85  54 

.3. 8  94  4 

-  14  1. 66  7  1 

10.  4  91  3 

1 .  1725 

) 

4.  8  8  88 

6.6010 

-  1  55.  42  1  1 

12.1725 

1.  6  358 

6 

4 . 8  3  81 

9.5  88  5 

-  1  3  3.  74  1  9 

13.8015 

2.280/ 

/ 

4.  88  82 

12.6156 

-8  1. 09  36 

15.  332  1 

3.0805 

8 

4. 8 179 

15.3584 

-  1  2.79  34 

16.7 59  5 

3.  972  9 

Q 

4.7842 

17.  531  1 

7  1.6313 

18.0891 

4.  8080 

10 

4.82  84 

1  9 . 0  86  5 

169.  7  326 

19.3291 

5.5786 

1  1 

4.  8  2  88 

1 9.912  3 

273.2322 

20.  4  24  5 

6.3234 

1  2 

5  .  15  14 

20.  167  3 

379. 9841 

21 .5004 

7.  0  627 

1  3 

5.5097 

19.9  30  3 

470.  0  7  14 

22  .542  6 

7  .  B  6  7  3 

14 

5.  9  0  64 

1  9.3  87  3 

5  2  2.  12  37 

23  .4  9CC 

8.4218 

15 

ii.  310  5 

18.8099 

524 . 83  64 

24.  34  14 

8.  9  6  3b 

1  6 

6  .  7  ->  29 

18.  29  66 

481 . 9015 

2  5  .1  5  3  2 

9.  4547 

17 

7.1  934 

17  .  9  20  6 

39  9.  581  8 

25.9676 

9.8151 

18 

/.  5748 

17.7712 

29  4.  0  5  5  4 

26.7285 

1  0.0  04  6 

18 

/.  84  5  3 

18.0382 

145.7069 

27.  4  105 

9. 9401 

2  n 

8.0613 

19.  1  4  02 

-  3  0.74  65 

28.0  6  29 

9.  4  3  3  6 

21 

3.  26  7  3 

21. 0809 

-.18.6522 

28.7  224 

8.  8 p  8  H 

22 

8. 4  6  12 

23.7958 

-  40  2.  5  3  31 

29  . ) 87  3 

7.9864 

2  1 

1.  65  35 

27.2  53  2 

-  57  4  .  5  3  20 

30.  0  54  5 

6.9890 

2  4 

8.8  >68 

31.  421  7 

-706.0094 

3  0 . 7  27  8 

5.  9448 

2  5 

3.97  75 

36.  195  2 

-  770.  79  1  3 

31  .4  C87 

4.960  1 

36 

9.  44  C'J 

53 .6633 

-  39  0.  8  4  28 

33.6  61  5 

1  .  3  805 

L  IN 

0 

ODT 

0  HOT 

NJO  . 

(RAO) 

(RAD/SEC) 

(R  AD/S  EC2 ) 

1 

-0.1  74  5 

0  .0 

0.0 

2 

-0.  10  38 

-0.6707 

39  .  229  5 

3 

-0.0  233 

-1.1605 

31  .974  3 

4 

-0.0491 

-1.5113 

11  .  71  0  0 

5 

0.0276 

-1  .753  1 

-8.7413 

b 

0.0  218 

-1  .83  55 

-  3  1 .  3  20  7 

7 

0.0  48  5 

-1  .884  8 

-41.3511 

8 

0 . 0  3  4  4 

-1.6112 

-57.4348 

) 

-0.  )  197 

-  '.  2762 

-72.934  3 

1  0 

-0.0502 

-0.8  09  1 

-90. 1 455 

1  1 

-0 . 1 520 

-0.2860 

-10  1  .0399 

1  > 

-0 . 2  84  5 

0.2049 

-  97  .5  9  7  3 

1  i 

-0 . J  318 

0 . 5723 

-72.9425 

1  4 

-n .4 8 38 

0.7113 

-27.8 930 

15 

-0. 5  n  4 

0.7087 

24  .  28  1  1 

16 

-0.4922 

0.5732 

73.3809 

1  / 

-0.4  9b 4 

0.3334 

1  1  4.  1 67  1 

18 

- 0  .  3  9  8  0 

0.0230 

1  4  3.5  2  36 

1  5 

-3) .  i  ll>1 

-1.2717 

1  95 .499  5 

l  n 

-■).  3  04  5 

-2.5136 

219.3491 

31 

-0.2  13 '1 

-  3  .677  2 

214.  3833 

2  2 

-8.2/54 

-4.7552 

1 8  1 . 2  50  6 

>  1 

-0 .26 27 

-5.7306 

116.1371 

24 

-0.2  110 

-6.579  3 

21  .1  02  0 

28 

-0.2  07  8 

-7.3006 

-  100 .8  480 

26 

-0. 2006 

-7.9069 

-2  14.7  00  1 

T  -  0 . 1 

05  s A T  1.0  A  0 

=  6  39.  7  0  T  RA  V=  18.2  L  AP  = 

0.0  MEAD 

A  NO  LE  =-0 

L  T  NK 

AXT  AL 

811  EAR 

MO  ME  NE 

"ACET 

no. 

-  on  OF  (Li!) 

FORCE  (LB) 

( IN  -  L  8  ) 

F  OR  C  p  (LB) 

1 

-  6  19.  7  24  2 

0.0 

-0.0 

0 

.0 

; 

-157.544 1 

3.  0  58  1 

8.  e4  51 

-7  0 

.645  9 

I 

-  1  3  7.  7457 

-1 .8808 

2.  514  3 

-67 

.1  C8C 

4 

-  1  zb.  527  4 

- 9  .4  854 

6 . 56 04 

-  56 

.  2456 

-  1  8 9.  08  8  8 
-282.  6296 
-118.  1 8  2  2 
-4  17.  >089 
-485.8828 
-572.  41  J7 
-648.  1043 
-726.5422 
-808. 1  785 
-88  2.  3182 
-9  17.  / U  2 5 

-  85 o.  8 885 

-  84 H. 1 6  85 

-  8  1  5.  5671 
-87  2  .  1  1  35 
-82b. 4721 
-7  11  .  0855 
-764.81  3  1 
-755.  5  172 
-777.  /  1b 5 
-8  3  8. 51  o 3 
-8  36.  5604 

-  1  16  4.  17  40 


.  2  4  3  K 0 U N T= 


2 1  4 


) 

-  12b.  4  SO  2 

-4.3274 

11  .546  4 

-47.8571 

(, 

- 1 1 0.  / SHI 

-14.1722 

27.  14  20 

-  3  3 . 4  7  82 

7 

-  14  0.  46  4  b 

-20.46  1  8 

1 1.  0854 

-15.5243 

3 

-  14  1.  01  66 

-32.3088 

5 1 .  38  16 

-5.  1274 

) 

-  1  2  7  .  4  36  S 

-  24.  7774 

90.  0444 

-  14.0055 

1 ') 

-  12  3.  7 >3 7  S 

-  26.  4009 

108. 9543 

-8.6325 

1  1 

-1113.  IMS') 

-27.4  2  1  1 

120. 6R62 

-5.6270 

1 

-113.  2242 

-20.982 1 

12  8.  3  384 

-  4.  4855 

1  1 

-  V'4.  7  20  4 

-4. 4H8S 

1  2  8.  62  38 

-7.0997 

1  ft 

-  1  1  1  .  3  3  2  7 

2.  1740 

106.  1741 

6  .  8  6  i  5 

15 

-44.2312 

1  .878  1 

97.  6  3  12 

8.774  1 

IS 

-7  9.  2M(3  3 
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